树论
LCA
#include <bits/stdc++.h>
using namespace std;
const int maxn=5*1e6+5;
int f[maxn][25],dep[maxn],head[maxn],nxt[maxn],to[maxn],cnt,N;
void add(int x,int y)
{
to[++cnt]=y;
nxt[cnt]=head[x];
head[x]=cnt;
}
int lca(int x,int y)
{
if(dep[x]<dep[y]) swap(x,y);
for(int i=20;i>=0;i--) if(dep[f[x][i]]>=dep[y]) x=f[x][i];
if(x==y) return x;
for(int i=20;i>=0;i--)
if(f[x][i]!=f[y][i])
x=f[x][i],y=f[y][i];
return f[x][0];
}
void dfs(int x,int fa)
{
dep[x]=dep[fa]+1;
f[x][0]=fa;
for(int i=head[x];i;i=nxt[i])
{
if(to[i]==fa) continue;
dfs(to[i],x);
}
}
int main()
{
ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
int M,S,x,y;cin>>N>>M>>S;
for(int i=1;i<N;i++) cin>>x>>y,add(x,y),add(y,x);
dfs(S,0);int a,b;
for(int j=1;j<=20;j++) for(int i=1;i<=N;i++) f[i][j]=f[f[i][j-1]][j-1];
for(int i=1;i<=M;i++) cin>>a>>b,cout<<lca(a,b)<<endl;
return 0;
}
笛卡尔树
#include <bits/stdc++.h>
using namespace std;
#define int long long
const int maxn=1e7+5;
int l[maxn],r[maxn],p[maxn],stk[maxn];
inline int read()
{
register int x(0),f(1);register char c(getchar());
while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
while(c>='0'&&c<='9'){x=(x<<1)+(x<<3)+(c^'0');c=getchar();}
return x*f;
}
signed main()
{
int n;cin>>n;
for(int i=1;i<=n;i++) p[i]=read();
stk[1]=0;
int top=1;
for(int i=1;i<=n;i++)
{
int k=top;
while(k&&p[stk[k]]>p[i]) k--;
if(k) r[stk[k]]=i;
if(k<top) l[i]=stk[k+1];
stk[++k]=i,top=k;
}
int ans1=0,ans2=0;
for(int i=1;i<=n;i++) ans1^=(i*(l[i]+1)),ans2^=(i*(r[i]+1));
cout<<ans1<<' '<<ans2;
return 0;
}
图论
传递闭包
#include <bits/stdc++.h>
using namespace std;
const int maxn=105;
bool a[maxn][maxn],b[maxn][maxn];
int main()
{
int n;cin>>n;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
cin>>a[i][j];
for(int k=1;k<=n;k++)
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++) a[i][j]|=(a[i][k]&a[k][j]);
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++) cout<<a[i][j]<<' ';
cout<<endl;
}
return 0;
}
拓扑排序
#include <bits/stdc++.h>
using namespace std;
#define int long long
const int maxn=105;
struct edge{int to,nxt;}e[maxn*maxn];
int head[maxn],cnt,deg[maxn],N;
queue<int> q;
void add(int x,int y){e[++cnt]=(edge){y,head[x]},head[x]=cnt;}
void topo()
{
for(int i=1;i<=N;i++)
if(!deg[i]) q.push(i),cout<<i<<' ';
while(!q.empty())
{
int x=q.front();q.pop();
for(int i=head[x];i;i=e[i].nxt)
{
deg[e[i].to]--;
if(!deg[e[i].to]) cout<<e[i].to<<' ',q.push(e[i].to);
}
}
}
signed main()
{
cin>>N;
for(int i=1,a;i<=N;i++)
while(cin>>a&&a) add(i,a),deg[a]++;
topo();
return 0;
}
差分约束
#include <bits/stdc++.h>
using namespace std;
const int maxn=1e4+5;
int head[maxn],to[maxn],nxt[maxn],w[maxn],dis[maxn],cnt,n,m,tot[maxn];
bool vis[maxn];
void add(int x,int y,int z)
{to[++cnt]=y,w[cnt]=z,nxt[cnt]=head[x],head[x]=cnt;}
queue<int> q;
bool spfa(int s)
{
dis[s]=0,vis[s]=1,q.push(s);tot[s]++;
while(!q.empty())
{
int x=q.front();q.pop();vis[x]=0;
for(int i=head[x];i;i=nxt[i])
if(dis[to[i]]>dis[x]+w[i])
{
dis[to[i]]=dis[x]+w[i];
if(!vis[to[i]])
{
vis[to[i]]=1,tot[to[i]]++;q.push(to[i]);
if(tot[to[i]]==n+1) return 1;
}
}
}
return 0;
}
int main()
{
memset(dis,0x3f3f3f,sizeof dis);
cin>>n>>m;
for(int i=1;i<=n;i++) add(0,i,0);
for(int i=1;i<=m;i++)
{
int c,cc,y;cin>>c>>cc>>y;
add(cc,c,y);
}
if(spfa(0)) cout<<"NO";
else for(int i=1;i<=n;i++) cout<<dis[i]<<' ';
return 0;
}
Tarjan
缩点
#include <bits/stdc++.h>
using namespace std;
const int maxn=1e4+5,maxm=1e5+5;
int head1[maxn],head2[maxn],a[maxn],low[maxn],scc[maxn]/*sizeof scc*/,dfscnt,dfn[maxn],f[maxn],s[maxn],top,scccnt,cnt1,cnt2,sccnum[maxn];
struct edge{int to,nxt;}e1[maxm],e2[maxm];
bool vis[maxn];
void add1(int x,int y){e1[++cnt1]=(edge){y,head1[x]},head1[x]=cnt1;}
void add2(int x,int y){e2[++cnt2]=(edge){y,head2[x]},head2[x]=cnt2;}
void tarjan(int x)
{
low[x]=dfn[x]=++dfscnt,s[++top]=x,vis[x]=1;
for(int i=head1[x];i;i=e1[i].nxt)
{
if(!dfn[e1[i].to]) tarjan(e1[i].to),low[x]=min(low[x],low[e1[i].to]);
else if(vis[e1[i].to]) low[x]=min(low[x],dfn[e1[i].to]);
}
if(dfn[x]==low[x])
{
sccnum[x]=++scccnt;
while(s[top]!=x) sccnum[s[top]]=scccnt,vis[s[top--]]=0;
vis[x]=0,top--;
}
}
int main()
{
int n,m;cin>>n>>m;
for(int i=1;i<=n;i++) cin>>a[i];
for(int i=1,u,v;i<=m;i++) cin>>u>>v,add1(u,v);
for(int i=1;i<=n;i++) if(!dfn[i]) tarjan(i);
for(int i=1;i<=n;i++)
{
scc[sccnum[i]]+=a[i];
for(int j=head1[i];j;j=e1[j].nxt) if(sccnum[i]!=sccnum[e1[j].to]) add2(sccnum[i],sccnum[e1[j].to]);
}
int ans=-1;
for(int ii=scccnt;ii;ii--)
{
f[ii]+=scc[ii],ans=max(ans,f[ii]);
for(int i=head2[ii];i;i=e2[i].nxt) f[e2[i].to]=max(f[e2[i].to],f[ii]);
}
cout<<ans;
return 0;
}
割点
#include <bits/stdc++.h>
using namespace std;
const int maxn=2e5+5;
int head[maxn],to[maxn],nxt[maxn],cnt,dfscnt,low[maxn],dfn[maxn];
bool cut[20005];
void add(int x,int y)
{
to[++cnt]=y;
nxt[cnt]=head[x];
head[x]=cnt;
}
void dfs(int u,int fa)
{
low[u]=dfn[u]=++dfscnt;
int c=0;
for(int i=head[u];i;i=nxt[i])
{
int v=to[i];
if(!dfn[v])
{
dfs(v,u);
low[u]=min(low[u],low[v]);//树边
if(low[v]>=dfn[u]&&u!=fa) cut[u]=1;
if(u==fa) c++;
}
else low[u]=min(low[u],dfn[v]);//其他
}
if (c>=2&&u==fa) cut[u]=1;
}
int main()
{
int n,m,ans=0;cin>>n>>m;
for(int i=1;i<=m;i++) {int y,x;cin>>x>>y,add(x,y),add(y,x);}
for(int i=1;i<=n;i++) if(!dfn[i]) dfs(i,i);
for(int i=1;i<=n;i++) if(cut[i]) ans++;
cout<<ans<<endl;
for(int i=1;i<=n;i++) if(cut[i]) cout<<i<<' ';
return 0;
}
2-SAT
#include <bits/stdc++.h>
using namespace std;
const int maxn=2*1e6+5;
int head[maxn],to[maxn],nxt[maxn],cnt;
int dfn[maxn],sta[maxn]/*栈*/,low[maxn],dfncnt,vis[maxn],top/*记录栈中元素编号*/,sccnum/*记录强连通分量编号*/,scc[maxn]/*编号为i的点所在的强连通分量编号*/;
void add(int x,int y)
{
to[++cnt]=y;
nxt[cnt]=head[x];
head[x]=cnt;
}
void tarjan(int u)
{
low[u]=dfn[u]=++dfncnt;
sta[++top]=u;vis[u]=1;
for(int i=head[u];i;i=nxt[i])
{
int v=to[i];
if(!dfn[v]) tarjan(v)/*继续往下搜*/,low[u]=min(low[u],low[v]);
else if(vis[v]) low[u]=min(low[u],dfn[v]);
}
if(dfn[u]==low[u])
{
++sccnum;
while(sta[top]!=u)//把它搜索树里的都弹完
scc[sta[top]]=sccnum,vis[sta[top]]=0,top--;//数组模拟弹栈
//处理它自己
scc[sta[top]]=sccnum;
vis[sta[top]]=0;
top--;
}
}
int main()
{
int n,m;cin>>n>>m;
for(int i=1;i<=m;i++)
{
int ii,a,jj,b;cin>>ii>>a>>jj>>b;
//实现建两个点，对于一个点k，k+n是true点，k是false点
add(ii+!a*n,jj+b*n);
add(jj+!b*n,ii+a*n);
}
for(int i=1;i<=2*n;i++) if(!dfn[i]) tarjan(i);
for(int i=1;i<=n;i++)
if(scc[i]==scc[i+n]) puts("IMPOSSIBLE"),exit(0);
puts("POSSIBLE");
for(int i=1;i<=n;i++)
{
if(scc[i]>scc[i+n]) cout<<1<<' ';
else cout<<0<<' ';
}
return 0;
}
边双
#include <bits/stdc++.h>
using namespace std;
#define int long long
const int maxn=4e6+5;
int head[maxn],to[maxn],nxt[maxn],cnt=1,low[maxn],dfn[maxn],scc[maxn],dfscnt,scccnt;//cnt从1开始，标记桥的时候成对标记
bool bridge[maxn],vis[maxn];
vector<int> ans[maxn];
void add(int x,int y)
{
to[++cnt]=y;
nxt[cnt]=head[x];
head[x]=cnt;
}
void tarjan(int u,int fa)
{
low[u]=dfn[u]=++dfscnt;
for(int i=head[u];i;i=nxt[i])
{
int v=to[i];
if(!dfn[v])
{
tarjan(v,i);
low[u]=min(low[u],low[v]);
if(low[v]>dfn[u]) bridge[i]=bridge[i^1]=1;
}
else if(i!=(fa^1)) low[u]=min(low[u],dfn[v]);
}
}
void dfs(int x)
{
scc[x]=scccnt;
if(x) ans[scccnt].push_back(x);
for(int i=head[x];i;i=nxt[i])
{
if(scc[to[i]]||bridge[i]) continue;
dfs(to[i]);
}
}
signed main()
{
ios::sync_with_stdio(false);
cin.tie(0);cout.tie(0);
int n,m;cin>>n>>m;
for(int i=1;i<=m;i++)
{
int u,v;cin>>u>>v;
add(u,v);add(v,u);
}
for(int i=1;i<=n;i++)
if(!dfn[i]) tarjan(i,0);
for(int i=1;i<=n;i++)
if(!scc[i]) ++scccnt,dfs(i);
cout<<scccnt<<endl;
for(int i=1;i<=scccnt;i++)
{
cout<<ans[i].size()<<' ';
for(int j=0;j<ans[i].size();j++) cout<<ans[i][j]<<' ';
cout<<endl;
}
return 0;
}
最小生成树
// 并查集维护是否在一个强连通分量中
#include <bits/stdc++.h>
using namespace std;
int to[500000],v[500000],father[500000],nxt[500000],cnt,head[500000],uu,vv,ans,M,N;
int find(int x)
{
if(father[x]!=x) father[x]=find(father[x]);
return father[x];
}
struct node
{
int to,from,v;
}gg[500000];
bool cmp(node a,node b)
{
return a.v<b.v;
}
int main()
{
cin>>N>>M;
int sum;
for(int i=1;i<=N;i++) father[i]=i;
for(int i=1;i<=M;i++) cin>>gg[i].from>>gg[i].to>>gg[i].v;
sort(gg+1,gg+1+M,cmp);
for(int i=1;i<=N;i++) father[i]=i;
for(int i=1;i<=M&&cnt<N-1;i++)
{
int uu=find(gg[i].from),vv=find(gg[i].to);
if(uu!=vv)
{
sum+=gg[i].v;
father[uu]=vv;
cnt++;
}
}
if(cnt!=N-1)
cout<<"orz";
else
cout<<sum;
return 0;
}
最短路
Dijkstra
#include <bits/stdc++.h>
using namespace std;
const int maxn=2*1e5+5;
int head[maxn],nxt[maxn],to[maxn],w[maxn],cnt,dis[maxn],vis[maxn];
void add(int x,int y,int z)
{
to[++cnt]=y;
w[cnt]=z;
nxt[cnt]=head[x];
head[x]=cnt;
}
struct node
{
int id,dis;
friend bool operator < (node a,node b)
{
return a.dis>b.dis;//小顶堆！
}
};
priority_queue<node> q;
void dij(int s)
{
memset(vis,0,sizeof(vis));
memset(dis,0x3f,sizeof(dis));
dis[s]=0;q.push(node{s,0});
while(!q.empty())
{
int u=q.top().id;q.pop();
if(vis[u]) continue;
vis[u]=1;
for(int i=head[u];i;i=nxt[i])
if(dis[to[i]]>dis[u]+w[i])
dis[to[i]]=dis[u]+w[i],q.push(node{to[i],dis[to[i]]});
}
}
int main()
{
int n,m,s,u,v,w;cin>>n>>m>>s;
for(int i=1;i<=m;i++) cin>>u>>v>>w,add(u,v,w);
dij(s);
for(int i=1;i<=n;i++) cout<<dis[i]<<' ';
return 0;
}
Johnson 全源最短路
#include <bits/stdc++.h>
using namespace std;
#define int long long
const int maxn=9005,maxx=1e9;//注意原有m条边+新建0节点n条边，数组小了会炸
int nxt[maxn],head[maxn],cnt,to[maxn],w[maxn],h[maxn],vis[3005],tim[3005],m,n,u,v,ww,dis[3005];
void add(int x,int y,int z)
{
to[++cnt]=y;
w[cnt]=z;
nxt[cnt]=head[x];
head[x]=cnt;
}
bool spfa()//SPFA判负环
{
queue<int> q;
memset(h,127/3,sizeof(h));
h[0]=0,vis[0]=1;
q.push(0);
while(!q.empty())
{
int u=q.front();q.pop();
vis[u]=0;
for(int i=head[u];i;i=nxt[i])
{
int v=to[i];
if(h[v]>h[u]+w[i])
{
h[v]=h[u]+w[i];
if(!vis[v])
{
q.push(v),vis[v]=1,tim[v]++;
if(tim[v]>n) return true;
}
}
}
}
return false;
}
struct node
{
int id,dis;
bool friend operator < (node a,node b)
{
return a.dis>b.dis;
}
};
void dij(int s)
{
priority_queue<node> q;
memset(vis,0,sizeof(vis));
for(int i=1;i<=n;i++) dis[i]=maxx;
dis[s]=0;q.push(node{s,0});
while(!q.empty())
{
int u=q.top().id;q.pop();
if(vis[u]) continue;
vis[u]=1;
for(int i=head[u];i;i=nxt[i])
if(dis[to[i]]>dis[u]+w[i])
dis[to[i]]=dis[u]+w[i],q.push(node{to[i],dis[to[i]]});
}
}
signed main()
{
cin>>n>>m;
for(int i=1;i<=m;i++) cin>>u>>v>>ww,add(u,v,ww);
for(int i=1;i<=n;i++) add(0,i,0);
if(spfa()) cout<<-1,exit(0);
for(int u=1;u<=n;u++) for(int i=head[u];i;i=nxt[i]) w[i]+=h[u]-h[to[i]];
for(int i=1;i<=n;i++)
{
dij(i);
int ans=0;
for(int j=1;j<=n;j++)
{
if(dis[j]==maxx) ans+=j*maxx;
else ans+=j*(dis[j]+h[j]-h[i]);
}
cout<<ans<<endl;
}
return 0;
}
二分图最大匹配
#include <bits/stdc++.h>
using namespace std;
const int maxn=5e4+5;
int head[maxn],nxt[maxn],to[maxn],vis[maxn],cnt,match[505]/*表示右边点对应左边的cp*/,m,n,e;
void add(int x,int y)
{
to[++cnt]=y;
nxt[cnt]=head[x];
head[x]=cnt;
}
bool dfs(int now)//找对象
{
for(int i=head[now];i;i=nxt[i])
{
if(vis[to[i]]) continue;
vis[to[i]]=1;
if(!match[to[i]]||dfs(match[to[i]]))//想找的对象没有男朋友或者想找的对象还有别的选择
{
match[to[i]]=now;//右侧元素找到了新的匹配
return 1;
}
}
return 0;
}
void gett()
{
int ans=0;
for(int i=1;i<=n;i++)
{
memset(vis,0,sizeof(vis));
if(dfs(i)) ans++;
}
cout<<ans;
}
int main()
{
cin>>n>>m>>e;
for(int i=1;i<=e;i++)
{
int u,v;cin>>u>>v;
add(u,v);
}
gett();
return 0;
}
负权环
#include <bits/stdc++.h>
using namespace std;
const int maxn=40005;
int nxt[maxn],to[maxn],head[maxn],val[maxn],dis[maxn],vis[maxn],rec[maxn],cnt,n,m;
queue<int> q;
void add(int x,int y,int z)
{
to[++cnt]=y;
val[cnt]=z;
nxt[cnt]=head[x];
head[x]=cnt;
}
bool spfa()
{
memset(dis,127,sizeof(dis));
memset(vis,0,sizeof(vis));
memset(rec,0,sizeof(rec));
while(!q.empty) q.pop();
q.push(1);
dis[1]=0,vis[1]=1,rec[1]++;
while(!q.empty())
{
int u=q.front();q.pop();
vis[u]=0;//队首已出队
for(int i=head[u];i;i=nxt[i])
{
if(dis[to[i]]>dis[u]+val[i])
{
dis[to[i]]=dis[u]+val[i];
//对于判断是否有负环，用数组rec记录点的入队次数，如果入队次数>n，就证明出现了负环导致没有最短路
if(!vis[to[i]]) vis[to[i]]=true,rec[to[i]]++,q.push(to[i]);//能更新，压入队列
if(rec[to[i]]>=n) return true;
}
}
}
return false;
}
int main()
{
int T,u,v,w;cin>>T;
while(T--)
{
cin>>n>>m;
memset(head,0,sizeof(head));
cnt=0;
for(int i=1;i<=m;i++)
{
cin>>u>>v>>w;
add(u,v,w);
if(w>=0) add(v,u,w);
}
if(spfa()) cout<<"YES"<<endl;
else cout<<"NO"<<endl;
}
return 0;
}
最大流
#include <bits/stdc++.h>
using namespace std;
#define int long long
typedef long long ll;
const int maxn=5005;
int n,m,s,t,cnt=1,head[maxn],now[maxn];
ll flow[maxn],ans,dis[205];
struct edge{int to,nxt;ll c;}e[maxn*2];
void add(int x,int y,ll z){e[++cnt]={y,head[x],z},head[x]=cnt;}
bool bfs()
{
for(int i=1;i<=n;i++) dis[i]=LLONG_MAX;
queue<int> q;
q.push(s);dis[s]=0,now[s]=head[s];
while(!q.empty())
{
int x=q.front();q.pop();
for(int i=head[x];i;i=e[i].nxt)
if(e[i].c>0&&dis[e[i].to]==LLONG_MAX)
{
q.push(e[i].to),now[e[i].to]=head[e[i].to],dis[e[i].to]=dis[x]+1;
if(e[i].to==t) return 1;
}
}
return 0;
}
int dfs(int x,ll mxf)
{
if(x==t) return mxf;
ll sum=0;
for(int i=now[x];i;i=e[i].nxt)
{
now[x]=i;
if(e[i].c>0&&dis[e[i].to]==dis[x]+1)
{
int ff=dfs(e[i].to,min(mxf-sum,e[i].c));
e[i].c-=ff,e[i^1].c+=ff,sum+=ff;
if(sum==mxf) break;
}
}
return sum;
}
signed main()
{
cin>>n>>m>>s>>t;
for(int i=1,u,v,w;i<=m;i++)
cin>>u>>v>>w,add(u,v,w),add(v,u,0);
while(bfs()) ans+=dfs(s,LLONG_MAX);
cout<<ans;
return 0;
}
费用流
最小费用最大流
#include <bits/stdc++.h>
using namespace std;
const int maxn=5005,maxm=1e5+5;
struct edge{int to,nxt,c,w;}e[maxm];
int head[maxn],pre[maxn],dis[maxn],cnt=1,n,m,s,t,mxf,minc,flow[maxn];
bool vis[maxn];
void add(int x,int y,int z,int q){e[++cnt]=(edge){y,head[x],z,q},head[x]=cnt;}
bool spfa()
{
queue<int> q;
for(int i=1;i<=n;i++) dis[i]=INT_MAX,vis[i]=0;
q.push(s),dis[s]=0,flow[s]=INT_MAX,vis[s]=1,pre[t]=-1;
while(!q.empty())
{
int x=q.front();q.pop(),vis[x]=0;
for(int i=head[x];i;i=e[i].nxt)
if(e[i].c&&dis[e[i].to]>dis[x]+e[i].w)
{
dis[e[i].to]=dis[x]+e[i].w,pre[e[i].to]=i,flow[e[i].to]=min(flow[x],e[i].c);
if(!vis[e[i].to]) q.push(e[i].to),vis[e[i].to]=1;
}
}
return pre[t]!=-1;
}
void ek()
{
while(spfa())
{
int x=t;
mxf+=flow[t],minc+=flow[t]*dis[t];
while(x!=s) e[pre[x]].c-=flow[t],e[pre[x]^1].c+=flow[t],x=e[pre[x]^1].to;
// cout<<flow[t]<<' '<<dis[t]<<endl;
}
}
int main()
{
cin>>n>>m>>s>>t;
for(int i=1,u,v,w,c;i<=m;i++) cin>>u>>v>>w>>c,add(u,v,w,c),add(v,u,0,-c);
ek();
cout<<mxf<<' '<<minc;
return 0;
}
差分约束
#include <bits/stdc++.h>
using namespace std;
const int maxn=1e4+5;
int head[maxn],tot[maxn],n,m,dis[maxn],cnt;
bool vis[maxn];
struct edge{int to,nxt,w;}e[maxn];
void add(int x,int y,int z){e[++cnt]={y,head[x],z},head[x]=cnt;}
bool spfa(int s)
{
memset(dis,0x3f3f3f,sizeof dis);
queue<int> q;
dis[s]=0,vis[s]=1,++tot[s],q.push(s);
while(!q.empty())
{
int x=q.front();q.pop(),vis[x]=0;
for(int i=head[x];i;i=e[i].nxt)
if(dis[e[i].to]>dis[x]+e[i].w)
{
dis[e[i].to]=dis[x]+e[i].w;
if(!vis[e[i].to])
{
q.push(e[i].to),vis[e[i].to]=1,++tot[e[i].to];
if(tot[e[i].to]==n+1) return 1;
}
if(tot[e[i].to]==n+1) return 1;
}
}
return 0;
}
int main()
{
cin>>n>>m;
for(int i=1;i<=n;i++) add(0,i,0);
for(int i=1,u,v,y;i<=m;i++) cin>>u>>v>>y,add(v,u,y);
if(spfa(0)) cout<<"NO";
else for(int i=1;i<=n;++i) cout<<dis[i]<<' ';
return 0;
}
数学
裴蜀定理
#include <bits/stdc++.h>
using namespace std;
int n,a,ans;
int gcd(int x,int y)
{
return !y?x:gcd(y,x%y);
}
int main()
{
cin>>n;
for(int i=1;i<=n;i++)
cin>>a,a=a>0?a:-a,ans=gcd(ans,a);
cout<<ans;
return 0;
}
乘法逆元
#include<bits/stdc++.h>
using namespace std;
#define int long long
int n,p,a[5000005];
signed main()
{
cin>>n>>p;
a[1]=1;
printf("1\n");
for(int i=2;i<=n;i++)
a[i]=(p-p/i)*a[p%i]%p,cout<<a[i]<<endl;
return 0;
}
线性筛
这里采用的是我惯用的写法。
#include <bits/stdc++.h>
using namespace std;
const int maxn=1e8+5;
int v[maxn],prime[maxn];
void primes(int n)
{
memset(v,0,sizeof(v));//存储最小质因子
int m=0;//质数数量
for(int i=2;i<=n;i++)
{
if(v[i]==0)
{
v[i]=i;
prime[++m]=i;
}
for(int j=1;j<=m;j++)
{
//i有比prime[j]更小的质因子或者超出n的范围
if(prime[j]>v[i]||prime[j]*i>n) break;
v[i*prime[j]]=prime[j];
}
}
}
int main()
{
int n,q;cin>>n>>q;
primes(n);
for(int i=1;i<=q;i++)
{
int k;cin>>k;
cout<<prime[k]<<endl;
}
return 0;
}
快速幂
#include <bits/stdc++.h>
using namespace std;
#define int long long
int quickpow(int a,int b,int n)
{
if(b==1) return a;
if(b%2==0)
{
int t=quickpow(a,b/2,n);
return t*t%n;
}
else
{
int t=quickpow(a,b/2,n);
t=t*t%n;
t=t*a%n;
return t;
}
}
signed main()
{
int a,b,p;
cin>>a>>b>>p;
printf("%lld^%lld mod %lld=%lld",a,b,p,quickpow(a,b,p));
}
线性数据结构
单调栈
#include <bits/stdc++.h>
using namespace std;
int a[3000005],f[3000005];
stack<int> s;
int main()
{
int n;cin>>n;
for(int i=1;i<=n;i++) cin>>a[i];
for(int i=n;i>=1;i--)
{
while(!s.empty()&&a[s.top()]<=a[i]) s.pop();
f[i]=s.empty()?0:s.top();
s.push(i);
}
for(int i=1;i<=n;i++) cout<<f[i]<<' ';
return 0;
}
单调队列
#include <bits/stdc++.h>
using namespace std;
#define int long long
const int maxn=1e6+5;
int q1[maxn],q2[maxn],a[maxn];
signed main()
{
int n,k;cin>>n>>k;
for(int i=1;i<=n;i++) cin>>a[i];
int h1=1,t1=0,h2=1,t2=0;
for(int i=1;i<=n;i++)
{
while(h1<=t1&&q1[h1]+k<=i) h1++;
while(h1<=t1&&a[q1[t1]]>a[i]) t1--;
q1[++t1]=i;
if(i>=k) cout<<a[q1[h1]]<<' ';
}
cout<<endl;
for(int i=1;i<=n;i++)
{
while(h2<=t2&&q2[h2]+k<=i) h2++;
while(h2<=t2&&a[q2[t2]]<a[i]) t2--;
q2[++t2]=i;
if(i>=k) cout<<a[q2[h2]]<<' ';
}
return 0;
}
ST 表
#include <bits/stdc++.h>
using namespace std;
const int maxn=1e5+5,maxm=2e6+5;
int a[maxn],lg[maxn],f[maxn][25];
inline int read()
{
int x=0,f=1;char ch=getchar();
while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
while (ch>='0'&&ch<='9'){x=x*10+ch-48;ch=getchar();}
return x*f;
}
void query(int l,int r)
{
int k=lg[r-l+1];
cout<<max(f[l][k],f[r-(1<<k)+1][k])<<'\n';
}
int main()
{
ios::sync_with_stdio(false);
int M,N;cin>>N>>M;
lg[1]=0;
for(int i=1;i<=N;i++) f[i][0]=read();
for(int i=2;i<=N;i++) lg[i]=lg[i/2]+1;
for(int j=1;j<=lg[N];j++)
for(int i=1;i<=N-(1<<j)+1;i++) f[i][j]=max(f[i][j-1],f[i+(1<<(j-1))][j-1]);
while(M--)
{
int l=read(),r=read();
query(l,r);
}
return 0;
}
树形数据结构
并查集
#include<bits/stdc++.h>
using namespace std;
int i,j,k,n,m,ans,f[10005],z,x,y;
int find(int k)
{
if(f[k]==k) return k;
return f[k]=find(f[k]);
}
int main()
{
cin>>n>>m;
for(i=1;i<=n;i++)
f[i]=i;
for(i=1;i<=m;i++)
{
cin>>z>>x>>y;
if(z==1)
f[find(x)]=find(y);
else if(find(x)==find(y))
printf("Y\n");
else
printf("N\n");
}
return 0;
}
树状数组
#include <bits/stdc++.h>
using namespace std;
const int maxn=5*1e5+5;
int c[maxn],a[maxn],n;
int lowbit(int x)
{
return x&(-x);
}
int ask(int x)//查询x的前缀和
{
int res=0;
while(x) res+=c[x],x-=lowbit(x);
return res;
}
void add(int x,int v)//把点x增加v
{
while(x<=n) c[x]+=v,x+=lowbit(x);
}
int main()
{
int qwq,m;cin>>n>>m;
for(int i=1;i<=n;i++) cin>>qwq,add(i,qwq);
for(int i=1;i<=m;i++)
{
int opt,x,y,k;
cin>>opt;
if(opt==1)
cin>>x>>k,add(x,k);
else
cin>>x>>y,cout<<(ask(y)-ask(x-1))<<endl;
//这里x要-1，可以考虑一下如果是要求[1,2]的区间和，应该用1到2的区间和减去ask(0)即0
}
return 0;
}
#include <bits/stdc++.h>
using namespace std;
const int maxn=5*1e5+5;
int c[maxn],N;
int lowbit(int x)
{
return x&(-x);
}
int ask(int x)//查询x的前缀和
{
int res=0;
while(x) res+=c[x],x-=lowbit(x);
return res;
}
void add(int x,int v)//把点x增加v
{
while(x<=N) c[x]+=v,x+=lowbit(x);
}
int main()
{
int M,now,last;cin>>N>>M;
for(int i=1;i<=N;i++) cin>>now,add(i,now-last),last=now;
for(int i=1;i<=M;i++)
{
int opt,x,y,k;cin>>opt;
if(opt==1) cin>>x>>y>>k,add(x,k),add(y+1,-k);//差分
else cin>>x,cout<<ask(x)<<endl;
}
return 0;
}
线段树
#include<bits/stdc++.h>
using namespace std;
#define int long long
const int N=1e5+5;
int n,m,a[N],s[N<<2],lazy[N<<2];
inline void buildtree(int k,int l,int r)
{
if(l>r) return;
if(l==r){s[k]=a[l];return;}
int mid=(l+r)>>1;
buildtree(k<<1,l,mid);
buildtree(k<<1|1,mid+1,r);
s[k]=s[k<<1]+s[k<<1|1];
}
inline void add(int k,int l,int r,int val)
{
s[k]+=(r-l+1)*val;
lazy[k]+=val;
}
inline void pushdown(int k,int l,int r)
{
int mid=(l+r)>>1;
add(k<<1,l,mid,lazy[k]);
add(k<<1|1,mid+1,r,lazy[k]);
lazy[k]=0;
}
inline void pushup(int k,int l,int r)
{
s[k]=s[k<<1]+s[k<<1|1];
}
inline void plu(int k,int l,int r,int x,int y,int val)
{
if(x<=l&&y>=r) return add(k,l,r,val);
pushdown(k,l,r);
int mid=l+r>>1;
if(x<=mid) plu(k<<1,l,mid,x,y,val);
if(y>mid) plu(k<<1|1,mid+1,r,x,y,val);
pushup(k,l,r);
}
inline int ask(int k,int l,int r,int x,int y)
{
if(x<=l&&y>=r) return s[k];
pushdown(k,l,r);
int mid=(l+r)>>1,res=0;
if(x<=mid) res+=ask(k<<1,l,mid,x,y);
if(y>mid) res+=ask(k<<1|1,mid+1,r,x,y);
return res;
}
signed main()
{
cin>>n>>m;
for(int i=1;i<=n;i++) cin>>a[i];
buildtree(1,1,n);
for(int i=1;i<=m;i++)
{
int opt,x,y,k;cin>>opt>>x>>y;
if(opt==1)
{
cin>>k;
plu(1,1,n,x,y,k);
}
else cout<<ask(1,1,n,x,y)<<endl;
}
return 0;
}
#include <bits/stdc++.h>
using namespace std;
#define int long long
const int maxn=4e5+5;
int a[maxn],s[maxn],laz1[maxn],laz2[maxn],m,k;
void build(int k,int l,int r)
{
laz1[k]=1;
if(l==r){s[k]=a[l]%m;return;}
int mid=(l+r)/2;
build(k*2,l,mid),build(k*2+1,mid+1,r),s[k]=(s[k*2]+s[k*2+1])%m;
}
void add(int k,int l,int r,int q1,int q2)
{
s[k]=(s[k]*q1%m+(r-l+1)*q2%m)%m;
laz1[k]=(laz1[k]*q1)%m,laz2[k]=(laz2[k]*q1+q2)%m;
}
void pushup(int k){s[k]=(s[k*2]+s[k*2+1])%m;}
void pushdown(int k,int l,int r)
{
int mid=(l+r)/2;
add(k*2,l,mid,laz1[k],laz2[k]),add(k*2+1,mid+1,r,laz1[k],laz2[k]);
laz1[k]=1,laz2[k]=0;
}
void update(int k,int l,int r,int x,int y,int laz1,int laz2)
{
if(x<=l&&y>=r) return add(k,l,r,laz1,laz2);
pushdown(k,l,r);
int mid=(l+r)/2;
if(x<=mid) update(k*2,l,mid,x,y,laz1,laz2);
if(y>mid) update(k*2+1,mid+1,r,x,y,laz1,laz2);
pushup(k);
}
int query(int k,int l,int r,int x,int y)
{
if(x<=l&&y>=r) return s[k];
int res=0,mid=(l+r)/2;
pushdown(k,l,r);
if(x<=mid) res=(res+query(k*2,l,mid,x,y))%m;
if(y>mid) res=(res+query(k*2+1,mid+1,r,x,y))%m;
return res;
}
signed main()
{
int n,q;cin>>n>>q>>m;
for(int i=1;i<=n;i++) cin>>a[i];
build(1,1,n);
for(int i=1,op,x,y;i<=q;i++)
{
cin>>op>>x>>y;
if(op==1) cin>>k,update(1,1,n,x,y,k,0);
else if(op==2) cin>>k,update(1,1,n,x,y,1,k);
else cout<<query(1,1,n,x,y)<<endl;
}
return 0;
}
可持久化线段树
#include <bits/stdc++.h>
using namespace std;
const int maxn=4e7+5;
int N,M,a[maxn],cnt,rr[maxn]/*存储版本i的根节点编号从0开始*/;
struct node
{
int ls,rs,v;
}t[maxn];
int build(int l,int r)
{
int id=++cnt;
if(l==r) t[id].v=a[l];
else
{
int mid=(l+r)/2;
t[id].ls=build(l,mid);
t[id].rs=build(mid+1,r);
}
return id;
}//记录左右儿子编号
int neww(int x)
{
cnt++,t[cnt]=t[x];
return cnt;
}//新建节点
int update(int now,int l,int r,int k,int v)
{
now=neww(now);
if(l==r) t[now].v=v;
else
{
int mid=(l+r)/2;
if(k<=mid) t[now].ls=update(t[now].ls,l,mid,k,v);
else t[now].rs=update(t[now].rs,mid+1,r,k,v);
}
return now;
}
int query(int now,int k,int l,int r)
{
if(l==r) return t[now].v;
else
{
int mid=(l+r)/2;
if(k<=mid) return query(t[now].ls,k,l,mid);
else return query(t[now].rs,k,mid+1,r);
}
}
signed main()
{
ios::sync_with_stdio(false);
cin.tie(0);cout.tie(0);
cin>>N>>M;
for(int i=1;i<=N;i++) cin>>a[i];
rr[0]=build(1,N);
for(int i=1;i<=M;i++)
{
int opt,v,loc,value;
cin>>v>>opt>>loc;
if(opt==1) cin>>value,rr[i]=update(rr[v],1,N,loc,value);
else cout<<query(rr[v],loc,1,N)<<endl,rr[i]=rr[v];
}
return 0;
}
#include <bits/stdc++.h>
using namespace std;
#define int long long
const int maxn=6e6+5;
int n,m,rr[maxn],s[maxn],cnt,a[maxn],b[maxn];
struct node
{
int ls,rs;
}t[maxn];
int build(int l,int r)
{
int id=++cnt;
s[id]=0;
if(l==r) return id;
int mid=(l+r)/2;
t[id].ls=build(l,mid);
t[id].rs=build(mid+1,r);
return id;
}//记录左右儿子编号
int update(int now,int l,int r,int k)
{
int rt=++cnt;
t[rt].ls=t[now].ls,t[rt].rs=t[now].rs;
s[rt]=s[now]+1;
if(l<r)
{
int mid=(l+r)/2;
if(k<=mid) t[rt].ls=update(t[now].ls,l,mid,k);
else t[rt].rs=update(t[now].rs,mid+1,r,k);
}
return rt;
}
int query(int u,int v,int k,int l,int r)
{
if(l>=r) return l;
else
{
int mid=(l+r)/2,x=s[t[v].ls]-s[t[u].ls];
if(x>=k) return query(t[u].ls,t[v].ls,k,l,mid);
else return query(t[u].rs,t[v].rs,k-x,mid+1,r);
}
}
signed main()
{
cin>>n>>m;
for(int i=1;i<=n;i++) cin>>a[i],b[i]=a[i];
sort(b+1,b+n+1);
int mm=unique(b+1,b+n+1)-b-1;
rr[0]=build(1,mm);
for(int i=1;i<=n;i++)
a[i]=lower_bound(b+1,b+mm+1,a[i])-b,rr[i]=update(rr[i-1],1,mm,a[i]);
for(int i=1;i<=m;i++)
{
int l,r,k;cin>>l>>r>>k;
int ans=query(rr[l-1],rr[r],k,1,mm);
cout<<b[ans]<<endl;
}
return 0;
}
三维偏序
#include <bits/stdc++.h>
using namespace std;
#define int long long
const int maxn=2e5+5;
int c[maxn],qwq[maxn],n,k;
int lowbit(int x){return x&(-x);}
struct node
{
int a,b,c,cnt,ans;
}p1[maxn],p2[maxn];
//---OK---
bool cmp1(node aa,node bb)
{
if(aa.a!=bb.a) return aa.a<bb.a;
else if(aa.b!=bb.b) return aa.b<bb.b;
return aa.c<bb.c;
}
bool cmp2(node aa,node bb)
{
if(aa.b!=bb.b) return aa.b<bb.b;
return aa.c<bb.c;
}
//---------
void add(int x,int v)
{
while(x<=k) c[x]+=v,x+=lowbit(x);
return;
}
int ask(int x)
{
int res=0;
while(x) res+=c[x],x-=lowbit(x);
return res;
}
void cdq(int l,int r)
{
if(l==r) return;
int mid=(l+r)/2;
cdq(l,mid);cdq(mid+1,r);
sort(p2+l,p2+mid+1,cmp2);sort(p2+mid+1,p2+r+1,cmp2);
int j=l;
for(int i=mid+1;i<=r;i++)
{
while(p2[i].b>=p2[j].b&&j<=mid) add(p2[j].c,p2[j].cnt),j++;
p2[i].ans+=ask(p2[i].c);
}
for(int i=l;i<j;i++) add(p2[i].c,-p2[i].cnt);
}
signed main()
{
cin>>n>>k;
for(int i=1;i<=n;i++) cin>>p1[i].a>>p1[i].b>>p1[i].c;
sort(p1+1,p1+n+1,cmp1);
int cntt=0,sum=0;
for(int i=1;i<=n;i++)
{
cntt++;
if(p1[i].a!=p1[i+1].a||p1[i].b!=p1[i+1].b||p1[i].c!=p1[i+1].c)
sum++,p2[sum].a=p1[i].a,p2[sum].b=p1[i].b,p2[sum].c=p1[i].c,p2[sum].cnt=cntt,cntt=0;
}
cdq(1,sum);
// for(int i=1;i<=sum;i++) cout<<p2[i].ans<<p2[i].a<<p2[i].b<<p2[i].c<<endl;
for(int i=1;i<=sum;i++) qwq[p2[i].ans+p2[i].cnt-1]+=p2[i].cnt;
for(int i=0;i<n;i++) cout<<qwq[i]<<endl;
return 0;
}
字符串
Hash
#include <bits/stdc++.h>
using namespace std;
typedef unsigned long long ull;
const int mmax=1505,maxn=10005;
ull base=131,prime=23317,mod=212370440130137957;
int N,a[maxn],ans=1;
char s[mmax];
ull hash[maxn],power[maxn];//自然溢出
ull hashh(char s[])//先构造哈希函数，递推求解
{
int len=strlen(s);
ull ans=0;
for(int i=0;i<len;i++)
ans=(ans*base+(ull)s[i])%mod+prime;
return ans;
}
int main()
{
cin>>N;
for(int i=1;i<=N;i++)
scanf("%s",s),a[i]=hashh(s);
//一个小技巧，统计一组数中不同数的个数，可以先排序然后统计前一个数和后一个数是否不相等，若不相等ans++
sort(a+1,a+N+1);
for(int i=1;i<N;i++)
if(a[i]!=a[i+1]) ans++;
cout<<ans;
return 0;
}
字典树
#include <bits/stdc++.h>
using namespace std;
const int maxn=2e6+5;
int t[maxn][65],cnt[maxn],tot;
char s[maxn];
int getn(char x)
{
if(x<='Z'&&x>='A') return x-'A';
else if(x<='z'&&x>='a') return x-'a'+26;
else return x-'0'+52;
}
void insert(char s[])
{
int p=0,len=strlen(s);//根节点为0
for(int i=0;i<len;i++)
{
int u=getn(s[i]);//当前字母子节点
if(!t[p][u]) t[p][u]=++tot;//如果当前子节点不存在就创造一个点来存储子节点
p=t[p][u];//让p走到子节点的位置
cnt[p]++;//结尾是p的字符串个数增加
}
}
int ask(char s[])
{
int p=0,len=strlen(s);
for(int i=0;i<len;i++)
{
int u=getn(s[i]);
if(!t[p][u]) return 0;
p=t[p][u];
}
return cnt[p];
}
int main()
{
int T;cin>>T;
while(T--)
{
for(int i=0;i<tot;i++)
for(int j=0;j<65;j++) t[i][j]=0;
for(int i=0;i<tot;i++) cnt[i]=0;
tot=0;
int n,q;cin>>n>>q;
for(int i=1;i<=n;i++)
cin>>s,insert(s);
while(q--)
cin>>s,cout<<ask(s)<<endl;
}
return 0;
}
KMP
#include <bits/stdc++.h>
using namespace std;
const int maxn=1e6+5;
char a[maxn],b[maxn];
int nxt[maxn],n,m;
int main()
{
scanf("%s%s",a+1,b+1);
n=strlen(a+1),m=strlen(b+1);
for(int i=1,j=0;i<=m;i++)
{
while(j&&b[i+1]!=b[j+1]) j=nxt[j];
if(b[i+1]==b[j+1]) j++;
nxt[i+1]=j;
}
for(int i=0,j=0;i<=n;i++)
{
while(j&&a[i+1]!=b[j+1]) j=nxt[j];
if(a[i+1]==b[j+1]) j++;//如果当前位能匹配
if(j==m) cout<<(i+2-m)<<endl;
}
for(int i=1;i<=m;i++) cout<<nxt[i]<<" ";
return 0;
}
Z 函数
#include <bits/stdc++.h>
using namespace std;
const int maxn=2*1e7+5;
char a[maxn],b[maxn];
long long za[maxn],zb[maxn];
int main()
{
scanf("%s%s",a,b);
int la=strlen(a),lb=strlen(b);
zb[0]=lb;
for(int i=1,l=0,r=0;i<lb;++i)
{
if(i<=r&&zb[i-l]<r-i+1) zb[i]=zb[i-l];//因为不确定下一位是否相等，所以直接继承z[i-l]
else
{
zb[i]=max(0,r-i+1);//0说明i>r，r-i+1说明超出匹配区间
while(i+zb[i]<lb&&b[zb[i]]==b[i+zb[i]]) ++zb[i];
}
if(i+zb[i]-1>r) l=i,r=i+zb[i]-1;
}
for(int i=0;i<la;++i)//求a关于模板串b的z[0]
{
if(b[i]==a[i]) ++za[0];
else break;
}
for(int i=1,l=0,r=0;i<la;++i)
{
if(i<=r&&zb[i-l]<r-i+1) za[i]=zb[i-l];
else
{
za[i]=max(0,r-i+1);
while(i+za[i]<la&&b[za[i]]==a[i+za[i]]) ++za[i];
}
if(i+za[i]-1>r) l=i,r=i+za[i]-1;
}
long long ans=0,bp=0;
for(int i=0;i<la;i++) ans^=(i+1)*(za[i]+1);
for(int i=0;i<lb;i++) bp^=(i+1)*(zb[i]+1);
cout<<bp<<endl<<ans;
return 0;
}
Manacher 算法
#include <bits/stdc++.h>
using namespace std;
const int maxn=4.2*1e7+5;
char s[maxn],ss[maxn];
int p[maxn];
int init()
{
int len=strlen(ss+1);
int j=3;
s[1]='#',s[2]='$';
for(int i=1;i<=len;i++)
s[j++]=ss[i],s[j++]='$';
s[j]='@';
return j;
}
int manacher()
{
int len=init(),mid=1,mx=1,ans=1;
for(int i=1;i<=len;i++)
{
if(i<mx) p[i]=min(mx-i,p[mid*2-i]);
else p[i]=1;//前面的mx一点儿用没用，直接赋为1开始暴力
while(s[i-p[i]]==s[i+p[i]]) p[i]++;
if(mx<i+p[i]) mid=i,mx=i+p[i];//更新mid和mx
ans=max(ans,p[i]-1);
}
return ans;
}
int main()
{
scanf("%s",ss+1);
cout<<manacher();
return 0;
}
DP
LCS
#include <bits/stdc++.h>
using namespace std;
int a[100005],b[100005],c[100005],dp[100005],mp[100005];
int main()
{
int n;cin>>n;
for(int i=1;i<=n;i++)
cin>>a[i],mp[a[i]]=i;
for(int i=1;i<=n;i++) cin>>b[i];
for(int i=1;i<=n;i++)
c[i]=mp[b[i]];
int tmp=0,ans=-0x7fffffff;
for(int i=1;i<=n;i++)
{
int j=upper_bound(dp+1,dp+tmp+1,c[i])-dp;
if(j==tmp+1) tmp++,dp[j]=c[i];
else dp[j]=c[i];
}
cout<<tmp;
}
线性代数
线性基
#include <bits/stdc++.h>
using namespace std;
#define int long long
int p[55];
void insert(int x)
{
for(int i=50;i>=0;--i)
if(x>>i&1)
{
if(!p[i]) {p[i]=x;break;}
else x^=p[i];
}
}
int pmax()
{
int ans=0;
for(int i=50;i>=0;--i)
if((ans^p[i])>ans) ans^=p[i];
return ans;
}
signed main()
{
int n,x;cin>>n;
for(int i=1;i<=n;i++) cin>>x,insert(x);
cout<<pmax();
return 0;
}
基础算法
三分
#include <bits/stdc++.h>
using namespace std;
const double eps=1e-6;
double a[15],l,r;
int N;
double f(double x)
{
double res=0;
for(int i=N;i>=0;i--) res=res*x+a[i];
return res;
}
int main()
{
cin>>N>>l>>r;
for(int i=N;i>=0;i--) cin>>a[i];
while(r-l>eps)
{
double mid=(l+r)/2.0;
if(f(mid-eps)>f(mid)) r=mid;
else l=mid;
}
printf("%.5lf",l);
return 0;
}
其他
整数快读
inline int read()
{
int x=0,f=1;
char ch=getchar();
while(ch>'9'||ch<'0'){if(ch=='-') f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
return x*f;
}
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
