【BZOJ】2321: [BeiJing2011集训]星器(数学+特殊的技巧)

完全想不到。。

第一眼以为是爆搜,看到数据范围果断放弃;第二眼以为是网络流(因为只有行列操作,我们将起点向重点连边然后设置容量等等麻烦的东西再跑最大流),觉得建模麻烦放弃

数学。。。

首先本题有个性质:

  答案与移动方法无关(因为有行列限制并且是2个同时在同行同列移动,这样,只要这行这列有个点是终点,无论是怎么累计起来的,能到达的点到这里的距离和相等)

这样就使得我们可以绕开答案。

我们来看矩形

0 0 0 1

0 0 0 0

1 0 0 0

我们要从左下角到右上角(本题的话一定左下角要到达左上角或右下角,然后再到右上角)

距离是一样的。。即这样合并的答案一样。

在矩形中,对角线的距离是不会变的。。。

而且我们要充分挖掘它的性质。

我们将这个矩形放到每个点来看,原矩形的左上角到每个点的矩阵

然后有网上的题解:

首先我们假设两个点(i,j),(i,k)向中间移动一格,且k>j+1,那么我们可以获得的价值为k-j,这样,我们定义每个点的每个星的能量为 a[(i,j)]=i*i+j*j,这样这两个点开始的能量为i*i+j*j+i*i+k*k,移动之后,两个点变为(i,j+1),(i,k-1),这 时的能量为i*i+(j+1)*(j+1)+i*i+(k-1),这时的能量差为2*k-2*j,为获得价值的2倍,因为对于所有的价值获得都可以采用这样的方法,所以我们可以算出开始局面的能量和,结束局面的能量和,相减>>1就是答案。

这样做的原因就是我上边说的。

#include <cstdio>
#include <cstring>
#include <cmath>
#include <string>
#include <iostream>
#include <algorithm>
using namespace std;
#define rep(i, n) for(int i=0; i<(n); ++i)
#define for1(i,a,n) for(int i=(a);i<=(n);++i)
#define for2(i,a,n) for(int i=(a);i<(n);++i)
#define for3(i,a,n) for(int i=(a);i>=(n);--i)
#define for4(i,a,n) for(int i=(a);i>(n);--i)
#define CC(i,a) memset(i,a,sizeof(i))
#define read(a) a=getint()
#define print(a) printf("%d", a)
#define dbg(x) cout << #x << " = " << x << endl
#define printarr(a, n, m) rep(aaa, n) { rep(bbb, m) cout << a[aaa][bbb]; cout << endl; }
inline const int getint() { int r=0, k=1; char c=getchar(); for(; c<‘0‘||c>‘9‘; c=getchar()) if(c==‘-‘) k=-1; for(; c>=‘0‘&&c<=‘9‘; c=getchar()) r=r*10+c-‘0‘; return k*r; }
inline const int max(const int &a, const int &b) { return a>b?a:b; }
inline const int min(const int &a, const int &b) { return a<b?a:b; }

int main() {
	int ans=0, n, m;
	read(n); read(m);
	for1(i, 1, n) for1(j, 1, m) ans+=(i*i+j*j)*getint();
	for1(i, 1, n) for1(j, 1, m) ans-=(i*i+j*j)*getint();
	print(ans>>1);
	return 0;
}

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

   
1


2


0


5


1


4

   
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。

所有数据保证了至少存在一个操作方法使得星器由初态变为终态,同时保证了初态与终态不是完全相同的。

Source

时间: 2024-11-05 23:34:18

【BZOJ】2321: [BeiJing2011集训]星器(数学+特殊的技巧)的相关文章

BZOJ 2321 BeiJing2011集训 星器

题目大意:给定一个矩阵,定义一个操作: 选择两个同一行或同一列不相邻的点,将这两个点上各一个星向中间移动一位,产生魔力为两点间距离-1,求始态到终态的产生魔力 定义一个星的势能为这个点到原点的欧几里得距离的平方 即一个在(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

【BZOJ】2321: [BeiJing2011集训]星器

Description Magic Land上的时间又过了若干世纪…… 现在,人们谈论着一个传说:从前,他们的祖先来到了一个位于东方的岛屿,那里简直就是另外一个世界.善于分析与构造的Magic Land上的人们总是不明白那里的人们是如何不借助精确的实验与计算驱动和操纵魔法. 偶然地,一个魔法使(Magician)来到了Magic Land,在临走的时候留下了一个神奇的盒子,叫做星器(Casket of star). 虽然不知道这个盒子是做什么的,但是经过了大量的实验和计算后,人们已经清楚它的一些

Bzoj2321 [BeiJing2011集训]星器

Time Limit: 1 Sec  Memory Limit: 128 MBSubmit: 479  Solved: 310 Description Magic Land上的时间又过了若干世纪…… 现在,人们谈论着一个传说:从前,他们的祖先来到了一个位于东方的岛屿,那里简直就是另外一个世界.善于分析与构造的Magic Land上的人们总是不明白那里的人们是如何不借助精确的实验与计算驱动和操纵魔法. 偶然地,一个魔法使(Magician)来到了Magic Land,在临走的时候留下了一个神奇的盒

BZOJ 1013: [JSOI2008]球形空间产生器sphere

二次联通门 : BZOJ 1013: [JSOI2008]球形空间产生器sphere /* BZOJ 1013: [JSOI2008]球形空间产生器sphere 高斯消元 QAQ SB的我也能终于能秒题了啊 设球心的坐标为(x,y,z...) 那么就可以列n+1个方程,化化式子高斯消元即可 */ #include <cstdio> #include <iostream> #include <cstring> #define rg register #define Max

BZOJ 4196: [Noi2015]软件包管理器

二次联通门 : BZOJ 4196: [Noi2015]软件包管理器 /* BZOJ 4196: [Noi2015]软件包管理器 树链剖分 树链剖分 , 用线段树维护链上的和与子树的和 对于1操作, 每次预先查询待修改点x点的深度 后查询1到x这条链上有多少已安装的, 然后上面的值减下面的值 对于2操作,每次查询以x为根的子树的大小 修改整颗子树即可.. 说白了其实就是板子题 */ #include <cstdio> #define Max 500090 inline int min (int

代码收藏系列--jquery--筛选器、事件绑定技巧

Jquery筛选器的一些常用技巧,比如过滤属性等 /* 过滤获取没有含data-xsui-grid-colspan的节点 */$(this).find('.xsui-grid-cell:not([data-xsui-grid-colspan])'); jQuery绑定input控件的回车事件 // input文本输入框按确定按键后,触发Tab按键,自动跳到下一个输入框 $('input[type="text"]').bind('keydown', function (e) { if (

BZOJ 2462: [BeiJing2011]矩阵模板

2462: [BeiJing2011]矩阵模板 Time Limit: 2 Sec  Memory Limit: 128 MBSubmit: 915  Solved: 432[Submit][Status][Discuss] Description 给定一个M行N列的01矩阵,以及Q个A行B列的01矩阵,你需要求出这Q个矩阵哪些在原矩阵中出现过.    所谓01矩阵,就是矩阵中所有元素不是0就是1. Input 输入文件的第一行为M.N.A.B,参见题目描述. 接下来M行,每行N个字符,非0即1

BZOJ 1013 [JSOI2008]球形空间产生器sphere 【高斯消元】

Description 有一个球形空间产生器能够在n维空间中产生一个坚硬的球体.现在,你被困在了这个n维球体中,你只知道球面上n+1个点的坐标,你需要以最快的速度确定这个n维球体的球心坐标,以便于摧毁这个球形空间产生器. HINT 1<=n<=10 提示:给出两个定义:1. 球心:到球面上任意一点距离都相等的点.2. 距离:设两个n为空间上的点A, B的坐标为(a1, a2, …, an), (b1, b2, …, bn),则AB的距离定义为:dist = sqrt( (a1-b1)^2 +

bzoj 1923 [Sdoi2010]外星千足虫(高斯消元+bitset)

1923: [Sdoi2010]外星千足虫 Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 634  Solved: 397[Submit][Status][Discuss] Description Input 第一行是两个正整数 N, M. 接下来 M行,按顺序给出 Charles 这M次使用“点足机”的统计结果.每行 包含一个“01”串和一个数字,用一个空格隔开.“01”串按位依次表示每只虫 子是否被放入机器:如果第 i 个字符是“0”则代表编号