結果

問題 No.990 N×Mマス計算（Kの倍数）
コンテスト
ユーザー MarcusAureliusAntoninus
提出日時 2020-02-14 23:02:59
言語 C++17
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=c++17 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 164 ms / 2,000 ms
コード長 1,260 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,780 ms
コンパイル使用メモリ 225,204 KB
実行使用メモリ 23,936 KB
最終ジャッジ日時 2026-06-08 19:16:48
合計ジャッジ時間 3,424 ms
ジャッジサーバーID
（参考情報） 		judge1_0 / judge2_1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>

int64_t calcGCD(int64_t a, int64_t b)
{
	while (a)
	{
		b %= a;
		std::swap(a, b);
	}
	return b;
}

int main()
{
	int N, M;
	int64_t K;
	char op;
	scanf("%d%d%lld %c", &N, &M, &K, &op);
	std::vector<int64_t> A(N), B(M);
	for (auto& e: B)
		scanf("%lld", &e);
	for (auto& e: A)
		scanf("%lld", &e);

	int64_t ans{};
	if (op == '+')
	{
		std::map<int64_t, int64_t> rest_a, rest_b;
		for (auto& e: A)
			rest_a[e % K]++;
		for (auto& e: B)
			rest_b[e % K]++;
		for (auto& e: rest_a)
			ans += e.second * rest_b[K - e.first];
		ans += rest_a[0] * rest_b[0];
	}
	else
	{
		std::vector<int> div, count;
		div.reserve(10'000);
		for (int64_t i{1}; i * i <= K; i++)
			if (K % i == 0)
			{
				div.push_back(i);
				if (i * i != K)
					div.push_back(K / i);
			}
		std::sort(div.begin(), div.end());
		count.resize(div.size());
		
		for (auto& e: B)
			count[std::lower_bound(div.begin(), div.end(), K / calcGCD(e, K)) - div.begin()]++;
		for (int i{(int)div.size() - 1}; i >= 0; i--)
			for (int j{i + 1}; j < (int)div.size(); j++)
				if (div[j] % div[i] == 0)
					count[j] += count[i];
		
		for (auto& e: A)
			ans += count[std::lower_bound(div.begin(), div.end(), calcGCD(e, K)) - div.begin()];
	}
	printf("%lld\n", ans);

	return 0;
}
0