[BZOJ1877]SDOI2009晨跑|费用流

       看到每个点只能经过一次就要想到网络流嘛。。把每个点拆成两个点x1x2,除1nx1都向x2连费用0流量1的边,1n连流量无穷。。然后若i可以到j,就i2j1连流量1费用为边权的边,注意这题是有向图。。建完图之后跑1-n的最小费用流即可,跑出来最大流就是天数,费用就是路程。。

#include
#include
#define N 405
#define M 20005
#define inf 0x3f3f3f3f
using namespace std;
struct edge{
	int s,e,f,q,next;
}ed[M*2+N*2];
int n,m,ne=1,s,i,e,q,nd,a[N],minf[N],d[N],inq[N],que[N],last[N],road[N];
void add(int s,int e,int f,int q)
{
	ed[++ne].e=e;ed[ne].f=f;
	ed[ne].q=q;ed[ne].next=a[s];a[s]=ne;
}
void ins(int s,int e,int f,int q)
{
	add(s,e,f,q);add(e,s,0,-q);
}
bool spfa(int s,int t)
{
	int head=1,tail=1,hh=1,tt=1,get,i,j,to;
	for (i=1;i<=nd;i++) inq[i]=0,d[i]=inf,last[i]=road[i]=-1;
	que[1]=s;inq[s]=1;d[s]=0;minf[s]=inf;
	while (hh<=tt)
	{
		get=que[head++];hh++;
		if (head>400) head=1;
		for(j=a[get];j;j=ed[j].next)
			if (ed[j].f&&d[get]+ed[j].q400) tail=1;
					que[tail]=to;inq[to]=1;
				}
			}
		inq[get]=0;
	}
	return last[t]!=-1;
}
void fyl(int s,int t)
{
	int flow=0,cost=0,i,j;
	while (spfa(s,t))
	{
		flow+=minf[t];cost+=minf[t]*d[t];
		for (i=t;last[i]!=-1;i=last[i])
		{
			j=road[i];
			ed[j].f-=minf[t];ed[j^1].f+=minf[t];
		}
	}
	printf("%d %d\n",flow,cost);
}
int main()
{
	freopen("1877.in","r",stdin);
	scanf("%d%d",&n,&m);
	nd=2*n;
	for (i=1;i<=n;i++) a[i]=0;
	for(i=1;i<=m;i++)
	{
		scanf("%d%d%d",&s,&e,&q);
		ins(s+n,e,1,q);
	}
	for(i=1;i<=n;i++) if (i==1||i==n) ins(i,i+n,inf,0);else ins(i,i+n,1,0);
	fyl(1,nd);
}


你可能感兴趣的:(BZOJ)