TYVJ 1741 [SDOI2010]星际竞速 最小费用流

这个题联系的时候没有做出来。最后还是lyd神犇教的。。。我还是太弱了。。

每个点拆成一个入点(右边一排)和一个出点(左边一排),源点S向每个出点连容量1费用0的边,每个入点向汇点连容量1费用0的边,如果有高速航道(x,y),注意x<y,那么从x的出点到y的入点连容量1费用为读入权值的边。这样就能够处理所有高速航道了……

空间跳跃处理:从S向所有入点连容量1费用为定位费用的边就行了。

流过一个点对应的出点不代表经过了这个点,只有流过了这个点对应的的入点才算真正流过了这个点!

View Code
 1 #include <iostream>
 2 #include <algorithm>
 3 #include <cstring>
 4 #include <cstdio>
 5 #include <cstdlib>
 6 
 7 #define N 10000
 8 #define M 1000000
 9 
10 using namespace std;
11 
12 int head[N],next[M],to[M],pr[M],len[M];
13 int q[M*5],dis[N],pre[N];
14 bool vis[N];
15 int n,m,cnt,S,T;
16 
17 inline void add(int u,int v,int r,int w)
18 {
19     to[cnt]=v; len[cnt]=r; pr[cnt]=w; next[cnt]=head[u]; head[u]=cnt++;
20     to[cnt]=u; len[cnt]=0; pr[cnt]=-w; next[cnt]=head[v]; head[v]=cnt++;
21 }
22 
23 inline void read()
24 {
25     memset(head,-1,sizeof head); cnt=0;
26     scanf("%d%d",&n,&m);
27     S=0; T=2*n+1;
28     for(int i=1,a;i<=n;i++)
29     {
30         scanf("%d",&a);
31         add(S,i+n,1,a);
32     }
33     for(int i=1,a,b,c;i<=m;i++)
34     {
35         scanf("%d%d%d",&a,&b,&c);
36         if(a>b) swap(a,b);
37         add(a,b+n,1,c);
38     }
39     for(int i=1;i<=n;i++) add(S,i,1,0),add(i+n,T,1,0);
40 }
41 
42 inline bool spfa()
43 {
44     memset(dis,0x3f,sizeof dis);
45     memset(pre,-1,sizeof pre);
46     int h=1,t=2,sta;
47     q[1]=S; vis[S]=true; dis[S]=0;
48     while(h<t)
49     {
50         sta=q[h++]; vis[sta]=false;
51         for(int i=head[sta];~i;i=next[i])
52             if(len[i]&&dis[to[i]]>dis[sta]+pr[i])
53             {
54                 dis[to[i]]=dis[sta]+pr[i];
55                 pre[to[i]]=i;
56                 if(!vis[to[i]]) vis[to[i]]=true,q[t++]=to[i];
57             }
58     }
59     return pre[T]!=-1;
60 }
61 
62 inline void updata()
63 {
64     for(int i=pre[T];~i;i=pre[to[i^1]])
65     {
66         len[i]-=1; len[i^1]+=1;
67     }
68 }
69 
70 inline void go()
71 {
72     int ans=0;
73     while(spfa()) ans+=dis[T],updata();
74     printf("%d\n",ans);
75 }
76 
77 int main()
78 {
79     read();
80     go();
81     return 0;
82 }
没有人能阻止我前进的步伐,除了我自己!
原文地址:https://www.cnblogs.com/proverbs/p/2848043.html