[HDU 2448] Mining Station on the Sea

题目描述：

有N个海港，N个船只，M个矿藏点，有K条海路连接矿藏点，每条路都有一个权值，并且海路双向，有P条路连接矿藏点和海港，只能从矿藏点驶向海港，每条路有一个权值，N条船初始在N个矿藏点上，现在要求你求出N个船进入N个海港的最小代价，保证有解。

题目分析：

不就是个费用流吗.
首先，我们从源点向每个有船只的矿藏点连一条容量为1，权值为0的边
然后根据矿藏点之间的路径信息，连接容量为INF，权值为C的双向边，
然后根据矿藏点和海港之间的路径信息，连接容量为1（INF也行）权值为C的单向边
然后每个海港向汇店连容量为1，权值为0的边。
跑最小费用流就Vans辣

题目链接：

HDU 2448

代码：

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <queue>
#include <string>
const int maxm=300000; 
const int N=400;
const int inf=0x7fffffff;
int head[N],net[maxm],to[maxm],cap[maxm],cost[maxm];
int tx[maxm];
int cnt=1;
std::queue <int> dl;
int flow[N],id[N],pre[N],dis[N],vis[N];
int n,m,k,p;
inline char fgc() 
{
    static char buf[100000],*p1 = buf, *p2 = buf;
    return p1==p2&&(p2 = (p1 = buf)+fread(buf, 1, 100000, stdin),p1==p2)?EOF:*p1++;
}
inline int read()
{
    register int res=0,neg=1;
    register char c=fgc();
    while(!isdigit(c))
    {
        if(c=='-') neg=-1;
        c=fgc();
    }
    while(isdigit(c))
    {
        res=(res << 1)+(res << 3)+c-'0';
        c=fgc();
    }
    return res*neg;
}
inline void addedge(int u,int v,int c1,int c2)
{
    cnt++;
    to[cnt]=v,cap[cnt]=c1,cost[cnt]=c2,net[cnt]=head[u],head[u]=cnt;
    cnt++;
    to[cnt]=u,cap[cnt]=0,cost[cnt]=-c2,net[cnt]=head[v],head[v]=cnt;
} 
inline bool SPFA(int s,int t)
{
    memset(pre,-1,sizeof(pre));
    memset(dis,0x3f,sizeof(dis));
    dis[s]=0,flow[s]=inf,pre[s]=0,vis[s]=1;
    dl.push(s);
    while(!dl.empty())
    {
        int now=dl.front();
        dl.pop();
        vis[now]=0;
        for(int i=head[now];i;i=net[i])
        if(cap[i]&&dis[to[i]]>dis[now]+cost[i])
        {
            dis[to[i]]=dis[now]+cost[i];
            flow[to[i]]=std::min(flow[now],cap[i]);
            pre[to[i]]=now;
            id[to[i]]=i;
            if(!vis[to[i]]) vis[to[i]]=1,dl.push(to[i]);
        }
    }
    return pre[t]!=-1;
}
inline void change_cap(int s,int t,int x)
{
    int now=t;
    int a,b;
    while(now!=s)
    {
        cap[id[now]]-=x,cap[id[now]^1]+=x;
        now=pre[now];
    }
}
inline int MCMF(int s,int t)
{
    int maxflow=0,mincost=0;
    while(SPFA(s,t))
    {
        maxflow+=flow[t],mincost+=flow[t]*dis[t];
        change_cap(s,t,flow[t]);
    }
    return mincost;
}
inline void work()
{
	memset(head,0,sizeof(head));
	int s,t;
    cnt=1;
	s=n+m+1,t=n+m+2;
	for(int i=1;i<=n;i++) 
	{
		int u;
		scanf("%d",&u);
		addedge(s,u,1,0);
	}
	for(int i=1;i<=k;i++)
	{
		int u,v,c;
		scanf("%d%d%d",&u,&v,&c);
		addedge(u,v,inf,c);
		addedge(v,u,inf,c);
	}
	for(int i=1;i<=p;i++)
	{
		int u,v,c;
		scanf("%d%d%d",&u,&v,&c);
		addedge(v,u+m,1,c);
	}
	for(int i=1;i<=n;i++) addedge(m+i,t,1,0);
	printf("%d\n",MCMF(s,t));
}
signed main()
{
    while(scanf("%d%d%d%d",&n,&m,&k,&p)!=EOF) work();
    return 0; 
}