洛谷P3179 [HAOI2015] 数组游戏

 [HAOI2015] 数组游戏

 题目描述

有一个长度为 n 的数组，甲乙两人在上面进行这样一个游戏：首先，数组上有一些格子是白的，有一些是黑的。然后两人轮流进行操作。

每次操作选择一个白色的格子，假设它的下标为 x。接着，选择一个大小在 1,2,...,n/x 之间的整数 k，然后将下标为x,2×x,…,k×x的格子都进行颜色翻转。不能操作的人输。

现在甲（先手）有一些询问。每次他会给你一个数组的初始状态，你要求出对于这种初始状态他是否有必胜策略。

 输入格式

第一行包含一个整数 n，表示数组的长度。

第二行包含一个整数 k，表示询问的个数。

接下来 2k 行，每两行表示一次询问。

在这两行中，第一行一个正整数 w，表示数组中有多少个格子是白色的，第二行则有w 个1 到n 之间的正整数，表示白色格子的对应下标。

## 输出格式

对于每个询问，输出一行一个字符串，若先手必胜输出 `Yes`，否则输出 `No`。

## 样例 #1

### 样例输入 #1
3
2
2
1 2
2
2 3
```

### 样例输出 #1

```
Yes
No
```

## 提示

#### 样例输入输出 1 解释

在第一个询问中，甲选择点1，然后将格子1×1和2 x 1 翻过来即可。

第二个询问中，无论甲选择哪个点，都只能翻掉一个格子。乙只需翻掉另一个格子就行了。

#### 数据规模与约定

对于100% 的数据，保证1 <= n <= 1e9，1 <= k,w <= 100 ，不会有格子在同一次询问中多次出现。

---

思路

看似起了这个好玩的题目，但本质是有一点点翻棋子游戏的意思的。

翻棋子游戏的解法是将初始状态的SG值异或和求出来，0则必败，反之必胜；

首先棋子对应的一个数，其SG值就是这个数的SG值

此时就要介绍一下SG值了

我们先假设一个状态的后继状态的SG值构成的集合为S，那么这个状态的SG值为mex{S}，其中mex{S}是指不属于集合S的最小自然数。我们现在要做的就是求每个数的SG值，那么就得先找出x的后继状态。

我们再次假设只有x是白色的。于是，我们可以选择翻x,2x,...,kx。翻完后x变黑，其他的会变白。

根据翻棋子游戏的解法，此时这个后继状态的SG值为现在的白色棋子的SG值异或和。由于只翻x时没有任何棋子是白色，故此后继状态SG值为0。。其他后继状态的SG值分别为SG(2x),SG(2x)xorSG(3x),...，故我们得到当前的Sx = {0,SG(2x),SG(2x)xorSG(3x),...},由于2x,3x均大于x，因此我们要从大往小推SG值。

下面对于floor(n/x)相同的x有SG(x)相同，假设t>x对于上述假设成立。

我们知道了Sx = {0,SG(2x),SG(2x)xorSG(3x),...}，对于floor(n/x) = floor(n/y)的x,y，容易发现Sx,Sy的值均为floor(n/x)，对于Sx,Sy中各个元素，由于，由于floor(n/x) = floor(n/y)时也会有floor(n/(tx))=floor(floor(n/x)/t)=floor(n/(ty))=floor(floor(n/y)/t),所以Sx与Sy完全相同，故mex{Sx}=mex{Sy};即SG(x)=SG(y).

又因为在tx中有一些满足floor(n/(x*t1))=floor(n/(x*t2)),所以再对t进行整除分块，设对于范围Li <=t<=Ri内的floor(n/(t*x))相同，则我们只需将此状态的SG值与当前维护的前缀异或和取异或得到的结果标记入Sx ，然后再判断区间长度奇偶性，奇数则将前缀异或和异或上此状态的SG值，否则前缀异或和不变。

我们对于>sqrt(n)​的x将其SG值存储在另一个数组的floor(n/x)位置就可以只用数组解决问题了。预处理出SG就可以求解了。这个最好开O2。

上AC代码

#include <iostream>
#include <bits/stdc++.h>
//#define register register
int n,t,w,a;
int root,cnt,res;
int sg[2][400007],rt[400007],mrk[400007];
inline int SG(int x){
	return ((x = n / (n / x)) > root)?sg[1][n / x]:sg[0][x];
}
inline void A(){
	for(register int i = 1,j;i <= n;i = j + 1){
		j = n / (n / i);
		rt[++ cnt] = j;
	}
	++ cnt;
	while(-- cnt){
		int x = rt[cnt],psg = 0,rsg = 1;mrk[psg] = cnt;
		for(register int i = x + x,j;i <= n;i = j + x){
			j = n / (n / i) / x * x ,mrk[psg ^ SG(j)] = cnt;
			((j - i) / x & 1 ^ 1) && (psg ^= SG(j));
		}
		while(mrk[rsg] == cnt) ++ rsg;
		(x>root)?sg[1][n / x]=rsg:sg[0][x] = rsg;
	}
}
int main(){
//	ios::sync_with_stdio(false);
//	cin.tie(0);
	scanf("%d %d",&n,&t);
	while(root * root <= n) ++ root;-- root;
	A();
	while(t --){
		scanf(" %d",&w);
		res = 0;
		while(w --) scanf(" %d",&a),res ^= SG(a);
		printf("%s\n",res?"Yes":"No");
	}
	
	return 0;
}