[ACM]CCF CSP [201509-5]E题 最佳文章【90分】

思路：AC自动机把所有可能的后缀处理出来(最多有100个不同后缀)，然后dp，F[i][u]表示i个字符后缀为AC自动机里的第u个后缀，最多匹配次数。

max(f[m][i],i=0..100)即答案，这样可以拿到【60分】，因为m可以非常大。

但是我们可以容易证明：当m趋于无限大时，最优情况下，一定出现重复串，而且是重复串连续放在一起。

这里，我们假设当m=1000时，一定会产生重复串，把字符序列输出，找出起重复串(并记录这串的匹配数量)，然后从m中减去重复串，最后m留下不超过1000的。

然后都是瞎搞【90分】 还有10分不知道为什么错，也许1000取太小？

#include
#include
#include
#include
#include
using namespace std;
#define ACMXNODE 110   //AC自动机节点数(最大为所有字符串字符个数和)
int charIdx(char ch){ return ch-'a'; }
#define CHARSIZE 30   //不同字符个数(128可以处理大部分情况，空间紧张可以减少这个值,但要保证charIdx每个字符都对应一个下标)
/*
  用法：
  ac.clearMap();
  ac.ins(str,val);//相同字符串会相互覆盖
  ac.getfail()
*/
  int cnt[ACMXNODE];       //cnt[u]：插入的字符串中，等于字符串u的某个后缀的字符串个数
class ACAUTO{
public:
  int sz;
  int ch[ACMXNODE][CHARSIZE];
  int val[ACMXNODE];       //节点标记值，0表示非叶子节点。叶子节点记录一个自定义值（如是第几个字符串）。（相同字符串会互相覆盖）
  int f[ACMXNODE];         //fail数组
  int last[ACMXNODE];      last[u]：上一个u结尾的字符串。如有两个字符串(abcd,bc)，节点u表示abc,则节点last[u]指向bc。
  int print(int i,int j);
public:
  ACAUTO();
  //初始化+建图
  int clearMap();          //初始化
  int ins(char *s,int v);  //插入一个字符串，v代表字符串结尾节点标记值。(0：非字符串结尾，其他：自定义)
  int getFail();           //建立失配指针。bfs查找
  int kmpFind(char *T);    //找出trie树中所有为T子串的字符串
};
ACAUTO::ACAUTO(){clearMap();}
int ACAUTO::clearMap(){
  sz=1;
  memset(ch[0],0,sizeof(ch[0]));
  memset(val,0,sizeof(val));
  memset(cnt,0,sizeof(cnt));
  return 0;
}
int ACAUTO::ins(char* s,int v=1){
  int u=0;
  int n=strlen(s);
  for (int i=0;i q;
  //先处理节点0的fail值
  f[0]=0;
  for (int c=0;c0);
    }
  }
}

int ACAUTO::kmpFind(char *T){
  int n=strlen(T);
  int j=0;
  for (int i=0;i path[101000][110];
vector g[110];
int nodeList[1100],nodeListSz=0;
int repList[1100],repListSz=0,repGain=0;
void outputPath(int i,int u){
  if (i==0)return ;
  outputPath(i-1,path[i][u].first);
  nodeList[++nodeListSz]=u;
  //putchar(path[i][u].second+'a');
}
int main(){
  int n;long long m;
  scanf("%d",&n);cin>>m;
  char s[110];
  ac.clearMap();
  for (int i=1;i<=n;i++){
    scanf("%s",s);
    ac.ins(s,i);
  }
  ac.getFail();
  //
  for (int u=0;uv (c)
      g[v].push_back(u);
    }
  }
  memset(f,255,sizeof(f));
  memset(path,0,sizeof(path));
  f[0][0]=0;
  int sz=ac.sz;
  int rep=0;
  #define FOR(i,s,t) for (int i=s;i<=t;i++)
  for (int i=1;i<=1000;i++){
    f[i][sz]=0;
    for (int u=0;uu
         if (f[i-1][v]==-1)continue;
         if (f[i-1][v]>f[i][u]){
            f[i][u]=f[i-1][v];
            path[i][u].first=v;
         }
       }
       //
      int v=path[i][u].first;
      FOR(c,0,CHARSIZE-1)if (ac.ch[v][c]==u){path[i][u].second=c;break;}
      //
      if (f[i][u]!=-1) f[i][u]+=cnt[u];
      if (f[i][sz]\n",nodeList[i]);
    repList[++repListSz]=nodeList[i];
    if (nodeList[i]==nodeList[399]) { repGain=f[i][sz]-f[399][sz];break;}
  }
  //FOR(i,1,repListSz) printf("%d ",repList[i]);putchar('\n');
  //printf("RepGain=%d\n",repGain);
  long long p=max(m/repListSz-6,0ll);//p段的重复串提取出来
  m-=p*repListSz;
  cout<