记录编号 |
281506 |
评测结果 |
AAAAAAAAAA |
题目名称 |
[SCOI 2007] 修车 |
最终得分 |
100 |
用户昵称 |
TenderRun |
是否通过 |
通过 |
代码语言 |
C++ |
运行时间 |
0.173 s |
提交时间 |
2016-07-11 20:07:35 |
内存使用 |
9.29 MiB |
显示代码纯文本
#include <iostream>
#include <cstring>
#include <cstdio>
#include <queue>
using namespace std;
const int INF=1000000000;
const int maxn=1010,maxm=400010;
int cnt=1,fir[maxn],nxt[maxm],to[maxm];
int cap[maxm],val[maxm],dis[maxn],path[maxn];
void add(int a,int b,int c,int v){
nxt[++cnt]=fir[a];to[cnt]=b;
cap[cnt]=c;val[cnt]=v;fir[a]=cnt;
}
void addedge(int a,int b,int c,int v){
add(a,b,c,v);
add(b,a,0,-v);
}
int S,T;
int vis[maxn];
int Spfa(){
deque<int>q;
memset(dis,127,sizeof(dis));
memset(vis,0,sizeof(vis));
q.push_front(S);
dis[S]=0;vis[S]=1;
while(!q.empty()){
int x=q.front();q.pop_front();vis[x]=0;
for(int i=fir[x];i;i=nxt[i])
if(cap[i]&&dis[x]+val[i]<dis[to[i]]){
dis[to[i]]=val[i]+dis[x];
path[to[i]]=i;
if(vis[to[i]])continue;
if(dis[to[i]]<dis[x])
q.push_front(to[i]);
else
q.push_back(to[i]);
vis[to[i]]=1;
}
}
return dis[T]==dis[T+1]?0:dis[T];
}
int Aug(){
int p=T,f=INF;
while(p!=S){
f=min(f,cap[path[p]]);
p=to[path[p]^1];
}
p=T;
while(p!=S){
cap[path[p]]-=f;
cap[path[p]^1]+=f;
p=to[path[p]^1];
}
return f;
}
int MCMF(){
int ret=0,d;
while(d=Spfa())
ret+=Aug()*d;
return ret;
}
int n,m;
int a[maxn][maxn];
int main(){
freopen("scoi2007_repair.in","r",stdin);
freopen("scoi2007_repair.out","w",stdout);
scanf("%d%d",&m,&n);
S=0;T=n*m+n+1;
for(int i=1;i<=n;i++){
addedge(n*m+i,T,1,0);
for(int j=1;j<=m;j++)
addedge(S,(i-1)*m+j,1,0);
}
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
scanf("%d",&a[i][j]);
for(int j=1;j<=m;j++)
for(int i=1;i<=n;i++)
for(int k=1;k<=n;k++)
addedge((i-1)*m+j,n*m+k,1,i*a[k][j]);
printf("%.2f\n",1.0*MCMF()/n);
return 0;
}