poj 3422 Kaka's Matrix Travels

给出一个n*n的矩阵，每个点有权值。一个人要从左上角点每次向下或向右走走到右下角点，每到一个点拿走该点权值。现在他要走k遍，每个点的权值只能取一次，问能收获的最大权值。

1<=n<=50.0<=k<=10.

最小费用最大流。

将每个点拆点。

源点向左上入点连一条容量为k费用为0的边。

每个点的入点向出点连一条容量为1费用为此点权值相反数的边。

每个点的入点向出点连一条容量为正无穷费用为0的边。

每个点的出点向右边入点和下边入点连一条容量为正无穷费用为0的边。

右下出点向汇点连一条容量为正无穷费用为0的边。

答案为最小费用最大流的相反数。

 1 #include<iostream>
 2 #include<cstring>
 3 #include<cstdio>
 4 #include<algorithm>
 5 #include<queue>
 6 using namespace std;
 7 const int dian=5005;
 8 const int bian=40005;
 9 const int INF=0x3f3f3f3f;
10 int h[dian],nxt[bian],ver[bian],val[bian],cos[bian],minn[dian],with[dian];
11 int v[dian],d[dian];
12 int n,k,a,tot;
13 int S,T;
14 int add(int a,int b,int c,int d){
15     tot++;ver[tot]=b;val[tot]=c;cos[tot]=d,nxt[tot]=h[a];h[a]=tot;
16     tot++;ver[tot]=a;val[tot]=0;cos[tot]=-d;nxt[tot]=h[b];h[b]=tot;
17 }
18 bool tell(){
19     memset(v,0,sizeof(v));
20     memset(d,0x3f,sizeof(d));
21     memset(with,0,sizeof(with));
22     memset(minn,0x3f,sizeof(minn));
23     queue<int>q;
24     q.push(S);
25     v[S]=1;
26     d[S]=0;
27     while(!q.empty()){
28         int x=q.front();
29         q.pop();
30         v[x]=0;
31         for(int i=h[x];i;i=nxt[i]){
32             int y=ver[i];
33             if(d[y]>d[x]+cos[i]&&val[i]){
34                 d[y]=d[x]+cos[i];
35                 minn[y]=min(minn[x],val[i]);
36                 with[y]=i;
37                 if(!v[y]){
38                     v[y]=1;
39                     q.push(y);
40                 }
41             }
42         }
43     }
44     if(d[T]!=0x3f3f3f3f)
45         return true;
46     return false;
47 }
48 int zeng(){
49     for(int i=T;i!=S;i=ver[with[i]^1]){
50         val[with[i]]-=minn[T];
51         val[with[i]^1]+=minn[T];
52     }
53     return d[T]*minn[T];
54 }
55 int dinic_cost(){
56     int r=0;
57     while(tell())
58         r+=zeng();
59     return r;
60 }
61 int main(){
62     memset(h,0,sizeof(h));
63     memset(nxt,0,sizeof(nxt));
64     tot=1;
65     scanf("%d%d",&n,&k);
66     S=n*n*2+1,T=n*n*2+2;
67     for(int i=1;i<=n;i++)
68         for(int j=1;j<=n;j++){
69             scanf("%d",&a);
70             add((i-1)*n+j,(i-1)*n+j+n*n,1,-a);
71             add((i-1)*n+j,(i-1)*n+j+n*n,INF,0);
72             if(j!=n)
73                 add((i-1)*n+j+n*n,(i-1)*n+j+1,INF,0);
74             if(i!=n)
75                 add((i-1)*n+j+n*n,i*n+j,INF,0);
76         }
77     add(S,1,k,0);
78     add(n*n*2,T,INF,0);
79     printf("%d",-dinic_cost());
80     return 0;
81 }

 

转载于:https://www.cnblogs.com/dugudashen/p/6242110.html