JZOJ 4802. 探险计划

---

题目：

传送门

---

题意：

给出山上每个点的危险系数，求在不经过同样的点和路、不经过同样的路但可以经过同样的点的两种情况的最小危险系数

---

分析：

爽快费用流
对于情况1，因为每个点只能走一次所以每个点对自己连一条容量为$1$的边。而因为每条边也只能经过$1$次，所以它们的容量也是$1$
对于情况2，因为每个点都可以经过多次，所以我们就把容量设为$+无穷$，而其他的按照情况1就好了

---

代码：

#pragma GCC optimize(2)
#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#define LL long long 
#define LZX LU
#define f(x,y,z) 2*((y-1)*n+x)-z 
inline LL read() {
    LL d=0,f=1;char s=getchar();
    while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}
    while(s>='0'&&s<='9'){d=d*10+s-'0';s=getchar();}
    return d*f;
}
using namespace std;
struct node{
	int to,next,w,c;
}a[320000];
int ans,n,m,s,e,tot=1,dan[200][400],t;
int f[80000],mf[80000],ls[80000],pre[80000];
bool v[80000];
queue<int> q;
void add(int x,int y,int w,int c)
{
	a[++tot]=(node){y,ls[x],w,c};
	ls[x]=tot;
	a[++tot]=(node){x,ls[y],0,-c};
	ls[y]=tot;
}
bool spfa()
{
	memset(f,0x3f,sizeof(f));
	mf[s]=2147483647/2;f[s]=0;
	v[s]=1;q.push(s);
	pre[e]=0;
	while(!q.empty())
	{
		int x=q.front();q.pop();v[x]=0;
		for(int i=ls[x];i;i=a[i].next)
		{
			if(!a[i].w) continue;
			int y=a[i].to;
			if(f[x]+a[i].c<f[y])
			{
				f[y]=f[x]+a[i].c;
				mf[y]=min(a[i].w,mf[x]);
				pre[y]=i;
				if(!v[y])
				{
					q.push(y);
					v[y]=1;
				}
			}
		}
	}
	return pre[e];
}
void over_path()
{
	int now=e,flow=mf[e];
	ans+=mf[e]*f[e];
	while(now!=s)
	{
		a[pre[now]].w-=flow;
		a[pre[now]^1].w+=flow;
		now=a[pre[now]^1].to;
	}
}
void Net_flow()
{
	while(spfa())
	  over_path();
}
int main()
{
	n=read();m=read();
	e=0;s=f(n,n+m-1,1)+1;t=s+1;
	for(int i=1;i<=n;i++)
	  for(int j=1;j<=m+i-1;j++)
	  {
	  	dan[i][j]=read();
		if(i==1) add(f(i,j,1),t,2147483647/2,0);
	  	if(i==n) add(s,f(i,j,0),2147483647/2,0);
	  	add(f(i,j,0),f(i,j,1),1,dan[i][j]);
	  	if(j<m+i-1&&i>1) add(f(i,j,1),f(i-1,j,0),1,0);
	  	if(j>1&&i>1) add(f(i,j,1),f(i-1,j-1,0),1,0);
	  }
	add(t,e,m,0);
	Net_flow();
	printf("%d\n",ans);
	ans=0;tot=1;memset(ls,0,sizeof(ls));
	for(int i=1;i<=n;i++)
	  for(int j=1;j<=m+i-1;j++)
	  {
		if(i==1) add(f(i,j,1),t,m,0);
	  	if(i==n) add(s,f(i,j,0),m,0);
	  	add(f(i,j,0),f(i,j,1),m,dan[i][j]);
	  	if(j<m+i-1&&i>1) add(f(i,j,1),f(i-1,j,0),1,0);
	  	if(j>1&&i>1) add(f(i,j,1),f(i-1,j-1,0),1,0);
	  }
	add(t,e,m,0);
	Net_flow();
	printf("%d",ans);
}