template

Author: ffdy

【模板】AC自动机.cpp

@@ -0,0 +1,77 @@
+#include<cstdio>
+#include<algorithm>
+#include<cstring>
+#include<queue>
+using namespace std;
+const int maxn=1e6+10;
+int rt,n;
+char str[maxn];
+struct fy
+{
+	int son[26],fail,sum;
+}ac[maxn];
+void add(char str[])
+{
+	int l=strlen(str),root=0;
+	for(int i=0;i<l;i++)
+	{
+		int a=str[i]-'a';
+		if(!ac[root].son[a]) ac[root].son[a]=++rt;
+		root=ac[root].son[a];
+	}
+	ac[root].sum++;
+}
+void make()
+{
+	queue<int> q;
+	for(int i=0;i<26;i++)
+	{
+		if(ac[0].son[i])
+		{
+			ac[ac[0].son[i]].fail=0;
+			q.push(ac[0].son[i]);
+		}
+	}
+	while(!q.empty())
+	{
+		int a=q.front();q.pop();
+		for(int i=0;i<26;i++)
+		{
+			if(ac[a].son[i])
+			{
+				ac[ac[a].son[i]].fail=ac[ac[a].fail].son[i];
+				q.push(ac[a].son[i]);
+			}
+			else ac[a].son[i]=ac[ac[a].fail].son[i];
+		}
+	}
+}
+int acm(char str[])
+{
+	int l=strlen(str);
+	int root=0,ans=0;
+	for(int i=0;i<l;i++)
+	{
+		root=ac[root].son[str[i]-'a'];
+		for(int j=root;j&&ac[j].sum!=-1;j=ac[j].fail)
+		{
+			ans+=ac[j].sum;
+			ac[j].sum=-1;
+		}
+	}
+	return ans;
+}
+int main()
+{
+	scanf("%d",&n);
+	while(n--)
+	{
+		scanf("%s",str);
+		add(str);
+	}
+	ac[0].fail=0;
+	make();
+	scanf("%s",str);
+	printf("%d\n",acm(str));
+	return 0;
+}

【模板】AC自动机（加强版）.cpp

@@ -0,0 +1,78 @@
+#include<cstdio>
+#include<algorithm>
+#include<iostream>
+#include<cstring>
+#include<queue>
+using namespace std;
+const int maxn=2e4+10,maxl=1e6+10;
+int n,rt,bi;
+char str[160][80],mu[maxl];
+int ans[maxl];
+struct fy
+{
+	int son[26],fail,w;
+}ac[maxl];
+void add(char str[],int v)
+{
+	int l=strlen(str),root=0;
+	for(int i=0;i<l;i++)
+	{
+		int a=str[i]-'a';
+		if(!ac[root].son[a]) ac[root].son[a]=++rt;
+		root=ac[root].son[a];
+	}
+	ac[root].w=v;
+}
+void make()
+{
+	queue<int>q;
+	for(int i=0;i<26;i++)
+	if(ac[0].son[i]) q.push(ac[0].son[i]);
+	while(!q.empty())
+	{
+		int a=q.front();q.pop();
+		for(int i=0;i<26;i++)
+		{
+			if(ac[a].son[i]) 
+			{
+				ac[ac[a].son[i]].fail=ac[ac[a].fail].son[i];
+				q.push(ac[a].son[i]);
+			}
+			else ac[a].son[i]=ac[ac[a].fail].son[i];
+		}
+	}
+}
+void acm(char str[])
+{
+	int l=strlen(str),root=0;
+	for(int i=0;i<l;i++)
+	{
+		root=ac[root].son[str[i]-'a'];
+		for(int j=root;j;j=ac[j].fail)
+		if(ac[j].w)
+		{
+			ans[ac[j].w]++;
+			bi=max(bi,ans[ac[j].w]);
+		}
+	}
+}
+int main()
+{
+	while(scanf("%d",&n)==1&&n!=0)
+	{
+		memset(&ac,0,sizeof ac);
+		bi=0;
+		memset(ans,0,sizeof ans);
+		for(int i=1;i<=n;i++)
+		{
+			scanf("%s",str[i]);
+			add(str[i],i);
+		}
+		make();scanf("%s",mu);
+		acm(mu);
+		printf("%d\n",bi);
+		for(int i=1;i<=n;i++)
+		if(ans[i]==bi) cout<<str[i]<<endl;
+	}
+	return 0;
+}

【模板】KM(二分图带权匹配).cpp

@@ -0,0 +1,43 @@
+#include<cstdio>
+#include<algorithm>
+#include<cstring>
+using namespace std;
+const int maxn=30,inf=0x7fffffff;
+int n,v[maxn][maxn],v1[maxn][maxn],lx[maxn],ly[maxn],p[maxn],ans,minx;
+bool usex[maxn],usey[maxn];
+bool dfs(int a)
+{
+	usex[a]=true;
+	for(int i=1;i<=n;i++)if(!usey[i])
+	{
+		int t=lx[a]+ly[i]-v[a][i];
+		if(t==0)
+		{
+			usey[i]=true;
+			if(!p[i]||dfs(p[i]))
+			{p[i]=a;return true;}
+		}
+		else if(t)minx=min(minx,t);
+	}
+	return false;
+}
+void km()
+{
+	for(int i=1;i<=n;i++)while(1)
+	{
+		minx=inf;memset(usex,false,sizeof usex);memset(usey,false,sizeof usey);
+		if(dfs(i)) break;
+		for(int j=1;j<=n;j++)if(usex[j])lx[j]-=minx;
+		for(int j=1;j<=n;j++)if(usey[j])ly[j]+=minx;
+	}
+}
+int main()
+{
+	scanf("%d",&n);
+	for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)scanf("%d",&v[i][j]);
+	// ×îСƥÅä³õÖµ·´×ªÎª¸º 
+	for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)lx[i]=max(lx[i],v[i][j]);
+	km();
+	for(int i=1;i<=n;i++)ans+=v[p[i]][i];printf("%d\n",ans);
+	return 0;
+}

【模板】KMP字符串匹配.cpp

@@ -0,0 +1,39 @@
+#include<cstdio>
+#include<cstring>
+#include<algorithm>
+using namespace std;
+const int maxn=1e6+10;
+int n,m,next[maxn];
+char str[maxn],stt[maxn];
+void makenext()
+{
+	int j=0;
+	for(int i=1;i<m;i++)
+	{
+		while(j&&stt[j]!=stt[i])
+		j=next[j-1];
+		if(stt[j]==stt[i]) j++;
+		next[i]=j;
+	}
+}
+void kmp()
+{
+	int j=0;
+	makenext();
+	for(int i=0;i<n;i++)
+	{
+		while(j&&stt[j]!=str[i])
+		j=next[j-1];
+		if(stt[j]==str[i]) j++;
+		if(j==m) printf("%d\n",i-j+2);
+	}
+}
+int main()
+{
+	scanf("%s%s",str,stt);
+	n=strlen(str);m=strlen(stt);
+	kmp();
+	for(int i=0;i<m;i++)
+	printf("%d ",next[i]);
+	return 0;
+}

【模板】ST表.cpp

@@ -0,0 +1,27 @@
+#include<cstdio>
+#include<algorithm>
+using namespace std;
+const int maxn=1e5+10;
+int n,m,ma[maxn][40],log[maxn];
+int main()
+{
+	scanf("%d%d",&n,&m);
+	log[0]=-1;
+	for(int i=1;i<=n;i++)
+	{
+		scanf("%d",&ma[i][0]);
+		if((i&(i-1))==0) log[i]=log[i-1]+1;
+		else log[i]=log[i-1];
+	}
+	for(int j=1;j<=log[n];j++)
+	for(int i=1;i+(1<<j)-1<=n;i++)
+	ma[i][j]=max(ma[i][j-1],ma[i+(1<<(j-1))][j-1]);
+	for(int i=1;i<=m;i++)
+	{
+		int x,y;
+		scanf("%d%d",&x,&y);
+		int k=log[y-x+1];
+		printf("%d\n",max(ma[x][k],ma[y-(1<<k)+1][k]));
+	}
+	return 0;
+}

【模板】floyd最小环.cpp

@@ -0,0 +1,58 @@
+#include<cstdio>
+#include<cstring>
+#include<algorithm>
+#define ll long long
+using namespace std;
+const int maxn=210,inf=0x3f3f3f3f;
+int n,m,xx[maxn],w,fa[maxn][maxn];
+ll map[maxn][maxn],dis[maxn][maxn],ans=inf;
+void floyd()
+{
+	ll he;
+	for(int k=1;k<=n;k++)
+	{
+		for(int i=1;i<k;i++)for(int j=i+1;j<k;j++)
+		{
+			he=dis[i][j]+map[j][k]+map[k][i];
+			if(he<ans)
+			{
+				ans=he;w=0;int a=j;
+				while(a!=i){xx[++w]=a;a=fa[i][a];}
+				xx[++w]=i;xx[++w]=k;
+			}
+		}
+		for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)
+		{
+			he=dis[i][k]+dis[k][j];
+			if(he<dis[i][j])
+			{
+				dis[i][j]=he;
+				fa[i][j]=fa[k][j];
+			}
+		}
+	}
+}
+int main()
+{
+//	freopen("trip7.in","r",stdin);
+	scanf("%d%d",&n,&m);int a,b;ll c;
+	for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)
+	{
+		map[i][j]=dis[i][j]=inf;
+		fa[i][j]=i;
+	}
+	for(int i=1;i<=m;i++)
+	{
+		scanf("%d%d%lld",&a,&b,&c);c=min(c,dis[a][b]);
+		dis[a][b]=dis[b][a]=c;
+		map[a][b]=map[b][a]=dis[a][b];
+	}
+	floyd();
+	if(ans==inf)
+	{
+		printf("No solution.\n");
+		return 0;
+	}
+	for(int i=1;i<=w;i++)printf("%d ",xx[i]);
+	return 0;
+}

【模板】manacher算法.cpp

@@ -0,0 +1,33 @@
+#include<cstdio>
+#include<algorithm>
+#include<cstring>
+using namespace std;
+const int maxn=11e6+10;
+char str[maxn],strr[maxn<<1];
+int n;
+int len[maxn<<1];
+int main()
+{
+	scanf("%s",str);
+	n=strlen(str);
+	strr[0]='@';
+	for(int i=0;i<=n;i++)
+	{
+		strr[(i<<1)+2]=str[i];
+		strr[(i<<1)+1]='#';
+	}
+	int mx=0,ans=0,mid=0;
+	for(int i=1;i<=(n<<1)+2;i++)
+	{
+		len[i]=(mx>i ? min(mx-i,len[mid*2-i]):1);
+		while(strr[i-len[i]]==strr[i+len[i]]) len[i]++;
+		if(i+len[i]>mx)
+		{
+			mx=len[i]+i;
+			mid=i;
+		} 
+		ans=max(ans,len[i]);
+	}
+	printf("%d",ans-1);
+	return 0;
+}