hz2016工作室

欢迎关注hz2016.cn博客搬迁了,有兴趣的可以去网站私聊我。

[SCOI2007]修车

就是纯粹的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;
}
阅读更多
版权声明:原创文章,请转载表明链接 https://blog.csdn.net/huangzhengdoc/article/details/52353739
文章标签: spfa
个人分类: BZOJ spfa
所属专栏: C++算法大全
相关热词: scoi2007
上一篇[SCOI2005]最大子矩阵
下一篇[SCOI2007]组队
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

关闭
关闭