2023牛客暑期多校训练营7 C

C Beautiful Sequence

 

最痛苦的一场，真真切切地感受到了坐牢的感觉。

C题看过好多人的题解，看起来像一种做法看起来又像不同的做法，而且很多都是“显然”、“易得”这些听不懂的话，我觉得这些“显然”一点也不“显然”！！！研究了一下午才学明白。

思路：

差分异或数组，A[i + 1] = A[i] ^ B[i],所以A[1]确定了的话整个序列也确定了，所以只要找A[1]的第k小就行了。

根据xor的一个性质，x ^ x = 0。

a1 ^ a2 = b1

a2 ^ a3 = b2

a3 ^ a4 = b3

将等式左右相加(或许应该叫“相异或“)可以得到a1 ^ a2 ^ a2 ^ a3 ^ a3 ^ a4 = b1 ^ b2 ^ b3，即a1 ^ a4 = b1 ^ b2 ^ b3。

于是我们可以得到一个结论：a1 ^ ai = b1 ^ ... ^ bi-1。

可以将右边的式子看作前缀和，用pre[i - 1]来表示。

根据xor的另一个性质：a ^ b = c 等价于 a = b ^ c， 可以得到a1 = ai ^ pre[i - 1]，根据题意a1和ai都是未知量，只有pre[i]是已知量。 

这样第一部分就结束了，还有一部分。

a[i] ^ a[i+1] = b[i] --> a[i] ^ b[i] = a[i + 1]

b[i]从二进制的高位到低位枚举，如果当前位为0，a[i + 1]的当前位与a[i]的当前位相同，单从当前这一位来说满足条件a[i] <= a[i + 1]， 但如果b[i]的当前位为1呢，0^1 = 1,1 ^ 1 = 0,所以a[i]和a[i + 1]一定不相同，又要保证a[i] <= a[i + 1]，所以a[i]的当前位一定为0，a[i + 1]的当前位一定为1。再往后a[i + 1]随便是0还是1都无所谓了，因为高位已经满足条件了。

到此为止第二部分也结束了，然后我们就要想办法将这两部分利用起来。

b[i]的最高非零位可以得到a[i]的这一位一定为0，而a1 = ai ^ pre[i - 1]，在这一位上ai是多少以及pre[i - 1]是多少我们都知道了，然后a1在这一位上固定是多少我们也就知道了。然后如果出现矛盾的话就输出-1。

然后将找第k小，这一步反而是最简单的，我们只要把（k - 1）的二进制数依次填进空位就可以了。

大体思路就是这样，具体代码因人而异，以下是我的版本

#include<bits/stdc++.h>
#define IOS ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
#define endl '\n'

using namespace std;

typedef pair<int, int> PII;
typedef long long ll;

const int N = 1000010;

int n, k;
int b[N], a[N], pre[N];
int w[50];

inline void solve()
{
	memset(w, -1, sizeof w);
	cin >> n >> k;
	for(int i = 1; i < n; i ++)
	{
		cin >> b[i];
		pre[i] = b[i] ^ pre[i - 1];
	}
	
	int f = 0;
	for(int i = 1; i < n; i ++)
	{
		for(int j = 29; j >= 0; j --)
		{
			if(b[i] >> j & 1)
			{
				int x = pre[i - 1] >> j & 1;
				if(w[j] == -1)
				{
					w[j] = x;
				}
				else
				{
					if(w[j] != x)f = 1;
				}
				break;
			}
		}
	}

	if(f)
	{
		cout << -1 << endl;
		return;
	}
	
	k --;
	for(int i = 0; i <= 29; i ++)
	{
		if(w[i] == -1)
		{
			if(k & 1)
			{
				w[i] = 1;
				k >>= 1;
			}
			else w[i] = 0, k >>= 1;
		}
	}

	if(k) 
	{
		cout << -1 << endl;
		return;
	}
	a[1] = 0;
	for(int i = 0; i <= 29; i ++)
	{
		a[1] += w[i] << i;
	}
	for(int i = 2; i <= n; i ++)a[i] = a[i - 1] ^ b[i - 1];
	
	for(int i = 1; i <= n; i ++)cout << a[i] << ' ';
	cout << endl;
}

int main()
{
	IOS
	int _;
	cin >> _;
	while(_ --)
	{
		solve();
	}
	
	return 0;
}

另外提醒一句：第一小的组合中a[1]不一定是0，举个例子：

1
4 1
8 7 2

这个测出来的结果是2 10 13 15。