讨论bellman

咳咳，今天要聊一聊bellman
在bellman中模板一般都是这个样子

for(int k=1;k<=n-1;k++)
{
    flag=0;
    for(int j=1;j<=m;j++)
    {
        if(dis[v[j]]>dis[u[j]]+w[j])
        {
            dis[v[j]]=dis[u[j]]+w[j];
            flag=1;
        }
    }
    if(flag==0)
    {
        break;
    }
}

但这组测试数据中所有边值都是负数：

6 5
2 4 -10
4 6 -10
2 5 -10
5 6 -10
3 6 -10
我们发现这个图根本不连通，但是普通的bellman算法运行结果是这两点可以连通 所以在有负边时应判断1号点到u号点连不连通，可以避免此情况；

http://bbs.codeaha.com/problem-22100.html
这坑题 不加判定过不了啊555
当然这题用bellman会超时（80分，加了判定30变80 ），spfa是正解
网络流标准题，不是求最小吗？把边改成负数呗

#include <cstdio>
#include <algorithm>
using namespace std;
long long  lev[160001],father[160001],cnt=0,k=0,minn,dis[160001],num[160001];
long long n,m,p,a,b,person,st=0,ed,mm,now,ls,flag=0,sum=0,ans=0;//初始化
struct node
{
    int start;
    int end;
    int flow;
    int money;
}edge[160001];
void init2()
{
    int ll;
    for(ll=0;ll<=ed;ll++)
    {
        father[ll]=-1;
        dis[ll]=999999999;
    }
    dis[0]=0;
    return;
}
int bellman()
{
    init2();
    int k,i;
    for(k=1;k<=ed;k++)
    {
        flag=0;
        for(i=1;i<=mm;i++)
        {
            if(edge[i].flow>0)
            {
                if(dis[edge[i].start]!=999999999)//这里很重要！！
                {
                    if(dis[edge[i].end]>dis[edge[i].start]+edge[i].money)
                    {
                        dis[edge[i].end]=dis[edge[i].start]+edge[i].money;
                        father[edge[i].end]=i;
                        flag=1;
                    }
                }
            }
        }
        if(flag==0)
        {
            break;
        }
    }
    if(dis[ed]>=999999999)
    {
        return 0;
    }
    return 1;
}
void dfs()
{
    while(bellman()==1)
    {
        minn=2147483647;
        ls=ed;
        while(father[ls]!=-1)
        {
            now=father[ls];
            edge[now].flow-=1;
            edge[num[now]].flow+=1;
            ans+=edge[now].money;
            ls=edge[now].start;
        }
    }
}
void addedge(int s,int e,int f,int m,int u,int v)
{
    edge[u].start=s;
    edge[u].end=e;
    edge[u].flow=f;
    edge[u].money=m;
    edge[u+v].start=e;
    edge[u+v].end=s;
    edge[u+v].flow=0;
    edge[u+v].money=-m;
    num[u]=u+v;
    num[u+v]=u;
    return;
}
int main()
{
    scanf("%lld%lld%lld",&n,&p,&m);
    for(int i=1;i<=n;i++)
    {
        scanf("%lld",&lev[i]);
    }
    ed=n+p+1;
    for(int i=1;i<=n;i++)
    {
        addedge(0,i,1,0,i,n);
    }
    mm=n*2;
    k=1;
    for(int i=mm+1;i<=mm+p;i++)
    {
        addedge(n+k,ed,1,0,i,p);
        k++;
    }
    mm=mm+p*2;
    for(int i=mm+1;i<=mm+m;i++)
    {
        scanf("%lld%lld%lld",&a,&b,&person);
        b=b+n;
        addedge(a,b,1,-person*lev[a],i,m);
    }
    mm=mm+2*m;
    dfs();
    printf("%lld",-ans);
    return 0;
}