BJOI2011 day3 casket

描述

1．星之器之中有 N×M 个区域，可看作分成 N 行和 M 列的格子，每个区域之中有若干单位的称为“星”的对象，这个对象的最小单位已经被确定，所以，这个数量总是整数。
2．魔法使可以驱动星之器中位于同一行或同一列的不相邻（有公共边的区域称为相邻的）两个区域中各 1 单位的“星”，使得它们分别向中心移动 1 格。
3．每一次使用 2 中的方法驱动“星”，将会产生魔力，魔法使会得到这一部分魔力。魔力的量等于这个两个区域之间所间隔的区域数。

在经过了进一步的研究之后，人们知道了这个星之器最初的状态（Ini）以及
最终被他们得到时的状态（Fin）。
你希望知道，星之器最多帮助它的拥有者提供了多少的魔力。即：经过一系
列上述操作由初态（Ini）变为终态（Fin），至多产生多少魔力。
需要注意的是，显然操作过程中每个区域内“星”的数量不能是负的，即：
如果那个区域已经没有“星”了，当然就不能继续操作了。

100%的数据中 1 ≤ N，M ≤ 200，Iniij，Finij ≤ 1000。
所有数据保证了至少存在一个操作方法使得星之器由初态变为终态，同时保
证了初态与终态不是完全相同的。

思路

法一:有几个数据是n=1的，我就联想到把很多行压成一行，然后模拟。列也同理。我用的就是这个方法。

法二:通过观看了《宇宙与人》，我们悟出了“有钱不去投资，再多的财富也将没有任何意义。”。所以给予每个点一个【势能】，计算出δE。

代码

#include <iostream>
using namespace std;
long long int solve(int s[],int e[],int n)
{
	long long int ans=0;
	int top,bot;
	top=0;bot=n-1;
	while (top<bot-1)
	{
		if (s[top]-e[top]<s[bot]-e[bot])
		{
			ans+=(bot-top-1)*(s[top]-e[top]);
			s[bot]-=s[top]-e[top];
			s[bot-1]+=s[top]-e[top];
			s[top+1]+=s[top]-e[top];
			top++;continue;
		}
		if (s[top]-e[top]>s[bot]-e[bot])
		{
			ans+=(bot-top-1)*(s[bot]-e[bot]);
			s[top]-=s[bot]-e[bot];
			s[top+1]+=s[bot]-e[bot];
			s[bot-1]+=s[bot]-e[bot];
			bot--;continue;
		}
		if (s[top]-e[top]==s[bot]-e[bot])
		{
			ans+=(bot-top-1)*(s[bot]-e[bot]);
			s[top+1]+=s[bot]-e[bot];
			s[bot-1]+=s[bot]-e[bot];
			bot--;top++;continue;
		}
	}
	return ans;
}
int main()
{
	int n,m;
	long long ret=0;
	scanf("%d%d",&n,&m);
	static int ins[201][201],ine[201][201],te[201],ts[201];
	for (int i=0;i<n;i++)
		for (int j=0;j<m;j++)
			scanf("%d",&ins[i][j]);
	for (int i=0;i<n;i++)
		for (int j=0;j<m;j++)
			scanf("%d",&ine[i][j]);
	memset(ts,0,sizeof(ts));
	memset(te,0,sizeof(te));
	for (int i=0;i<n;i++)
		for (int j=0;j<m;j++)
			ts[i]+=ins[i][j],te[i]+=ine[i][j];
	ret+=solve(ts,te,n);
	memset(ts,0,sizeof(ts));
	memset(te,0,sizeof(te));
	for (int j=0;j<m;j++)
		for (int i=0;i<n;i++)
			ts[j]+=ins[i][j],te[j]+=ine[i][j];
	ret+=solve(ts,te,m);
	printf("%I64d\n",ret);
	return 0;
}