[hdu 2686]Matrix

网上说这道题的题解是费用流

我粗粗看了一下数据范围，觉得出题者似乎是让我们用 “大（d）屁（p）” 的样子，为了尊重出题人，我还是谢谢吧喵~

首先，一条回路可以看做是两条路齐头并进，这是 大屁 和 费永柳 的共同思想

故我们可用一条条次对角线（即 x+y=k  的对角线）来划分状态

考虑枚举对角线上两个点 x1，y1，x2，y2
很大力地用 f[x1][y1][x2][y2] 表示两条路分别走到 点（x1, y1）和 点（x2, y2）时的最佳方案

这样 f[x1][y1][x2][y2] 就可以由 f[x1-1][y1][x2-1][y2], f[x1-1][y1][x2][y2-1], f[x1][y1-1][x2-1][y2], f[x1][y1-1][x2][y2-1] 来得到（丧心病狂）

还有当 点（x1, y1）和 点（x2, y2）在对角线上相邻时情况特殊，因为 点（x1-1, y1）和 点（x2, y2-1）重合了，需要特别考虑（真是丧心病狂）

虽然看上去很恶心，但我觉得再怎么说短也是一大优势，你的费用流有这么短喵？

#include <cstdio>
#include <cstring>
#define max(x, y) ((x)>(y) ? (x):(y))
const int size=32;

int n;
int a[size][size];
int f[size][size][size][size];
inline int getint();
inline void putint(int);

int main()
{
	while (scanf("%d", &n)!=EOF)
	{
		for (int i=1;i<=n;i++)
			for (int j=1;j<=n;j++)
				a[i][j]=getint();
		memset(f, 0, sizeof f);
		f[2][1][1][2]=a[1][1]+a[1][2]+a[2][1];
		for (int i=4;i<2*n;i++)
			for (int y1=i<=n+1?1:i-n;y1<i-1 && y1<n;y1++)
				for (int y2=y1+1;y2<i && y2<=n;y2++)
				{
					int x1=i-y1, x2=i-y2;
					if (y2-y1==1) f[x1][y1][x2][y2]=max(max(f[x1][y1-1][x2][y1], f[x2][y1][x2-1][y2]), f[x1][y1-1][x2-1][y2]);
					else f[x1][y1][x2][y2]=max(max(f[x1-1][y1][x2-1][y2], f[x1-1][y1][x2][y2-1]), max(f[x1][y1-1][x2-1][y2], f[x1][y1-1][x2][y2-1]));
					f[x1][y1][x2][y2]+=a[x1][y1]+a[x2][y2];
				}
		putint(a[n][n]+f[n][n-1][n-1][n]);
	}

	return 0;
}
inline int getint()
{
	register int num=0;
	register char ch;
	do ch=getchar(); while (ch<‘0‘ || ch>‘9‘);
	do num=num*10+ch-‘0‘, ch=getchar(); while (ch>=‘0‘ && ch<=‘9‘);
	return num;
}
inline void putint(int num)
{
	char stack[15];
	register int top=0;
	for ( ;num;num/=10) stack[++top]=num%10+‘0‘;
	for ( ;top;top--) putchar(stack[top]);
	putchar(‘\n‘);
}