dinic模板:
struct edge{
int from,to,w,next;
}e[N<<1];
int cnt;
int dep[1111],head[2111];
//queue<int> que;
void init(){
memset(head,-1,sizeof(head));
cnt=0;
}
void add(int u,int v,int w)
{
e[cnt].from=u;
e[cnt].to=v;
e[cnt].w=w;
e[cnt].next=head[u];
head[u]=cnt++;
e[cnt].from=v;
e[cnt].to=u;
e[cnt].w=0;
e[cnt].next=head[v];
head[v]=cnt++;
}
int bfs(int s,int t)
{
// while(!que.empty()) que.pop();
queue<int> que;
memset(dep,-1,sizeof(dep));
dep[s]=0;
que.push(s);
while(!que.empty()){
int u=que.front();
que.pop();
// if(u==t) return 1;
for(int i=head[u];i!=-1;i=e[i].next){
int w=e[i].w;
int v=e[i].to;
if(w&&dep[v]==-1){
dep[v]=dep[u]+1;
que.push(v);
if(v==t) return 1;
}
}
}
return 0;
}
int dfs(int u,int mi,int t){
if(u==t||mi==0) return mi;
int tmp;
int ret=0;
// for(int v=1;v<=t;v++){
for(int i=head[u];i!=-1;i=e[i].next){
int w=e[i].w;
int v=e[i].to;
if(w&&dep[v]==dep[u]+1&&mi-ret>0){ ///mi-ret>0少一次递归,优化
tmp=dfs(v,min(mi-ret,w),t) ;
if(!tmp) dep[v]=0;
e[i].w-=tmp;
e[i^1].w+=tmp;
ret+=tmp;
if(ret==mi) return ret;
}
}
return ret;
}
int dinic(int s,int t)
{
int ans=0;
while(bfs(s,t)){
ans+=dfs(s,INF,t);
}
return ans;
}
最小费用最大流模板:
const int N=100;
const int INF=0x3f3f3f3f;
struct edge{
int u,v,flow,cost,next;
edge(){}
edge(int _u,int _v,int _flow,int _cost,int _next){
u=_u;v=_v;flow=_flow;cost=_cost;next=_next;
}
}e[30000];
int head[N],cnt;
queue<int> que;
void init()
{
memset(head,-1,sizeof(head));
cnt=0;
}
void add(int x,int y,int z,int c)
{
e[cnt]=edge(x,y,z,c,head[x]);
head[x]=cnt++;
e[cnt]=edge(y,x,0,-c,head[y]);
head[y]=cnt++;
}
bool vis[N];
int dis[N],pre[N];
bool spfa(int s,int t)
{
memset(dis,INF,sizeof(dis));
memset(vis,0,sizeof(vis));
memset(pre,-1,sizeof(pre));
dis[s]=0;
vis[s]=true;
while(!que.empty()) que.pop();
que.push(s);
while(!que.empty())
{
int x=que.front();
que.pop();
vis[x]=false;
for(int i=head[x];i!=-1;i=e[i].next){
if(e[i].flow){
int y=e[i].v;
if(dis[x]+e[i].cost<dis[y]){
dis[y]=dis[x]+e[i].cost;
pre[y]=i;
if(!vis[y]){
que.push(y);
vis[y]=true;
}
}
}
}
}
return pre[t]!=-1;
}
int costflow(int s,int t,int x)///源点x,汇点t,需要的总流x
{
int ret=0;
int max_flow=0;
while(spfa(s,t))
{
int min_flow=INF;
for(int i=pre[t];i!=-1;i=pre[e[i].u]){
min_flow=min(e[i].flow,min_flow);
}
for(int i=pre[t];i!=-1;i=pre[e[i].u]){
e[i].flow-=min_flow;
e[i^1].flow+=min_flow;
}
ret+=dis[t]*min_flow;
max_flow+=min_flow;
}
if(max_flow<x) return -1;///达不到满流
return ret;
}