bzoj 1070: [SCOI2007]修车

Posted 2020-06-17

 1 #include<cstdio>
 2 #include<iostream>
 3 #include<cstring>
 4 #define M 10000
 5 #define inf 2139062143
 6 using namespace std;
 7 int cnt=1,n,m,T,d[M],q[2*M],f[M],head[M],next[10*M],u[10*M],v[10*M],w[10*M],fro[10*M],fr[M];
 8 int mp[100][100];
 9 int ans;
10 void jia1(int a1,int a2,int a3,int a4)
11 {
12     cnt++;
13     next[cnt]=head[a1];
14     head[a1]=cnt;
15     fro[cnt]=a1;
16     u[cnt]=a2;
17     v[cnt]=a3;
18     w[cnt]=a4;
19 }
20 void jia(int a1,int a2,int a3,int a4)
21 {
22     jia1(a1,a2,a3,a4);
23     jia1(a2,a1,0,-a4);
24     return;
25 }
26 bool spfa()
27 {
28     memset(d,127,sizeof(int)*(T+1));
29     d[0]=0;
30     f[0]=1;
31     q[1]=0;
32     int h=0,t=1;
33     for(;h<t;)
34       {
35         h++;
36         int p=q[h];
37         f[p]=0;
38           for(int i=head[p];i;i=next[i])
39           if(v[i]&&d[u[i]]>d[p]+w[i])
40             {
41                 d[u[i]]=d[p]+w[i];
42                 fr[u[i]]=i;
43                 if(!f[u[i]])
44                   {
45                     f[u[i]]=1;
46                     t++;
47                     q[t]=u[i];
48                     }
49             }
50       }
51     if(d[T]!=inf)
52       return 1;
53     return 0;
54 }
55 void mcf()
56 {
57     int mx=inf;
58     for(int i=fr[T];i;i=fr[fro[i]])
59       mx=min(mx,v[i]);
60     for(int i=fr[T];i;i=fr[fro[i]])
61       {
62         v[i]-=mx;
63         v[i^1]+=mx;
64         ans+=mx*w[i];
65       }
66     return;
67 }
68 int main()
69 {
70     scanf("%d%d",&n,&m);
71     T=m+m*n+1;
72     for(int i=1;i<=m;i++)
73       jia(0,i,1,0);
74     for(int i=1;i<=m;i++)
75       for(int j=1;j<=n;j++)
76         scanf("%d",&mp[i][j]);
77     for(int i=1;i<=m;i++)
78       for(int j=1;j<=n;j++)
79         for(int k=1;k<=m;k++)
80           jia(i,(j-1)*m+k+m,1,mp[i][j]*k);
81     for(int i=1;i<=n;i++)
82       for(int j=1;j<=m;j++)
83         jia((i-1)*m+j+m,T,1,0);
84     for(;spfa();)
85       mcf();
86     printf("%.2lf\n",ans/(m*1.0));
87     return 0;
88 }

拆点网络流，将每个技术人员拆成顾客数个，分别代表倒数第一个修，倒数第二个修，。。。。；顾客向他们连边权值为他及以后的他消耗的等待时间。