[HDU3247]Resource Archiver

时间:2022-09-18 15:26:52

  AC自动机+状压DP

  首先对所有串建AC自动机,然后对于每个资源串,算出从串末走到其他资源串末所需的距离(中途避开非法点)

  也就是算出两两间的距离。。。然后就变成旅行商问题了。

  计算距离的时候要考虑一下。。直接从末尾走的话可能求出来的并不是最短路

  所以要从fail和ch(就是虚边和实边)一起转移,走实边的话加进队列,虚边只能更新距离(不然可能非法)

  但似乎正确性存疑?。。可能存在 三个资源串,分别按照最短距离拼在一起后会包含病毒串 的情况。。。

  反正我的程序会炸,但网上一些标程没事= =。。。另外本题数据极弱,随机数据就能艹翻一大坨标程系列

 #include<cstdio>
#include<iostream>
#include<cstring>
#define next NEXT
using namespace std;
const int maxn=,inf=;
int ch[maxn][],next[maxn][],fail[maxn],tot;
bool gg[maxn],uu[maxn];
int f1[maxn];
char s[maxn],mp[][];int len[],ed[];
int d[][],dis[maxn],d1[maxn];
int f[][],two[];
int dl[maxn],l,r;
int i,j,k,n,m; inline void getf1(int x){
int i,p=;f1[]=;
for(i=;i<=len[x];i++){
while(p&&mp[x][i]!=mp[x][p+])p=f1[p];
p+=mp[x][i]==mp[x][p+],f1[i]=p;
}
}
inline bool check(int x,int y){
int i,p=;
for(i=;i<=len[y];i++){
while(p&&mp[y][i]!=mp[x][p+])p=f1[p];
p+=mp[y][i]==mp[x][p+];
if(p==len[x])return ;
}
return ;
} inline void trie(int n,int id){
int i,p=;
for(i=;i<=n;i++){
s[i]-='';
if(!ch[p][s[i]])ch[p][s[i]]=++tot,p=tot;
else p=ch[p][s[i]];
}
if(!id)gg[p]=;else ed[id]=p;//printf(" gg:%d\n",p);
}
inline void getfail(){
int i,now,p,j;l=,r=,dl[]=;
while(l<r){
now=dl[++l];
for(i=;i<;/*printf("%d %d next:%d\n",now,i,next[now][i]),*/i++)if(ch[now][i]){
dl[++r]=j=next[now][i]=ch[now][i];
for(p=fail[now];p&&!ch[p][i];p=fail[p]);
if(!now)fail[j]=;else fail[j]=ch[p][i];
gg[j]|=gg[fail[j]];
}else{
for(p=fail[now];p&&!ch[p][i];p=fail[p]);
next[now][i]=ch[p][i];
}
}
for(i=;i<=tot;i++)if(gg[i])
for(j=/*,printf("!! gg:%d\n",i)*/;j<;j++)if(ch[i][j])gg[ch[i][j]]=;
}
inline void bfs(int s){
memset(dis,,(tot+)<<);memset(uu,,tot+);
int l=,r=,i,now,j;
if(!gg[s])dl[++r]=s,dis[s]=d1[s]=,uu[s]=;
for(i=fail[s];i;i=fail[i])if(!gg[i])dis[i]=;
while(l<r){
now=dl[++l];
for(i=;i<;i++)if(!gg[next[now][i]]&&!uu[next[now][i]]){
dl[++r]=next[now][i];uu[dl[r]]=;d1[dl[r]]=d1[now]+;
if(dis[dl[r]]==-)dis[dl[r]]=d1[dl[r]];
for(j=fail[dl[r]];j;j=fail[j])if(dis[j]==-)dis[j]=dis[dl[r]];
// if(s==16)printf("walk:%d->%d\n",now,next[now][i]);
}
}
// if(s==16)
// for(i=0;i<=tot;i++)printf(" %d->%d %d\n",s,i,dis[i]);
}
int main(){
freopen("in.txt","r",stdin);
for(i=two[]=;i<=;i++)two[i]=two[i-]<<;
for(scanf("%d%d",&n,&m);n&&m;scanf("%d%d",&n,&m)){
memset(ch,,(tot+)<<),memset(next,,(tot+)<<),memset(fail,,(tot+)<<),
memset(gg,,tot+),tot=; for(i=;i<=n;i++)
scanf("%s",mp[i]+),mp[i][]='%',len[i]=strlen(mp[i])-;
for(i=;i<n;i++)for(j=i+;j<=n;j++)if(len[i]>len[j])
swap(len[i],len[j]),swap(mp[i],mp[j]); for(i=;i<n;i++){
getf1(i);
for(j=i+;j<=n;j++)if(check(i,j)){gg[i]=;/*printf("%d in %d\n",i,j);*/break;}
}
int n1=;
for(i=;i<=n;i++)if(!gg[i])n1++,memcpy(mp[n1],mp[i],len[i]+),len[n1]=len[i];//,printf("%s\n",mp[n1]);
memset(gg,,n+); for(i=;i<=n1;i++)memcpy(s,mp[i],len[i]+),trie(len[i],i);
for(i=;i<=m;i++)scanf("%s",s+),s[]='%',trie(strlen(s)-,);
getfail(); for(i=;i<=n1;puts(""),i++){
bfs(ed[i]);
for(j=;j<=n1;/*printf("%s-->%s %d\n",mp[i]+1,mp[j]+1,d[i][j]),*/j++)if(dis[ed[j]]!=-)d[i][j]=dis[ed[j]];else d[i][j]=inf;
} int mx=<<n1;
for(i=;i<mx;i++)memset(f[i],,(n1+)<<);
for(i=;i<=n1;i++)f[two[i-]][i]=len[i];
for(i=;i<mx;i++){
for(j=;j<=n1;j++)if(i&two[j-])
for(k=;k<=n1;k++)if(k!=j&&(i&two[k-])&&d[k][j]<inf)
f[i][j]=min(f[i][j],f[i^two[j-]][k]+d[k][j]);
}
int ans=inf;
for(i=;i<=n1;i++)ans=min(ans,f[mx-][i]);
}
return ;
}