【bzoj2055】80人环游世界 上下界有源汇费用流

Description

Input

第一行两个正整数N，M。第二行有N个不大于M正整数，分别表示V1，V2……VN。接下来有N ¡ 1行。第i行有N ¡ i个整数，该行的第j个数表示从第i个国家到第i + j个国家的机票费（如果该值等于¡1则表示这两个国家间没有通航）。

Output

在第一行输出最少的总费用。

Sample Input

6 3
2 1 3 1 2 1
2 6 8 5 0
8 2 4 1
6 1 0
4 -1
4
Sample Output

27
HINT

1<= N < =100 1<= M <= 79

题解
懒
http://blog.csdn.net/clove_unique/article/details/54884437
http://blog.csdn.net/Clove_unique/article/details/54864211

代码

#include<cstdio>
#include<iostream>
#include<queue>
#define ll long long
#define inf 1000000000
#define mod 1000000007
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}   
int s,t,S,SS,TT;
queue<int>q;
int n,m,cnt=1,ans;
int Head[205],in[205],dis[205],from[205];
bool inq[205];
struct data{int to,from,next,v,c;}e[500005];
void ins(int u,int v,int w,int c)
{
    e[++cnt].to=v;e[cnt].from=u;
    e[cnt].next=Head[u];Head[u]=cnt;
    e[cnt].v=w;e[cnt].c=c;
}
void insert(int u,int v,int w,int c)
{ins(u,v,w,c);ins(v,u,0,-c);}
bool spfa()
{
    for (int i=1;i<=TT;i++)dis[i]=inf;
    q.push(SS);dis[SS]=0;inq[SS]=1;
    while (!q.empty())
    {
        int now=q.front();q.pop();
        for (int i=Head[now];i;i=e[i].next)
            if (e[i].v&&e[i].c+dis[now]<dis[e[i].to])
            {
                dis[e[i].to]=e[i].c+dis[now];
                from[e[i].to]=i;
                if (!inq[e[i].to])
                {
                    q.push(e[i].to);
                    inq[e[i].to]=1;
                }
            }
        inq[now]=0;
    }
    return dis[TT]!=inf;
}
void mcf()
{
    int x=inf;
    for (int i=from[TT];i;i=from[e[i].from])
        x=min(x,e[i].v);
    for (int i=from[TT];i;i=from[e[i].from])
        ans+=x*e[i].c,e[i].v-=x,e[i^1].v+=x;
}
int main()
{
    n=read();m=read();
    S=n+n+1;s=S+1;t=s+1;SS=t+1;TT=SS+1;
    in[s]-=m;in[S]+=m;
    for (int i=1;i<=n;i++)
    {
        int x=read();
        insert(S,i,inf,0);
        insert(i+n,t,inf,0);
        in[i]-=x;in[i+n]+=x;
    }
    for (int i=1;i<n;i++)
        for (int j=i+1;j<=n;j++)
        {
            int x=read();if (x==-1) continue;
            insert(n+i,j,inf,x);
        }
    for (int i=1;i<=t;i++)
        if (in[i]>0) insert(SS,i,in[i],0);
        else if (in[i]<0) insert(i,TT,-in[i],0);
    insert(t,s,inf,0);
    while (spfa()) mcf();
    printf("%d",ans);
    return 0;
}