模板汇总

%%会不断更新……

1.网络流

（1）最大流最小割模板

#include <iostream>
#include <queue>
#include <cstring>
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;
const int inf=0x7fffffff;
const int N=3e5+5;
const int E=5e6+5;
typedef struct{
    int from,to,c,next;
}edge;
edge e[E];
int ans,cnt=1;
int cur[N],dist[N],head[N],b[505][505],c[505];
//最大流=最小割(定理)
//基于dinic算法的最大流(最小割)求解模板开始/*
//插边函数
inline void add(int u,int v,int c){
    e[++cnt]=(edge){u,v,c,head[u]};head[u]=cnt;
    e[++cnt]=(edge){v,u,0,head[v]};head[v]=cnt;
};
//构建层次网络
inline bool bfs(int st,int ed)
{
    mem(dist,-1);
    queue<int> q;
    q.push(0);
    dist[0] = 1;
    while (!q.empty())
    {
        int vTop = q.front();q.pop();
        for (int i = head[vTop] ; i; i = e[i].next)
        {
            //目标节点还未被确立层次，并且可以到达目标节点，就进行层次更新
            if (e[i].c>0 && dist[e[i].to]==-1)
            {
                dist[e[i].to] = dist[vTop] + 1;
               if(e[i].to==ed) return true;
                q.push(e[i].to);
            }
        }
    }
    return false;
}
//dfs查找所有增广路并做流量调整
inline int dfs(int u ,int end , int delta)
{
    if (u == end||!delta) return delta;
    int res = 0;
    for (int &i = cur[u] ; i   ; i = e[i].next)
    {
        if (e[i].c>0&&dist[e[i].to] == dist[u] + 1)
        {
            int dd = dfs(end, e[i].to,e[i].c<(delta-res)? e[i].c:(delta-res));
            e[i].c  -= dd;	//前向弧流量增加(capacity为可调整流量)
            e[i ^ 1].c += dd; //反向弧流量减少
            delta -= dd;
            res += dd;
            if (res==delta) break;
        }
    }
    return res;
}

inline int dinic(int st,int end)
{
    int i,ret = 0;
    while (bfs(st,end))
    {
        for(i=st;i<=end;i++) cur[i]=head[i];
        ret += dfs(st, end, inf);
    }
    return ret;
}
//模板结束 */
//开始解题
int main(){
    ios::sync_with_stdio(0);
    cin.tie(0);cout.tie(0);
    //输入数据，建图……
    /*
    int n,i,j,sum=0,nn,s=0,t;
    cin>>n;
    nn=n*n;t=nn+n+1;
    for(i=1;i<=n;i++){
        for(j=1;j<=n;j++){
           cin>>b[i][j];
           sum+=b[i][j];
        }
    }
    for(i=1;i<=n;i++){
        cin>>c[i];
    }
    //插入边
    for(i=1;i<=n;i++){
        for(j=1;j<=n;j++){
            add(s,n*(i-1)+j,b[i][j]);
            add(n*(i-1)+j,nn+i,inf);
            add(n*(i-1)+j,nn+j,inf);
        }
    }
    for(i=1;i<=n;i++){
        add(nn+i,t,c[i]);
    }
    //求最小割，ans=ans-最小割；
    sum-=dinic(s,t);
    cout<<sum<<'\n';
    */
    return 0;
}

（2）最小费用模板

#include <iostream>
#include <vector>
#include <cstdio>
#include <cstring>
#define inf 0x7fffffff
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;
const int N=1e5+5;
const int MAXN=8e7+5;
int n,n2,ans,cnt=0,st,ed;
int head,tail;
int list[N],pre[N],vi[N],first[N],c[N];
bool v[N];
typedef struct {//节点类型
    int from,to,next,f,w,other;
}Node;//f为流量限制，w为花费
typedef struct {
    int x,y;
}point;
point p1[N],p2[N];
Node e[MAXN];
//模板开始/*
inline void insert(int u,int v,int f,int w){//插入算法
    int kx,ky;
    cnt++;
    kx=cnt;
    e[cnt].from=u;
    e[cnt].to=v;
    e[cnt].f=f;
    e[cnt].w=w;
    e[cnt].next=first[u];
    first[u]=cnt;
    cnt++;ky=cnt;
    e[cnt].from=v;
    e[cnt].to=u;
    e[cnt].f=0;
    e[cnt].w=-w;
    e[cnt].next=first[v];
    first[v]=cnt;
    e[kx].other=ky;
    e[ky].other=kx;
}
inline bool spfa(){//求解算法
    mem(c,0);c[st]=999999999;
    //memset(v,false,sizeof(v));
    mem(v,false);
    v[st]=true;
    mem(vi,0x3F);vi[st]=0;
    list[1]=st;head=1;tail=2;pre[st]=0;
    while(head!=tail){
        int x=list[head];
        for(int k=first[x];k>0;k=e[k].next){
            int y=e[k].to;
            if(e[k].f>0&&vi[y]>vi[x]+e[k].w){
                vi[y]=vi[x]+e[k].w;
                c[y]=min(c[x],e[k].f);
                pre[y]=k;
                if(v[y]==false){
                    v[y]=true;
                    list[tail++]=y;
                    if(tail==ed+1)tail=1;
                }
            }
        }
        head++;if(head==ed+1)head=1;
        v[x]=false;
    }
    if(vi[ed]>999999999)return false;
    ans+=vi[ed]*c[ed];
    int x=ed;
    while(x!=st){
        int k=pre[x];
        e[k].f-=c[ed];e[e[k].other].f+=c[ed];
        x=e[k].from;
    }
    return true;
}
inline int dis(point p1,point p2){
    if(p2.x<p1.x||p2.y>p1.y) return 0;
    else return p2.x-p1.x+p1.y-p2.y;
}
//模板结束*/
int main(){
    ios::sync_with_stdio(0);
    cin.tie(0);cout.tie(0);
    //输入数据，建图……
    /*
    cin>>n;
    st=2*n+1;ed=2*n+2;
    cnt=0;mem(first,0);
    int i,j,k;
    for(i=1;i<=n;i++){
        cin>>p1[i].x>>p1[i].y;
    }
    for(i=1;i<=n;i++){
       cin>>p2[i].x>>p2[i].y;
    }
    n2=2*n;
    //插入节点
    for(i=1;i<=n;i++) insert(st,i,1,0);//源点S到每个天然气井
    for(j=n+1;j<=n2;j++) insert(j,ed,1,0);//每个中转站到汇点T
    //每个天然气井到每个中转站
    for(i=1;i<=n;i++){
        for(j=1;j<=n;j++){
            k=dis(p1[i],p2[j]);
            if(k) insert(i,n+j,1,k);
        }
    }
    ans=0;
    while(spfa());
    cout<<ans<<'\n';*/
    return 0;
}