各种简单的困难的模板，持续更新

大部分模板为代码片段，可能有少部分的完整程序

---

邻接表

struct os{    int fa,son,next,v;}a[100000];//边数int first[10000]//点数//定义void add(int x,int y,int z){    a[++tot].fa=x;    a[tot].son=y;    a[tot].v=z;    a[tot].next=first[x];    first[x]=tot;}……    for (int i=1;i<=m;i++)    scanf("%d%d%d",&x,&y,&z),    add(x,y,z),    //add(y,x,z);有向图时不用这条//建边for (int i=first[x];i;i=a[i].next)……//调用所有以x为出点的边

矩阵乘法（mod为取模数）

struct matrix{    int map[130][130];}；matrix mul(matrix a,matrix b){    matrix c;    memset(c.map,0,sizeof(c.map));    for (int i=1;i<=n;i++)    for (int j=1;j<=n;j++)    for (int k=1;k<=n;k++)    c.map[i][j]=(c.map[i][j]+(a.map[i][k]*b.map[k][j]%mod))%mod;    return c;}

矩阵快速幂（配合上面使用）

matrix qr(matrix a,int b){    matrix c;    c=a;    b--;    while (b)    {        if (b&1) c=mul(c,a);        a=mul(a,a);        b>>=1;    }    return c;}

普通快速幂

int qr(int x,int y,int z){    int ans=1;    x%=z;    while (y)    {        if (y&1) ans=ans*x%z;        x=x*x%z;        y>>=1;    }    return ans;}

欧拉筛（求1-n的phi和质数）

    int prime[1000000],phi[1000000]    ……    phi[1]=1;    for (int i=2;i<=n;i++)    {        if (!pd[i])        prime[++prime[0]]=i,phi[i]=i-1;        for (int j=1;j<=prime[0];j++)        {            if (prime[j]*i>n) break;            pd[prime[j]*i]=1;            if (i%prime[j]==0) {phi[i*prime[j]]=phi[i]*prime[j];break;}            else phi[i*prime[j]]=phi[i]*phi[prime[j]];        }    }

求逆元

void exgcd(int a,int b,int &x,int &y){    if (b==0) {x=1;y=0;return;}    exgcd(b,a%b,x,y);    int x1=y,y1=x-a/b*y;    x=x1;y=y1;}int getinv(int p){    int x,y;    exgcd(p,r,x,y);    return (x%r+r)%r;}

位运算小技巧：http://blog.csdn.net/xym_csdn/article/details/50725540
Lucas定理

……#define LL long long……LL lucas(LL p,LL q,LL num){    LL sum=1;    while (p&&q)    {        LL x=p%prime[num],y=q%prime[num];        if (y>x) return 0;        sum=(sum*(fac[x]*qr(fac[y]*fac[x-y],prime[num]-2,prime[num]))%prime[num])%prime[num];        p/=prime[num];        q/=prime[num];    }    return sum;}……for (LL j=1;j<=prime[i];j++) fac[j]=fac[j-1]*j%mod;

SPFA（邻接表）

#include<bits/stdc++.h>using namespace std; struct os{    int fa,son,next,w;}a[100010];bool flag[10010];int first[10010],tot,dis[10010],n,m,x,y,z;queue <int>q;void add(int x,int y,int z){    a[++tot].fa=x;    a[tot].son=y;    a[tot].w=z;    a[tot].next=first[x];    first[x]=tot;}main(){    scanf("%d%d",&n,&m);    for (int i=1;i<=m;i++)    scanf("%d%d%d",&x,&y,&z),    add(x,y,z);//  add(y,x,z);//是无向图就去掉注释    memset(dis,63,sizeof(dis));    dis[1]=0;    q.push(1);    flag[1]=1;    while (!q.empty())    {        int t=q.front();        flag[t]=0;        q.pop();        for (int i=first[t];i;i=a[i].next)        {            int p=a[i].son;            if (dis[t]+a[i].w<dis[p])            {                dis[p]=dis[t]+a[i].w;                if (!flag[p]) q.push(p),flag[p]=1;            }        }    }    printf("%d",dis[n]);}

快速读入

int in(){    int f=1,t=0;    char ch=getchar();    while (ch>'9'||ch<'0')    {        if (ch=='-') f=-1;        ch=getchar();    }    while (ch>='0'&&ch<='9') t=t*10+ch-'0',ch=getchar();    return f*t;}

秦九韶算法

#include<bits/stdc++.h>#define LL long longusing namespace std;int n,x,a[100000];LL ans;main(){    scanf("%d%d",&n,&x);    for(int i=0;i<=n;i++) scanf("%d",&a[i]);    for (int i=n;i>=0;i--)    ans=ans*x+a[i];    printf("%lld",ans);}

单调队列
给定一个长度为N的整数数列a(i),i=0,1,…,N-1和窗长度k.

要求：

f(i)=max(a(i−k+1),a(i−k+2),...,a(i))i=1,2,3,..n

#include<bits/stdc++.h>using namespace std;int x,n,k,head=1,tail;struct os{    int data,num;}q[500000];main(){    scanf("%d%d",&n,&k);    for (int i=1;i<=n;i++)    {        scanf("%d",&x);        if (head<=tail&&i-q[head].num>=k) head++;        while (head<=tail)        {            if (q[tail].data>x) break;            tail--;        }        q[++tail].data=x;        q[tail].num=i;        if (i>=k) printf("%d ",q[head].data);    }}

邻接表+堆优化dijskra求最短路
点对中的优先级比较first>second，所以把dis当成first

#include<iostream>#include<cstdio>#include<cstring>#include<queue>using namespace std;int tot,n,m;int first[10010],dis[10010];bool flag[10010];struct os{    int w,u,v,next;}e[200010];void add(int x,int y,int z){    e[++tot].u=x;    e[tot].v=y;    e[tot].w=z;    e[tot].next=first[x];    first[x]=tot;}typedef pair<int,int> xy;priority_queue<xy,vector<xy>,greater<xy> >q; main(){    scanf("%d%d",&n,&m);    int x,y,z;    for (int i=1;i<=m;i++)    scanf("%d%d%d",&x,&y,&z),    add(x,y,z),add(y,x,z);//无向图    memset(dis,63,sizeof(dis));    dis[1]=0;    q.push(make_pair(dis[1],1));    while (!q.empty())    {        xy now=q.top();        q.pop();        if (flag[now.second]) continue;        flag[now.second]=1;        for (int i=first[now.second];i;i=e[i].next)        if (dis[now.second]+e[i].w<dis[e[i].v])        {            dis[e[i].v]=dis[now.second]+e[i].w;            if (!flag[e[i].v]) q.push(make_pair(dis[e[i].v],e[i].v));        }    }    printf("%d",dis[n]);}

kmp算法（求ch在s中第一次出现的位置）

#include<cstdio>#include<cstdlib>#include<cstring>#include<iostream>char s[100002],ch[102];int l,next[102];main(){    scanf("%s%s",s,ch);    for (int i=1;i<strlen(ch);i++)    {        int k=i;        while (ch[i]!=ch[next[k]]&&next[k]) k=next[k];        next[i+1]=next[k]+(ch[i]==ch[next[k]]);    }    int now=0;    for (int i=0;i<strlen(s);i++)    {        while (ch[now]!=s[i]&&now) now=next[now];        now+=(ch[now]==s[i]);        if (now==strlen(ch)) {printf("%d\n",i-now+1);break;}    }}

AC自动机：AC自动机练习
LCT：弹飞绵羊
莫队：小Z的袜子
Splay（这个特别丑，可以看别的题解）：普通平衡树
背包：背包小例题
欧拉筛求φ，μ，d

#include<cstdio>using namespace std;int n,prime[1000010];bool flag[1000010];int t[1000010],d[1000010],mu[1000010],phi[1000010];main(){    scanf("%d",&n);    mu[1]=1;d[1]=1;    for (int i=2;i<=n;i++)    {        if (!flag[i])            prime[++prime[0]]=i,            mu[i]=-1,            phi[i]=i-1,            t[i]=i,            d[i]=2;        for (int j=1;j<=prime[0];j++)        {            if (i*prime[j]>n) break;            flag[prime[j]*i]=1;            if (i%prime[j])            {                phi[i*prime[j]]=phi[i]*(prime[j]-1);                mu[i*prime[j]]=-mu[i];                t[i*prime[j]]=1;                d[i*prime[j]]=d[i]<<1;            }            else            {                phi[i*prime[j]]=phi[i]*prime[j];                mu[i*prime[j]]=0;                t[i*prime[j]]=t[i]+1;                d[i*prime[j]]=d[i]/(t[i]+1)*(t[i]+2);                break;            }        }    }}

O(n)求1-n的逆元（mod为质数）
证明

#include<cstdio>using namespace std;int n,mod,inv[100000];main(){    scanf("%d%d",&n,&mod);    inv[1]=1;    printf("1 ");    for (int i=2;i<=n;i++)    inv[i]=(mod-mod/i)*inv[mod%i]%mod,    printf("%d ",inv[i]);}

匈牙利

bool find(int x){    for (int i=first[x];i;i=e[i].next        if (!vis[e[i].v])        {            vis[e[i].v]=1;            if (!belong[e[i].v]||find(belong[e[i].v]))            {                belong[e[i].v]=x;                return 1;            }        }    return 0;}

tarjan
定义：

int low[],dfn[],cnt,first[];bool vis[];stack<int>sstruct edge{int v,next;}

操作：

void dfs(int x){    low[x]=dfn[x]=++cnt;    vis[x]=1;s.push(x);    for (int i=first[x];i;i=e[i].next)        if (!dfn[e[i].v])            dfs(e[i].v),            low[x]=min(low[x],low[e[i].v]);        else if (vis[e[i].v])            low[x]=min(low[x],dfn[e[i].v];    if (dfn[x]==low[x])        for (int y=-1;y!=x;y=s.top(),s.pop())            vis[s.top()]=1,            //...}