結果

問題 No.826 連絡網
ユーザー Bwambocos
提出日時 2019-05-03 22:28:04
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,297 bytes
コンパイル時間 1,831 ms
コンパイル使用メモリ 177,756 KB
実行使用メモリ 6,824 KB
最終ジャッジ日時 2024-12-31 18:11:42
合計ジャッジ時間 8,263 ms
ジャッジサーバーID
(参考情報)
judge5 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 1 WA * 1
other AC * 11 WA * 19
権限があれば一括ダウンロードができます

ソースコード

diff #

#include "bits/stdc++.h"
#define in std::cin
#define out std::cout
#define rep(i,N) for(LL i=0;i<N;++i)
typedef long long int LL;

std::map<LL, LL> decomposit_prime(LL n)
{
	LL d = 2;
	std::map<LL, LL>res;
	while (d * d <= n)
	{
		if (n % d == 0)
		{
			++res[d];
			n /= d;
		}
		else ++d;
	}
	++res[n];
	return res;
}

LL gcd(LL a, LL b)
{
	if (a < b) std::swap(a, b);
	LL r = a % b;
	while (r != 0)
	{
		a = b;
		b = r;
		r = a % b;
	}
	return b;
}

LL lcm(LL a, LL b)
{
	return a / gcd(a, b) * b;
}

int main()
{
	LL N, P;
	in >> N >> P;

	std::vector<LL>primes, nums;
	for (LL i = 2; i <= N; ++i)
	{
		bool flag = true;
		for (LL j = 2; j * j <= i; ++j)
		{
			if (i % j == 0)
			{
				flag = false;
				break;
			}
		}
		if (flag) primes.push_back(i);
	}
	auto temp = decomposit_prime(P);
	for (auto t : temp) nums.push_back(t.first);

	std::vector<LL>ok;
	for (LL i : primes)
	{
		for (LL j : nums)
		{
			if (lcm(i, j) <= N)
			{
				ok.push_back(i);
				break;
			}
		}
	}

	LL alllcm = 1;
	for (LL i : ok) alllcm = lcm(alllcm, i);
	LL ans = 0;
	rep(i, primes.size())
	{
		bool flag = false;
		for (LL j : nums)
		{
			if (lcm(j, primes[i]) <= N)
			{
				flag = true;
				break;
			}
		}
		if (!flag)
		{
			ans = N - (primes.size() - i) - 1;
			break;
		}
	}
	out << ans << std::endl;
}
0