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;}