秒

宏定义（下面部分用旧版本）：

#define rep(i,l,r) for(int i=(l),_=(r);i<=_;i++)
#define per(i,r,l) for(int i=(r),_=(l);i>=_;i--)
#define MS(arr,x) memset(arr,x,sizeof(arr))
#define INE(i,u) for(int i=head[u];~i;i=e[i].next)
#define LL long long

 

01背包：25s

for(int i=1;i<=n;i++)
for(int j=0;j<=v;j++)
{
   opt[i][j]=opt[i-1][j];
   if(j-v[i]>0) opt[i][j]=max(opt[i-1][j],opt[i-1][j-c[i]]+w[i]);
}

for(int i=1;i<=n;i++)
for(int j=v;j>=0;j--)
    if(j-c[i]>=0) opt[j]=max(opt[j],opt[j-c[i]]+w[i]);

 

完全背包：15s

rep(i,1,n) rep(j,w[i],m)
    opt[j]=max(opt[j],opt[j-w[i]]+q[i]);

 

 

LIS：40s

下降：两个>改成<    不降：两个>改成>=

for(int i=1;i<=n;i++)
{
   if(a[i]>b[cnt]) b[++cnt]=a[i];
   else
   {
      l=0,r=cnt;
      while(l<r)
      {
        mid=(l+r)>>1;
        if(a[i]>b[mid]) l=mid+1;
        else            r=mid;
      }
      b[l]=a[i];
   }
}

归并&逆序对：60s

void merge_sort(int l,int r)
{
     if(l>=r)return;
     int mid=(l+r)>>1;
     int p=l,q=mid+1,i=l;
     merge_sort(l,mid);
     merge_sort(mid+1,r);
     while(p<=mid||q<=r)
     {
         if(q>r||(p<=mid&&a[p]<a[q]))
             t[i++]=a[p++];
         else
         {
             t[i++]=a[q++];
             cnt+=mid-p+1;
         }
     }
     rep(i,l,r) a[i]=t[i];
}

 

快速幂：25s

现在熟练了，以后就一行了

 怕a*a或者t*a爆了的话见下面的

LL pow_mod(LL a,LL b,LL p){LL t=1;for(;b;b>>=1,a=a*a%p)if(b&1)t=t*a%p;return t;}

（快速乘）：

这个也升级啦，尽量把小的放前面

ll mul(ll a,ll b,ll p){ll t=0;for(;b;b>>=1,a=(a<<1)%p)if(b&1)t=(t+a)%p;return t;}

扩展欧几里得：30s

神奇的代码缩减

LL exgcd(LL a,LL b,LL &x,LL &y)
{
    if(!b) return x=1,y=0,a;
    LL d=exgcd(b,a%b,x,y);
    return swap(x,y),y-=a/b*x,d;
}

floyd+最小环：40s

rep(k,1,n)
{
    rep(i,1,k-1) rep(j,i+1,k-1)
        ans=min(ans,dis[i][j]+map[i][k]+map[k][j]);
    rep(i,1,n) rep(j,1,n)
        dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
}

 

 Phi：

LL phi(LL n)
{
    LL res=n;
    LL sqn=(LL)sqrt(n+0.5);
    for(LL i=2;i<=sqn;i==2?i++:i+=2) if(n%i==0)
    {
        res=res/i*(i-1);
	while(n%i==0)n/=i;
    }
    if(n>1)res=res/n*(n-1);
    return res;
}

线性筛素数&欧拉&莫比乌斯：100s

void get_prim_phi_mu()
{
    cnt=0; MS(flag,0); phi[1]=1; mu[1]=1;
    rep(i,2,n)
    {
        if(!flag[i])
            prim[++cnt]=i,phi[i]=i-1,mu[i]=-1;
        for(int j=1;j<=cnt && i*prim[j]<=n;j++)
        {
            flag[i*prim[j]]=1;
            if(i%prim[j]==0)
            {
                phi[i*prim[j]]=phi[i]*prim[j];
                mu[i*prim[j]]=0;
                break;
            }
            else
            {
            	phi[i*prim[j]]=phi[i]*(prim[j]-1);
            	mu[i*prim[j]]=-mu[i];
            }
        }
    }
}

 bsgs：学了再放

输入数据

inline const int read()
{
    int r=0,k=1;char c=getchar();
    for(;c<'0'||c>'9';c=getchar())if(c=='-')k=-1;
    for(;c>='0'&&c<='9';c=getchar())r=r*10+c-'0';
    return k*r;
}

组合数取模：30s

MS(C,0);
rep(i,0,n) C[i][0]=1;
rep(i,1,n) rep(j,1,i)
   C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;

求一个的话需要逆元：40s

inline LL C(int n,int m)
{
    if(m+m>n) m=n-m;
    LL s1=1,s2=1;
    rep(i,1,m) s1=(s1*(n-i+1))%mod,
               s2=(s2*i)*mod;
    return (s1*pow_mod(s2,mod-2,mod))%mod;
}

如果范围较小而p较大，可以预处理阶乘、阶乘逆元

据说逆元用pow_mod会被卡，特意收集了VFK的模板

inline int modinv(const int &a)
{
    int x1 = 1, x2 = 0, x3 = Mod;
    int y1 = 0, y2 = 1, y3 = a;
    while (y3 != 1)
    {
        int k = x3 / y3;
        x1 -= y1 * k, x2 -= y2 * k, x3 -= y3 * k;
        swap(x1, y1), swap(x2, y2), swap(x3, y3);
    }
    return y2 >= 0 ? y2 : y2 + Mod;
}

KMP:

rep(i,2,n)
{
    while(j&&s[i]!=s[j+1]) j=next[j];
    next[i]=j+=s[i]==s[j+1];
}

EXKMP:

int t=0; ext[1]=n;
rep(i,2,n)
{
	ext[i]=max(0,min(ext[i-t+1],t+ext[t]-i));
	while(s[ext[i]+1] == s[i+ext[i]]) ext[i]++;
	if(i+ext[i] > t+ext[t]) t=i;
}