ZKW费用流

bool spfa()
{
for(int i=0;i<=T;i++)dis[i]=inf;
int head=0,tail=1;
dis[0]=0;q[0]=0;inq[0]=1;
while(head!=tail)
{
int now=q[head++];if(head==1601)head=0;
for(int i=last[now];i;i=e[i].next)
if(e[i].v&&e[i].c+dis[now]<dis[e[i].to])
{
dis[e[i].to]=e[i].c+dis[now];
from[e[i].to]=i;
if(!inq[e[i].to])
{
inq[e[i].to]=1;
if(dis[e[i].to]<dis[q[head]])
{
head--;if(head==-1)head=1600;
q[head]=e[i].to;
}
else
{
q[tail++]=e[i].to;
if(tail==1601)tail=0;
}
}
}
inq[now]=0;
}
if(dis[T]==inf)return 0;
return 1;
}
void mcf()
{
int x=inf;
for(int i=from[T];i;i=from[e[i].from])
x=min(e[i].v,x);
for(int i=from[T];i;i=from[e[i].from])
{
ans+=x*e[i].c;
e[i].v-=x;e[i^1].v+=x;
}
}

zkw费用流
C++
bool spfa()
{
memset(mark,0,sizeof(mark));
for(int i=0;i<=T;i++)d[i]=-1;
int head=0,tail=1;
q[0]=T;mark[T]=1;d[T]=0;
while(head!=tail)
{
int now=q[head];head++;if(head==605)head=0;
for(int i=last[now];i;i=e[i].next)
if(e[i^1].v&&d[now]+e[i^1].c>d[e[i].to])
{
d[e[i].to]=d[now]+e[i^1].c;
if(!mark[e[i].to])
{
mark[e[i].to]=1;
q[tail++]=e[i].to;
if(tail==605)tail=0;
}
}
mark[now]=0;
}
return d[0]!=-1;
}
int dfs(int x,int f)
{
mark[x]=1;
if(x==T)return f;
int w,used=0;
for(int i=last[x];i;i=e[i].next)
if(d[e[i].to]==d[x]-e[i].c&&e[i].v&&!mark[e[i].to])
{
w=f-used;
w=dfs(e[i].to,min(w,e[i].v));
ans+=w*e[i].c;
e[i].v-=w;e[i^1].v+=w;
used+=w;if(used==f)return f;
}
return used;
}
void zkw()
{
while(spfa())
{
mark[T]=1;
while(mark[T])
{
memset(mark,0,sizeof(mark));
dfs(0,inf);
}
}
}
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52

bool spfa()
{
memset(mark,0,sizeof(mark));
for(int i=0;i<=T;i++)d[i]=-1;
int head=0,tail=1;
q[0]=T;mark[T]=1;d[T]=0;
while(head!=tail)
{
int now=q[head];head++;if(head==605)head=0;
for(int i=last[now];i;i=e[i].next)
if(e[i^1].v&&d[now]+e[i^1].c>d[e[i].to])
{
d[e[i].to]=d[now]+e[i^1].c;
if(!mark[e[i].to])
{
mark[e[i].to]=1;
q[tail++]=e[i].to;
if(tail==605)tail=0;
}
}
mark[now]=0;
}
return d[0]!=-1;
}
int dfs(int x,int f)
{
mark[x]=1;
if(x==T)return f;
int w,used=0;
for(int i=last[x];i;i=e[i].next)
if(d[e[i].to]==d[x]-e[i].c&&e[i].v&&!mark[e[i].to])
{
w=f-used;
w=dfs(e[i].to,min(w,e[i].v));
ans+=w*e[i].c;
e[i].v-=w;e[i^1].v+=w;
used+=w;if(used==f)return f;
}
return used;
}
void zkw()
{
while(spfa())
{
mark[T]=1;
while(mark[T])
{
memset(mark,0,sizeof(mark));
dfs(0,inf);
}
}
}