結果

問題 No.371 ぼく悪いプライムじゃないよ
ユーザー fiord
提出日時 2016-07-16 18:23:33
言語 C++11(廃止可能性あり)
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 1,221 bytes
コンパイル時間 532 ms
コンパイル使用メモリ 67,720 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-15 13:49:58
合計ジャッジ時間 1,956 ms
ジャッジサーバーID
(参考情報)
judge4 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 4
other AC * 40 WA * 2
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <iostream>
#include <cmath>
#include <vector>
#include <algorithm>
using namespace std;
typedef long long ll;
vector<ll> p;
bool isprime[100001];
ll findprime(ll a) {
	for (int i = 0; i < p.size()&&p[i]*p[i]<=a; i++) {
		if (a%p[i] == 0)	return p[i];
	}
	return a;
}
int main() {
	for (int i = 0; i <= 100000; i++)	isprime[i] = true;
	isprime[0] = isprime[1] = false;
	for (int i = 2; i <= 100000; i++) {
		if (!isprime[i])	continue;
		p.push_back(i);
		for (int j = 2 * i; j <= 100000; j += i)	isprime[j] = false;
	}
	ll l, h;	cin >> l >> h;
	ll ans = 1, out = 1;
	//bigを一番答えの中の大きな素数と仮定。l<=another*big<=hを満たす範囲でanotherを回す
	//その中でanotherを素因数分解し、一番小さな素数でansを更新
	//big<ansのとき、それ以上の更新はありえない(残りの因数は全てbig以下のため)
	for (ll big = sqrt(h); big >= ans; big--) {
		if (!isprime[big])	continue;
		for (ll another = h / big; another*big >= l; another--) {
			ll small=findprime(another);
			if (small > ans) {
				ans = small;
				out = another*big;
			}
			else if (small == ans) {
				out = max(out, another*big);
			}
		}
	}
	cout << out << endl;
	return 0;
}
0