模板汇总
本博客所有模板都经过测试，保证正确。
归并排序
#includeLCS
#include#include
using namespace std;
const int N=1e5+77;
int n,f[N],a[N],b[N];
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++) scanf("%d",&a[i]);
for(int i=1;i<=n;i++) scanf("%d",&b[i]);
for(int i=1;i<=n;i++) f[b[i]]=i,b[i]=0;
for(int i=1;i<=n;i++) a[i]=f[a[i]];
int len=0;
for(int i=1;i<=n;i++)
{
if(a[i]>b[len]) b[++len]=a[i];
else b[upper_bound(b+1,b+len+1,a[i])-b]=a[i];
}
printf("%d\n",len);
return 0;
} 小B的询问（莫队）
#include
using namespace std;
const int N=5e4+77;
struct Q{
int qv,id,l,r;
}q[N];
int n,m,k,a[N],ans[N],cnt[N],sum;
bool cmp(Q x,Q y)
{
if(x.qv==y.qv) return x.r<y.r;
return x.qv<y.qv;
}
void add(int x)
{
if(!x) return;
sum-=cnt[x]*cnt[x];
cnt[x]++;
sum+=cnt[x]*cnt[x];
}
void del(int x)
{
if(!x) return;
sum-=cnt[x]*cnt[x];
cnt[x]--;
sum+=cnt[x]*cnt[x];
}
int main()
{
scanf("%d%d%d",&n,&m,&k);
for(int i=1;i<=n;i++) scanf("%d",&a[i]);
int limit=sqrt(n+10);
for(int i=1;i<=m;i++)
{
int l,r;
scanf("%d%d",&l,&r);
q[i].l=l,q[i].r=r;
q[i].qv=q[i].l/limit;
q[i].id=i;
}
sort(q+1,q+m+1,cmp);
int l=0,r=0;
for(int i=1;i<=m;i++)
{
while(r<q[i].r)
{
add(a[++r]);
}
while(r>q[i].r)
{
del(a[r--]);
}
while(l<q[i].l)
{
del(a[l++]);
}
while(l>q[i].l)
{
add(a[--l]);
}
ans[q[i].id]=sum;
}
for(int i=1;i<=m;i++)
printf("%d\n",ans[i]);
return 0;
}并查集
#include
using namespace std;
const int N=1e4+77;
int n,m,f[N];
int find(int x)
{
return f[x]==x?x:f[x]=find(f[x]);
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++) f[i]=i;
for(int i=1;i<=m;i++)
{
int opt,x,y;
scanf("%d%d%d",&opt,&x,&y);
int fx=find(x),fy=find(y);
if(opt==1)
{
if(fx==fy) continue;
f[fx]=fy;
}
else
{
if(fx==fy) printf("Y\n");
else printf("N\n");
}
}
return 0;
}最小生成树
#includeint find(int x)
{
return f[x]==x?x:f[x]=find(f[x]);
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++) f[i]=i;
for(int i=1;i<=m;i++)
{
scanf("%d%d%d",&e[i].x,&e[i].y,&e[i].z);
}
sort(e+1,e+m+1,cmp);
int cnt=n;
for(int i=1;i<=m;i++)
{
int fx=find(e[i].x),fy=find(e[i].y);
if(fx==fy) continue;
f[fx]=fy;
ans+=e[i].z;
cnt--;
if(cnt==1) break;
}
if(cnt==1) printf("%d\n",ans);
else printf("orz\n");
return 0;
} 字符串hash
#include*b1%mo+s[i]-'0')%mo;
return h;
}
LL Hash()
{
LL h=0;
int len=strlen(s);
for(int i=0;i*b2%moo+s[i]-'0')%moo;
return h;
}
int main()
{
scanf("%d\n",&n);
for(int i=1;i<=n;i++)
{
gets(s);
q[i].a=hash();
q[i].b=Hash();
}
q[0].a=q[0].b=-1;
sort(q+1,q+n+1,cmp);
for(int i=1;i<=n;i++)
if(q[i].a!=q[i-1].a||q[i].b!=q[i-1].b)
ans++;
printf("%d\n",ans);
return 0;
} 最短路spfa O(K*m/RP)
#include最短路Dijsktra O(n^n)
#includeint minn=inf,x=0;
for(int j=1;j<=n;j++)
{
if(!vis[j]&&minn>dis[j])
minn=dis[j],x=j;
}
if(!x) break;
vis[x]=1;
for(int j=head[x];j;j=nxt[j])
if(dis[to[j]]>dis[x]+cost[j])
dis[to[j]]=dis[x]+cost[j];
}
for(int i=1;i<=n;i++)
printf("%d ",dis[i]);
return 0;
} 树状数组（单点修改+区间求和）
#include线段树 1 区间修改+区间求和
#includeST表
#includefor(int i=1;i<=n;i++)
{
if(i+(1<1<=n)
MAX[i][j]=max(MAX[i][j-1],MAX[i+(1<1)][j-1]);
}
}
for(int i=1;i<=m;i++)
{
int l,r;
scanf("%d%d",&l,&r);
int t=log2(r-l+1);
ans=max(MAX[l][t],MAX[r-(1<1][t]);
printf("%d\n",ans);
}
return 0;
} KMP
#includeif(j==-1||b[i]==b[j])
nxt[++i]=++j;
else
j=nxt[j];
}
}
void match(int i,int j)
{
while(iwhile(jif(j==-1||a[i]==b[j])
i++,j++;
else j=nxt[j];
}
if(j==len2) ans[++cnt]=i-len2+1,j=nxt[j];
}
}
int main()
{
cin>>a;len1=strlen(a);
cin>>b;len2=strlen(b);
next_make();
match(0,0);
for(int i=1;i<=cnt;i++)
printf("%d\n",ans[i]);
for(int i=1;i<=len2;i++)
printf("%d ",nxt[i]);
return 0;
} 线性乘法逆元
证明
#includeTarjan求LCA
#includeint x,y;
scanf("%d%d",&x,&y);
add(x,y);add(y,x);
}
for(int i=1;i<=m;i++)
{
int x,y;
scanf("%d%d",&x,&y);
qadd(x,y);qadd(y,x);
}
dfs(s);
for(int i=1;i<=m;i++)
{
printf("%d\n",lca[2*i]);
}
return 0;
} 倍增求LCA
#includeint i,j;
for(i=20;i>=0;i--)
{
if(dep[x]==dep[y]) break;
if(f[x][i])
if(dep[f[x][i]]>=dep[y])
x=f[x][i];
}
if(x==y) return x;
for(i=20;i>=0;i--)
{
if(f[x][i]&&f[y][i]&&f[x][i]!=f[y][i])//一直跳到下一步就能汇合
x=f[x][i],y=f[y][i];
}
return f[x][0];
}
int main()
{
scanf("%d%d%d",&n,&m,&s);
for(int i=1;iint x,y;
scanf("%d%d",&x,&y);
add(x,y);add(y,x);
}
dfs(s,0);
for(int j=1;(1<for(int i=1;i<=n;i++)
{
if(f[f[i][j-1]][j-1])
f[i][j]=f[f[i][j-1]][j-1];
}
}
for(int i=1;i<=m;i++)
{
int x,y,ans;
scanf("%d%d",&x,&y);
ans=find(x,y);
printf("%d\n",ans);
}
return 0;
} Tarjan缩点+树形dp（也可以用spfa）
#include二分图最大匹配（匈牙利算法）
#includedfs_spfa判断负环
其思想是：采用广度优先的思想，每当我们扩展出一个新的节点，总是把它放到队列的末尾，其缺点是中断了迭代的连续性。而实际上如果采用深度优先的思想，我们可以直接从这个新节点继续往下扩展。于是算法的实现方式可以改成不断从新节点往下递归进行求解。 而对于负环的判断则显得更为简单，因为假如存在负环a1->a2->….ak->a1，那么算法运行时，会从某个点a1开始Dfs，最后又回到了这个点。所以只需用一个辅助数组记录当前节点是否在递归栈中便可及时检测出负环。
#include欧拉函数的线性筛（附带素数线性筛）
#includeif(!np[i])
{
prime[++cnt]=i;
phi[i]=i-1;
}
for(int j=1;j<=cnt&&1ll*i*prime[j]1;
if(i%prime[j]==0)
{
phi[i*prime[j]]=phi[i]*prime[j];
break;
}
else phi[i*prime[j]]=phi[i]*(prime[j]-1);
}
}
}
int main()
{
PHI();
int n;
scanf("%d",&n);
for(int i=1;i<=n;i++)
printf("%d\n",phi[i]);
return 0;
} 同余方程
#include高精乘高精
#include低精除高精
#include高精乘高精
#include矩阵快速幂（多次询问离线）
蒟蒻不会重载运算符qwq
#includex;
}
LL q_pow(LL p)
{
p--;
for(int i=1;i<=3;i++) for(int j=1;j<=3;j++) B[i][j]=C[i][j]=A[i][j];
while(p)
{
if(p&1)
{
for(int i=1;i<=3;i++)
for(int j=1;j<=3;j++)
tmp[i][j]=B[i][j],B[i][j]=0;
for(int i=1;i<=3;i++)
for(int j=1;j<=3;j++)
for(int k=1;k<=3;k++)
(B[i][j]+=(tmp[i][k]*C[k][j])%MOD)%=MOD;
}
for(int i=1;i<=3;i++)
for(int j=1;j<=3;j++)
tmp[i][j]=C[i][j],C[i][j]=0;
for(int i=1;i<=3;i++)
for(int j=1;j<=3;j++)
for(int k=1;k<=3;k++)
(C[i][j]+=(tmp[i][k]*tmp[k][j])%MOD)%=MOD;
p>>=1;
}
for(int i=1;i<=3;i++)
for(int j=1;j<=3;j++)
tmp[i][j]=ans[i][j],ans[i][j]=0;
for(int i=1;i<=3;i++)
for(int j=1;j<=3;j++)
for(int k=1;k<=3;k++)
(ans[i][j]+=(tmp[i][k]*B[k][j])%MOD)%=MOD;
return ans[1][1];
}
int main()
{
scanf("%d",&n);
ans[1][1]=1,ans[1][2]=1,ans[1][3]=1;
for(int i=1;i<=n;i++)
{
scanf("%d",&q[i].x);
q[i].id=i;
}
sort(q+1,q+n+1,cmp);
bool first=1;
for(int i=1;i<=n;i++)
{
if(q[i].x<=3)
{
b[q[i].id]=1;
continue;
}
if(first) b[q[i].id]=q_pow(q[i].x-3),first=0;
else
{
if(q[i].x==q[i-1].x) b[q[i].id]=b[q[i-1].id];
else
b[q[i].id]=q_pow(q[i].x-q[i-1].x);
}
}
for(int i=1;i<=n;i++)
printf("%lld\n",b[i]);
return 0;
}