网络流&费用流模版

来源：互联网发布：淘宝购物省钱技巧编辑：程序博客网时间：2024/05/16 17:57

我写的太烂了，大家看这个人写的吧QAQ

有时间就来看看QAQ

网络流dinic

C++
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
bool bfs()
{
int head=0,tail=1;
for(int i=0;i<=T;i++)h[i]=-1;
q[0]=0;h[0]=0;
while(head!=tail)
{
int now=q[head];head++;
for(int i=last[now];i;i=e[i].next)
if(e[i].v&&h[e[i].to]==-1)
{
h[e[i].to]=h[now]+1;
q[tail++]=e[i].to;
}
}
return h[T]!=-1;
}
int dfs(int x,int f)
{
if(x==T)return f;
int w,used=0;
for(int i=cur[x];i;i=e[i].next)
if(h[e[i].to]==h[x]+1)
{
w=f-used;
w=dfs(e[i].to,min(w,e[i].v));
e[i].v-=w;e[i^1].v+=w;
if(e[i].v)cur[x]=i;
used+=w;if(used==f)return f;
}
if(!used)h[x]=-1;
return used;
}
void dinic()
{
while(bfs())
{
for(int i=0;i<=T;i++)
cur[i]=last[i];
ans+=dfs(0,inf);
}
}

最小费用最大流spfa

C++
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
boolspfa()
{
for(inti=0;i<=T;i++)dis[i]=inf;
inthead=0,tail=1;
dis[0]=0;q[0]=0;inq[0]=1;
while(head!=tail)
{
intnow=q[head++];if(head==1601)head=0;
for(inti=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)return0;
return1;
}
voidmcf()
{
intx=inf;
for(inti=from[T];i;i=from[e[i].from])
x=min(e[i].v,x);
for(inti=from[T];i;i=from[e[i].from])
{
ans+=x*e[i].c;
e[i].v-=x;e[i^1].v+=x;
}
}

zkw费用流

C++
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
boolspfa()
{
    memset(mark,0,sizeof(mark));
    for(inti=0;i<=T;i++)d[i]=-1;
    inthead=0,tail=1;
    q[0]=T;mark[T]=1;d[T]=0;
    while(head!=tail)
    {
intnow=q[head];head++;if(head==605)head=0;
for(inti=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;
    }
    returnd[0]!=-1;
}
intdfs(intx,intf)
{
    mark[x]=1;
    if(x==T)returnf;
    intw,used=0;
    for(inti=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)returnf;
}
    returnused;
}
voidzkw()
{
    while(spfa())
    {
mark[T]=1;
while(mark[T])
{
memset(mark,0,sizeof(mark));
dfs(0,inf);
}
    }
}

阅读全文

0 0