P1006传纸条

题解：
本题（昆虫般的窸窣声）既然要求一个最大值，并且因为是路径式，我们不能使用贪心通过此题，考虑dp；
显然，这是个矩阵dp，那么回忆下矩阵dp通式：

dp[i][j]=max(dp[i][k])+决策

对于此题，我们可以自然的想到正反同时走，而n的范围并不大（n<=50），考虑高复杂度（n⁴）,得出最终状态转移方程：

dp[i][j][p][q]=max(dp[i-1][j][p-1][q],max(dp[i-1][j][p][q-1],max(dp[i][j-1][p-1][q],dp[i][j-1][p][q-1])))+a[i][j]+a[p][q];

考虑优化
对于每个点，由于只能向对角线方向移动，则x+y=step-1
那么可得出新的状态转移方程，时间复杂度较低（n³）：

f[k][i][j] = max_ele(f[k - 1][i][j],f[k - 1][i - 1][j - 1],f[k - 1][i][j - 1],f[k - 1][i - 1][j]) + a[i][k - i + 1] + a[j][k - j + 1];

考虑最短路
最大费用最大流
我也不会
代码：

//n四次方
/*#include<bits/stdc++.h>
using namespace std;
int n,m,a[105][105],dp[55][55][55][55];
int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			scanf("%d",&a[i][j]);
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			for(int p=1;p<=n;p++){
				for(int q=1;q<=m;q++){
					dp[i][j][p][q]=max(dp[i-1][j][p-1][q],max(dp[i-1][j][p][q-1],max(dp[i][j-1][p-1][q],dp[i][j-1][p][q-1])))+a[i][j]+a[p][q];
					if(i==p&&j==q)
						dp[i][j][p][q]-=a[i][j];
				}
			}
		}
	}
	printf("%d",dp[n][m][n][m]);
	return 0;	
}*/
//n三次方
#include<bits/stdc++.h>
#define N 50020
#define INF 2147483647
const int maxn = 55;
using namespace std;
/*int f[2 * maxn][maxn][maxn];
int a[maxn][maxn];
int n,m;
int max_ele(int a,int b,int c,int d)
{
	if(b > a)
    	a = b;
	if(c > a)
		a = c;
	if(d > a)
		a = d;
	return a;
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			scanf("%d",&a[i][j]);
	for(int k = 1;k <= n + m - 1;k++)
	{
		for(int i = 1;i <= n;i++)
		{
			for(int j = 1;j <= n;j++)
			{
				if(k - i + 1 < 1 || k - j + 1 < 1)
				{
                    continue;
                }
				f[k][i][j] = max_ele(f[k - 1][i][j],f[k - 1][i - 1][j - 1],f[k - 1][i][j - 1],f[k - 1][i - 1][j]) + a[i][k - i + 1] + a[j][k - j + 1];
                if (i == j)
                {
                    f[k][i][j] -= a[i][k - i + 1];
                }
            }
        }
    }
    printf("%d",f[n + m - 1][n][n]);
    return 0;
}*/
//最大费用最大流
const int S=50001;
const int T=50002;
inline int read(){
    int x=0,f=1;char c;
    do c=getchar(),f=c=='-'?-1:f; while(!isdigit(c));
    do x=(x<<3)+(x<<1)+c-'0',c=getchar(); while(isdigit(c));
    return x*f;
}
queue<int>q;
int n,m,x,cost,top=1;
int d[N],fir[N];
bool b[N];
struct Edge{
    int to,nex,k,v;
    Edge(int _=0,int __=0,int ___=0,int ____=0):to(_),nex(__),k(___),v(____){}
}nex[250005];
inline int GetNum(int x,int y){
    return (x-1)*n+y;
}
inline bool spfa(){
    for(int i=0;i<=n*m*2;i++)
        d[i]=-INF,b[i]=false;
    d[T]=-INF;
    d[S]=0;q.push(S);
    while(!q.empty()){
        int x=q.front();q.pop();
        b[x]=false;
        for(int i=fir[x];i;i=nex[i].nex)
            if(nex[i].k && d[nex[i].to]<d[x]+nex[i].v){///最大费用
                d[nex[i].to]=d[x]+nex[i].v;
                if(!b[nex[i].to]) b[nex[i].to]=true,q.push(nex[i].to);
            }
    }
    return d[T]!=-INF;
}
int dfs(int x,int v){
    if(x==T || !v){
        cost=cost+v*d[T];
        return v;
    }
    b[x]=true;
    int tmp=0;
    for(int i=fir[x];i;i=nex[i].nex)
        if(!b[nex[i].to] && d[nex[i].to]==d[x]+nex[i].v && nex[i].k){
            int f=dfs(nex[i].to,min(v,nex[i].k));
            v-=f;nex[i].k-=f;nex[i^1].k+=f;tmp+=f;
            if(!v) break;
        }
    if(!tmp) d[x]=-INF;
    return tmp;
}
inline void Dinic(){
    while(spfa()) dfs(S,INF);
}
inline void add(int x,int y,int k,int v){
    nex[++top]=Edge(y,fir[x],k,v);
    fir[x]=top;
}
int main(){
    m=read();n=read();
    for(int i=1;i<=m;i++)
        for(int j=1;j<=n;j++){
            x=read();
            add(GetNum(i,j),GetNum(i,j)+n*m,1,x);
            add(GetNum(i,j)+n*m,GetNum(i,j),0,-x);
            add(GetNum(i,j),GetNum(i,j)+n*m,INF,0);
            add(GetNum(i,j)+n*m,GetNum(i,j),0,0);
            if(i+1<=m){
                add(GetNum(i,j)+m*n,GetNum(i+1,j),INF,0);
                add(GetNum(i+1,j),GetNum(i,j)+n*m,0,0);
            }
            if(j+1<=n){
                add(GetNum(i,j)+m*n,GetNum(i,j+1),INF,0);
                add(GetNum(i,j+1),GetNum(i,j)+n*m,0,0);
            }
        }
    add(S,GetNum(1,1),2,0);
	add(GetNum(1,1),S,0,0);
    add(GetNum(m,n)+n*m,T,INF,0);
	add(T,GetNum(m,n)+n*m,0,0);
    Dinic();
    printf("%d",cost);
	return 0;
}

亿些问题：
1.当点重合时应减去点权
2.n³做法应考虑边界