luogu P4013 数字梯形问题

关于构图：

这是一道最大费用最大流的好题，包括拆点和不拆点。

我们需要理解题目中的“相交”是指走的边不重复（本题最难理解），那很容易，即流量为1.

对于第一小问：

1拆点，超级源点和汇点,st=0,ed=tot*2+1（其中tot指梯形节点的个数）。

2.从i向i'建一条连接一条容量为1，费用为点i值的边。 

3.从源点向梯形顶层每个节点建一条容量为1，费用为0的边。 

4.从梯形底层每个节点向汇点建一条容量为1，费用为0的边。 

5.对于每个点i和左下、右下方向两个点j，分别建一条从i'到j（这样它就可以流回去，进行下一次选择）容量为1，费用为0的边，表示有且仅可走一次。 

最后跑最大费用最大流即可。

对于第二小问：

1超级源点和汇点,st=0,ed=tot+1（其中tot指梯形节点的个数）。

2.从源点向梯形顶层每个节点建一条容量为1，费用为0的边。 

3.从梯形底层每个节点i向汇点建一条容量为1，费用为点i值的边。 

4.对于每个点i和左下、右下方向两个点j，分别建一条从i 到j 容量为1 ，费用为点i值的边，表示有且仅可走一次。 

最后跑最大费用最大流即可。

对于第三小问：

1超级源点和汇点,st=0,ed=tot+1（其中tot指梯形节点的个数）。

2.从源点向梯形顶层每个节点建一条容量为1，费用为0的边。 

3.从梯形底层每个节点i向汇点建一条容量为INF，费用为点i值的边，表示可以走多次。 

4.对于每个点i和左下、右下方向两个点j，分别建一条从i到j容量为INF，费用为点i值的边，表示可以走多次。 

最后跑最大费用最大流即可。

最后，推荐两篇好的博客

（1.点击打开链接

2.点击打开链接）

关于代码：

#include<cstdio>
#include<queue>
#include<cstring>
#include<queue>
#include<algorithm>
#define INF 2147483647
using namespace std;
queue<int> f;
	int n,m,tot,len,st,ed;
	struct node{int x,y,c,d,next;} a[10000];
	int b[30][30],id[30][30],last[10000],dis[10000],pre[10000],pos[10000],p[10000];
	bool bz[10000];
void ins(int x,int y,int c,int d)
{
	a[++len].x=x;a[len].y=y;a[len].c=c;a[len].d=d;a[len].next=last[x];last[x]=len;
}
bool spfa()
{
	memset(bz,true,sizeof(bz));
	bz[st]=false;
	memset(dis,63,sizeof(dis));
	dis[st]=0;
	p[st]=INF;
	f.push(st);
	while(!f.empty())
	{
		int x=f.front();
		bz[x]=true;
		for(int i=last[x];i>=0;i=a[i].next)
		{
			int y=a[i].y;
			if(a[i].c>0&&dis[y]>dis[x]+a[i].d)
			{
				dis[y]=dis[x]+a[i].d;
				pos[y]=x;
				pre[y]=i;
				p[y]=min(p[x],a[i].c);
				if(bz[y])
				{
					f.push(y);
					bz[y]=false;
				}
			}
		}
		f.pop();
	}
	return dis[ed]<1061109567;
}
int flow()
{
	int ans=0;
	while(spfa())
	{
		ans+=p[ed]*dis[ed];
		for(int i=ed;i!=st;i=pos[i])
		{
			a[pre[i]].c-=p[ed];
			a[pre[i]^1].c+=p[ed];
		}
	}
	return ans;
}
void work1()
{
	st=0,ed=tot*2+1;
	len=-1;
	memset(last,-1,sizeof(last));
	for(int i=1;i<=n;i++)
		for(int j=1;j<i+m;j++)
		{
			ins(id[i][j],id[i][j]+tot,1,b[i][j]),ins(id[i][j]+tot,id[i][j],0,-b[i][j]);
			if(i==1) ins(st,id[i][j],1,0),ins(id[i][j],st,0,0);
			if(i==n) ins(id[i][j]+tot,ed,1,0),ins(ed,id[i][j]+tot,0,0);
			else{
				ins(id[i][j]+tot,id[i+1][j],1,0),ins(id[i+1][j],id[i][j]+tot,0,0);
				if(j<i+m) ins(id[i][j]+tot,id[i+1][j+1],1,0),ins(id[i+1][j+1],id[i][j]+tot,0,0);
			}
		}
	printf("%d\n",-flow());
}
void work2()
{
	st=0,ed=tot+1;
	len=-1;
	memset(last,-1,sizeof(last));
	for(int i=1;i<=n;i++)
		for(int j=1;j<i+m;j++)
		{
			if(i==1) ins(st,id[i][j],1,0),ins(id[i][j],st,0,0);
			if(i==n) ins(id[i][j],ed,INF,b[i][j]),ins(ed,id[i][j],0,-b[i][j]);
			else{
				ins(id[i][j],id[i+1][j],1,b[i][j]),ins(id[i+1][j],id[i][j],0,-b[i][j]);
				if(j<i+m) ins(id[i][j],id[i+1][j+1],1,b[i][j]),ins(id[i+1][j+1],id[i][j],0,-b[i][j]);
			}
		}
	printf("%d\n",-flow());
}
void work3()
{
	st=0,ed=tot+1;
	len=-1;
	memset(last,-1,sizeof(last));
	for(int i=1;i<=n;i++)
		for(int j=1;j<i+m;j++)
		{
			if(i==1) ins(st,id[i][j],1,0),ins(id[i][j],st,0,0);
			if(i==n) ins(id[i][j],ed,INF,b[i][j]),ins(ed,id[i][j],0,-b[i][j]);
			else{
				ins(id[i][j],id[i+1][j],INF,b[i][j]),ins(id[i+1][j],id[i][j],0,-b[i][j]);
				if(j<i+m) ins(id[i][j],id[i+1][j+1],INF,b[i][j]),ins(id[i+1][j+1],id[i][j],0,-b[i][j]);
			}
		}
	printf("%d\n",-flow());
}
int main()
{
	int t=0;
	scanf("%d %d",&m,&n);
	for(int i=1;i<=n;i++)
		for(int j=1;j<i+m;j++)
		{
			scanf("%d",&b[i][j]);
			b[i][j]=-b[i][j];
			id[i][j]=++t;
		}
	tot=(m+n+m-1)*n/2;
	work1();
	work2();
	work3();
}