标签：
第一行两个整数N，M表示待检查的作文数量，和小强的标准作文库
的行数
接下来M行的01串，表示标准作文库
接下来N行的01串，表示N篇作文
N行，每行一个整数，表示这篇作文的Lo 值。
输入文件不超过1100000字节
注意：题目有改动,可识别的长度不小于90%即可，而不是大于90%
先将作文库建后缀自动机，多篇文章可以通过在中间加入分隔符来完成。对于每组询问，预处理出每一个位置向前最多匹配多长g[]。
二分答案len，dp[]表示匹配到当前位置的最多匹配数，对于i位置，dp[i]由一段通过len，与g确定出的区间[l,r]转移，本来用了一个线段树维护，但是由于时间复杂度O(n*log^2n)，TLE了，观察发现[l,r]是单调的，故可直接用单调队列。
网上一半题解过不了数据：1 1 1 1
另外，对于0.9的问题确实说明了以后能用int就不要用double，实在不行要加eps。
#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> using namespace std; #define MAXN 4100000 #define MAXT MAXN*4 #define INF 0x3f3f3f3f #define lch (now<<1) #define rch (now<<1^1) #define smid ((l+r)>>1) char buf[MAXN]; char *bufnow(buf); int len[MAXN]; char *str[MAXN]; char buf2[MAXN]; char *bufnow2(buf2); struct sam_node { int nxt[3]; int pnt,len; void Print() { for (int i=0;i<3;i++) printf("%d[%d] ",i,nxt[i]); printf("\n"); printf("Pnt:%d\n",pnt); } }sam[MAXN]; int tops=1; int last=1; void Add_item(int w) { int p=last; int np=++tops; sam[np].len=sam[p].len+1; while (p && !sam[p].nxt[w]) sam[p].nxt[w]=np,p=sam[p].pnt; if (!p) { last=np; sam[np].pnt=1; }else { int q=sam[p].nxt[w]; if (sam[p].len+1==sam[q].len) { sam[np].pnt=q; }else { int nq=++tops; sam[nq]=sam[q]; sam[nq].len=sam[p].len+1; sam[nq].pnt=sam[q].pnt; sam[q].pnt=nq; sam[np].pnt=nq; while (p && sam[p].nxt[w]==q) { sam[p].nxt[w]=nq; p=sam[p].pnt; } } } last=np; } int g[MAXN]; int dp[MAXN]; int seq[MAXN]; int main() { freopen("input.txt","r",stdin); int n,m; int x,y; scanf("%d%d\n",&n,&m); for (int i=0;i<m;i++) { scanf("%s\n",bufnow2); bufnow2+=strlen(bufnow2); *(bufnow2++)=‘2‘; } for (int i=0;i<n;i++) { scanf("%s\n",bufnow); str[i]=bufnow; bufnow+=len[i]=strlen(bufnow); bufnow++; } for (char *i=buf2;i!=bufnow2;i++) { Add_item(*i-‘0‘); } for (int i=1;i<=tops;i++) { // printf("SAM<%d>:\n",i); // sam[i].Print(); } for (int i=0;i<n;i++) { int now=1; int clen=0; for (int j=0;j<len[i];j++) { int w=str[i][j]-‘0‘; if (sam[now].nxt[w]) { now=sam[now].nxt[w]; clen++; }else { while (now && !sam[now].nxt[w]) now=sam[now].pnt; if (!now) { now=1; clen=0; }else { clen=sam[now].len+1; now=sam[now].nxt[w]; } } g[j]=clen; // printf("%d\n",clen); } for (int j=len[i];j>=1;j--) g[j]=g[j-1]; g[0]=0; int l=0,r=len[i]+1; int mid; int head,tail=0; int ny; while (l+1<r) { mid=(l+r)>>1; int t; for (int j=1;j<=len[i];j++) dp[j]=-INF; dp[0]=0; head=0,tail=-1; ny=0; for (int j=1;j<=len[i];j++) { x=j-g[j]; y=j-mid; while (ny<=y) { while (head<=tail && dp[seq[tail]]<=dp[ny]) tail--; seq[++tail]=ny++; } while (head<=tail && seq[head]<x) head++; dp[j]=dp[j-1]+(j-1); if (head<=tail) dp[j]=max(dp[j],dp[seq[head]]+j); dp[j]-=j; } t=dp[len[i]]+len[i]; if (t*10>=len[i]*9) l=mid; else r=mid; } printf("%d\n",l); } }
bzoj 2806: [Ctsc2012]Cheat 后缀自动机DP
标签：
原文地址：http://www.cnblogs.com/mhy12345/p/4433683.html
