[인터넷 흐름 24문제] 운송 문제.

(인터넷 흐름 24문제는 대부분 spj가 필요하기 때문에 spj가 있는 oj가 필요합니다. 이 시리즈 코드는 모두 www.oj.swust.edu.cn 테스트를 통과했습니다) 매우 누드적인 비용 흐름입니다. 단지 최대 흐름과 최소 흐름을 한 번 뛰어야 하기 때문에 여기서 군말하지 말고 코드를 보십시오.
#include
#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;
#define INF 100000000
struct bian
{
    int l,r,f,v;
}a[1000000];
int fir[1000000];
int nex[1000000];
int tot=1;
void add_edge(int l,int r,int f,int v)
{
    a[++tot].l=l;
    a[tot].r=r;
    a[tot].f=f;
    a[tot].v=v;
    nex[tot]=fir[l];
    fir[l]=tot;
}
int n,m;
int v1[1001];
int v2[1001];
int val[1001][1001];
int s=0,t=9999;
void jianbian()
{
    memset(fir,0,sizeof(fir));
    memset(nex,0,sizeof(nex));
    tot=1;
    for(int i=1;i<=n;i++)
    {
        add_edge(s,i,v1[i],0);
        add_edge(i,s,0,0);
    }
    for(int i=1;i<=m;i++)
    {
        add_edge(n+i,t,v2[i],0);
        add_edge(t,n+i,0,0);
    }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            add_edge(i,n+j,INF,val[i][j]);
            add_edge(n+j,i,0,-val[i][j]);
        }
}
int dis[10000];
bool pd[10000];
int fro[10000];
bool spfamin()
{
    static int dui[1000000];
    memset(dis,0x1f,sizeof(dis));
    int top=1,my_final=2;
    dui[top]=s;
    dis[s]=0;
    pd[s]=true;
    while(topint u=dui[top++];
        for(int o=fir[u];o;o=nex[o])
        {
            if(dis[a[o].r]>dis[u]+a[o].v && a[o].f)
            {
                dis[a[o].r]=dis[u]+a[o].v;
                fro[a[o].r]=o;
                if(!pd[a[o].r]) pd[a[o].r]=true,dui[my_final++]=a[o].r;
            }
        }
        pd[u]=false;
    }
    if(dis[t]==0x1f1f1f1f) return false;
    return true;
}
bool spfamax()
{
    static int dui[1000000];
    memset(dis,0x8f,sizeof(dis));
    int top=1,my_final=2;
    dui[top]=s;
    dis[s]=0;
    pd[s]=true;
    while(topint u=dui[top++];
        for(int o=fir[u];o;o=nex[o])
        {
            if(dis[a[o].r]if(!pd[a[o].r]) pd[a[o].r]=true,dui[my_final++]=a[o].r;
            }
        }
        pd[u]=false;
    }
    if(dis[t]==0x8f8f8f8f) return false;
    return true;
}
int cost;
void add_flow()
{
    int mid=t;
    int temp=2147483647;
    while(mid!=s)
    {
        temp=min(temp,a[fro[mid]].f);
        mid=a[fro[mid]^1].r;
    }
    cost+=temp*dis[t];
    mid=t;
    while(mid!=s)
    {
        a[fro[mid]].f-=temp;
        a[fro[mid]^1].f+=temp;
        mid=a[fro[mid]^1].r;
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&v1[i]);
    for(int i=1;i<=m;i++) scanf("%d",&v2[i]);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            scanf("%d",&val[i][j]);
    jianbian();
    while(spfamin()) add_flow();
    cout<0;
    while(spfamax()) add_flow();
    cout<return 0;
}

좋은 웹페이지 즐겨찾기