BZOJ 2321 星器 能量守恒定律

2321: [BeiJing2011集训]星器

Time Limit: 1 Sec Memory Limit: 128 MB
Description

Magic Land上的时间又过了若干世纪……
现在，人们谈论着一个传说：从前，他们的祖先来到了一个位于东方的岛屿，那里简直就是另外一个世界。善于分析与构造的Magic Land上的人们总是不明白那里的人们是如何不借助精确的实验与计算驱动和操纵魔法。
偶然地，一个魔法使（Magician）来到了Magic Land，在临走的时候留下了一个神奇的盒子，叫做星器（Casket of star）。
虽然不知道这个盒子是做什么的，但是经过了大量的实验和计算后，人们已经清楚它的一些事实：
1．星器之中有N×M个区域，可看作分成N行和M列的格子，每个区域之中有若干单位的称为“星”的对象，这个对象的最小单位已经被确定，所以，这个数量总是整数。
2．魔法使可以驱动星器中位于同一行或同一列的不相邻（有公共边的区域称为相邻的）两个区域中各1单位的“星”，使得它们分别向中心移动1格。
3．每一次使用2中的方法驱动“星”，将会产生魔力，魔法使会得到这一部分魔力。魔力的量等于这个两个区域之间所间隔的区域数。
这样，我们可以用一个N×M的数表来表示星器的状态，比如N=2，M=3时：

2 0 1
5 1 4

1 2 0
5 1 4

当星器为上图的状态时，通过操纵第一行的第1和3个区域中的“星”（加粗的数字对应的区域），变为右图所示的状态，同时，将产生1单位的魔力（因为这两个区域之间恰好隔了1个区域）。
在经过了进一步的研究之后，人们知道了这个星器最初的状态（Ini）以及最终被他们得到时的状态（Fin）。
你希望知道，星器最多帮助它的拥有者提供了多少的魔力。即：经过一系列上述操作由初态（Ini）变为终态（Fin），至多产生多少魔力。
需要注意的是，显然操作过程中每个区域内“星”的数量不能是负的，即：如果那个区域已经没有“星”了，当然就不能继续操作了。

Input

第一行包含两个正整数N、M表示星器的大小。
接下来的N行，每行包含M个自然数：Iniij，描绘了初态（Ini）。
在一个空行后的N行，每行包含M个自然数：Finij，描绘了终态（Fin）。
Output

输出一个正整数，表示至多产生的魔力。
Sample Input

【输入样例1】

5 5
1 0 0 0 1
0 0 0 0 0
0 0 0 0 0
0 1 0 1 1
1 0 0 0 0
0 0 0 0 0
0 0 0 0 1
2 0 0 0 1
0 0 2 0 0
0 0 0 0 0

【输入样例2】

1 4

10 20 30 40

0 0 100 0

Sample Output

【输出样例1】

7

【样例1解释】

唯一的一种操作方法是：
对第5列的两个“星”进行一次操作，产生魔力2；
对第1列的两个“星”进行两次操作，产生魔力3+1；
对第4行的两个“星”进行一次操作，产生魔力1；
一共产生7单位的魔力。

【输出样例2】

50

HINT

【数据规模和约定】

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

题目大意：给定一个矩阵，定义一个操作：
选择两个同一行或同一列不相邻的点，将这两个点上各一个星向中间移动一位，产生魔力为两点间距离-1，求始态到终态的产生魔力

思路：
摘自 < popoqqq > 链接
定义一个星的势能为这个点到原点的欧几里得距离的平方
即一个在(i,j)位置上的星的势能为i*i+j*j
假如一次操作之前两个星的位置为(i,j)和(i,k)，其中j+2<=k
那么操作之前两个星的势能和为i*i+j*j+i*i+k*k
操作后两个星的位置为(i,j+1)和(i,k-1)
势能和为i * i+(j+1) * (j+1)+i * i+(k-1) * (k-1)
E前-E后=2 * k-2 * j-2=2 * (k-j-1)
而k-j-1恰好是此次操作所释放的势能
于是我们计算初末势能之差，除以2即是答案

#include <cstdio>  
#include <cstring>  
#include <iostream>  
#include <algorithm>  
#define LL long long
using namespace std;  

int m, n;  
LL E1, E2;  

int main(){  
    scanf("%d%d", &m, &n);  
    for(int i=1; i<=m; i++){
        for(int j=1; j<=n; j++){
            int num; scanf("%d", &num);
            E1 += (i * i + j * j) * num;  
        }  
    }
    for(int i=1; i<=m; i++){
        for(int j=1; j<=n; j++){
            int num; scanf("%d", &num);
            E2 += (i * i + j * j) * num;  
        }  
    }
    cout << ( (E1 - E2) >> 1 ) << endl;  
    return 0;
}