結果

問題 No.1881 Everything is the same...
コンテスト
ユーザー cologne
提出日時 2022-03-18 11:59:25
言語 C++17(clang)
(17.0.6 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 2,920 bytes
コンパイル時間 4,006 ms
コンパイル使用メモリ 135,120 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-02 14:54:58
合計ジャッジ時間 6,037 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 46 WA * 6
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <algorithm>
#include <cstring>
#include <map>
#include <utility>
#include <vector>

#include <cstdio>
#include <cinttypes>
#include <unistd.h>
#include <sys/stat.h>
#include <sys/mman.h>
char *p;
struct Input
{
	Input()
	{
		struct stat st;
		fstat(0, &st);
		p = (char *)mmap(0, st.st_size, PROT_READ, MAP_SHARED, 0, 0);
	}
} inp;

inline uint32_t readU32()
{
	int t;
	uint32_t r;
	for (r = *p++ - 48; (t = *p++ - 48) >= 0; r = r * 10 + t)
		(void)0;
	return r;
}

using namespace std;

const int MAXN = 100000;
int F[MAXN + 1];
void factor_init()
{
	for (int i = 2; i <= MAXN; i++)
		if (!F[i])
		{
			F[i] = i;
			for (long j = 1L * i * i; j <= MAXN; j += i)
			{
				F[j] = i;
			}
		}
}

#include <iostream>
vector<pair<int, int>> factorize(int N)
{
	vector<pair<int, int>> ret;
	while (N != 1)
	{
		int d = F[N];
		if (ret.empty() || ret.back().first != d)
			ret.emplace_back(d, 1);
		else
			ret.back().second++;
		N /= d;
	}
	return ret;
}

vector<int> T(vector<int> A)
{
	if (A.empty())
		return {};
	vector<int> ret(*max_element(A.begin(), A.end()));
	for (int a : A)
		for (int i = 0; i < a; i++)
			ret[i]++;
	return ret;
}

vector<vector<int>> bktk(vector<int> A)
{
	int tot = 1;
	for (int a : A)
		tot *= a;
	vector<vector<int>> R;
	for (int i = 0; i < tot; ++i)
	{
		int c = i;
		vector<int> X;
		for (int a : A)
		{
			X.push_back(c % a);
			c /= a;
		}
		R.push_back(X);
	}
	return R;
}

int calc(vector<int> A)
{
	static map<vector<int>, int> memo;
	if (memo.count(A))
		return memo[A];

	vector<int> grundy;

	vector<int> cA(A.size());
	for (int i = 0; i < (int)A.size(); ++i)
		cA[i] = 1 + A[i] - (i == 0 ? 0 : A[i - 1]);
	for (auto dA : bktk(cA))
	{
		vector<int> B(A.size());
		for (int i = 0; i < (int)A.size(); ++i)
			B[i] = A[i] - dA[i];
		if (!B.empty() && B.front() == 0)
			B.erase(B.begin());
		if (B != A)
			grundy.push_back(calc(B));
	}

	sort(grundy.begin(), grundy.end());
	grundy.erase(unique(grundy.begin(), grundy.end()), grundy.end());
	for (int i = 0;; ++i)
		if (i >= (int)grundy.size() || grundy[i] != i)
			return memo[A] = i;
}

int smemo[MAXN + 1];
void solve_init()
{
	memset(smemo, -1, sizeof smemo);
}
int solve(int x)
{
	if (~smemo[x])
		return smemo[x];
	map<int, vector<int>> C;
	auto add = [&](int p, int k)
	{
		if (!C.count(p))
			C[p] = {};
		C[p].push_back(k);
	};
	for (auto [p, k] : factorize(x))
		if (p == 2)
		{
			if (k >= 2)
				add(p, 1);
			if (k >= 3)
				add(p, k - 2);
		}
		else
		{
			if (k >= 2)
				add(p, k - 1);
			for (auto [q, t] : factorize(p - 1))
				add(q, t);
		}
	vector<int> part;
	for (auto [p, v] : C)
	{
		auto Tv = T(v);
		part.insert(part.end(), Tv.begin(), Tv.end());
	}
	part = T(part);
	reverse(part.begin(), part.end());
	return smemo[x] = calc(part);
}

int main()
{
	factor_init();
	solve_init();
	int N = readU32();
	int G = 0;
	for (int i = 0; i < N; i++)
		G ^= solve(readU32());
	printf("%c\n", G ? 'Y' : 'X');
}
0