注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

一个蒟蒻的代码回收站

最后一次省选求rp

 
 
 

日志

 
 

bzoj1195【HNOI2006】最短母串解题报告  

2014-04-17 11:05:00|  分类: bzoj |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

题目链接

我们先把这么大堆的串建AC自动机

然后在AC自动机上跑dp,由于要求的不是长度,而是串,所以要附加一点东西

dp[i][j]表示跑到自动机上节点i,匹配状态为j的最短长度(j为二进制串表示状态)。

f[i][j]表示dp[i][j]是从哪里转移过来的(f[i][j]应该是一个pair<int,int>的类型,由于我不太会用pair,所以自己写了个struct)

然后为了保证字典序最小,我们从A到Z开始加字母,然后就没有然后了……

这样在bzoj上就可以AC了,但是本地测似乎是会TLE的……

我们可以把每个节点放字母可以走到哪个位置预处理出来,然后就不用每次跑fail了……

优化挺明显的:3712ms->856ms

贴代码

#include<cstdio>
#include<cstring>
#include<string>
#include<queue>
#include<iostream>
using namespace std;

#define maxn 620
#define maxm 26
const char base='A';
int root,tot,n,m;
struct trie
{
int num;
int fail;
char c;
int ch[maxm];
//bool isend;
} t[maxn];

void clear()
{
root=tot=2;
}
void add(char *s,int x)
{
int n=strlen(s);
int q=root;
for (int i=0;i<n;i++)
{
int tmp=s[i]-base;
if (!t[q].ch[tmp]) t[q].ch[tmp]=++tot;
q=t[q].ch[tmp];
t[q].c=s[i];
if (i==n-1) t[q].num|=x;
}
}
void get_fail()
{
queue<int> Q;
t[root].fail=root;
for (int i=0;i<maxm;i++)
{
int u=t[root].ch[i];
if (u) {t[u].fail=root;Q.push(u);} else t[root].ch[i]=root;
}
while (!Q.empty())
{
int now=Q.front();
Q.pop();
for (int i=0;i<maxm;i++)
{
int u=t[now].ch[i];
if (!u) {t[now].ch[i]=t[t[now].fail].ch[i];continue;}
int v=t[now].fail;
while (v!=root && !t[v].ch[i]) v=t[v].fail;
t[u].fail=t[v].ch[i] ? t[v].ch[i] : root;
Q.push(u);
}
}
}
void train_fail()
{
for (int i=root+1;i<=tot;i++)
{
int j=i;
while (j!=root)
{
j=t[j].fail;
t[i].num|=t[j].num;
}
}
}
struct node
{
int first,second;
} f[maxn][1<<12];
char ch[50];
int dp[maxn][1<<12];
queue<node> Q;
string s,s1;
int ans;
int main()
{
scanf("%d",&n);
clear();
for (int i=0;i<n;i++)
{
scanf("%s",&ch);
add(ch,1<<i);
}
get_fail();
train_fail();
dp[root][0]=0;
Q.push((node){root,0});
while (!Q.empty())
{
node now;
now=Q.front();
for (int i=0;i<26;i++)
{
int u=now.first;
while (u!=root && !t[u].ch[i]) u=t[u].fail;
u=t[u].ch[i] ? t[u].ch[i] : root;
int q=now.second|t[u].num;
if (dp[u][q]>dp[now.first][now.second]+1 || dp[u][q]==0 && (u!=root || q!=0))
dp[u][q]=dp[now.first][now.second]+1,f[u][q]=now,Q.push((node){u,q});
}
Q.pop();
}
ans=1<<30;
for (int i=root;i<=tot;i++)
if (t[i].num!=0)
if (dp[i][(1<<n)-1]<=ans && dp[i][(1<<n)-1]!=0)
{
s1="";
int tmp=(1<<n)-1;
int j=i;
while (j!=root)
{
s1=t[j].c+s1;
int q=j;
j=f[q][tmp].first;
tmp=f[q][tmp].second;
}
if (ans>dp[i][(1<<n)-1]) {ans=dp[i][(1<<n)-1];s=s1;continue;}
if (s1<s) s=s1;
}
cout << s <<endl;
}


  评论这张
 
阅读(7)| 评论(0)
推荐 转载

历史上的今天

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2018