就是纯粹的spfa啦。没森么好说的啦。
就是建边有点难而已。
其实我建边想了一天……
祝你Ak~QaQ
#include<cmath>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
using namespace std;
struct node{
int x,y,c,d,other,next;
}a[210000];int len,first[210000];
void ins(int x,int y,int c,int d){
int kx,ky;
len++;kx=len;
a[len].x=x;a[len].y=y;a[len].c=c;a[len].d=d;
a[len].next=first[x];first[x]=len;
len++;ky=len;
a[len].x=y;a[len].y=x;a[len].c=0;a[len].d=-d;
a[len].next=first[y];first[y]=len;
a[kx].other=ky;a[ky].other=kx;
}
int list[110000],d[110000],c[110000],pre[110000],head,tail,st,ed,n,m,xx,ans;
bool v[110000];
bool spfa(){
memset(c,0,sizeof(c));c[st]=999999999;
memset(v,false,sizeof(v));v[st]=true;
memset(d,0x3F,sizeof(d));d[st]=0;
list[1]=st;head=1;tail=2;pre[st]=0;
while(head!=tail){
int x=list[head];
for(int k=first[x];k>0;k=a[k].next){
int y=a[k].y;
if(a[k].c>0&&d[y]>d[x]+a[k].d){
d[y]=d[x]+a[k].d;
c[y]=min(c[x],a[k].c);
pre[y]=k;
if(v[y]==false){
v[y]=true;
list[tail++]=y;
if(tail==ed+1)tail=1;
}
}
}
head++;if(head==ed+1)head=1;
v[x]=false;
}
if(d[ed]>999999999)return false;
ans+=d[ed]*c[ed];
int x=ed;
while(x!=st){
int k=pre[x];
a[k].c-=c[ed];a[a[k].other].c+=c[ed];
x=a[k].x;
}
return true;
}
int map[1100][1100];
int main(){
scanf("%d%d",&m,&n);
st=n*m+n+1;ed=n*m+n+2;
len=0;memset(first,0,sizeof(first));
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++){
scanf("%d",&map[i][j]);
}
}
for(int i=1;i<=n;i++)ins(st,i,1,0);
for(int i=n+1;i<=n+n*m;i++)ins(i,ed,1,0);
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++){
for(int k=1;k<=n;k++){
ins(i,j*n+k,1,map[i][j]*(n-k+1));
}
}
}
ans=0;
while(spfa()==true);
printf("%.2lf\n",double(ans)/double(n));
return 0;
}