結果

問題 No.371 ぼく悪いプライムじゃないよ
コンテスト
ユーザー chocorusk
提出日時 2018-09-29 00:01:42
言語 C++11
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=gnu++11 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 288 ms / 1,000 ms
コード長 1,468 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 630 ms
コンパイル使用メモリ 112,928 KB
実行使用メモリ 178,688 KB
最終ジャッジ日時 2026-04-28 19:24:34
合計ジャッジ時間 3,684 ms
ジャッジサーバーID
（参考情報） 		judge2_1 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 42
権限があれば一括ダウンロードができます
コンパイルメッセージ
In file included from /home/linuxbrew/.linuxbrew/Cellar/gcc/15.2.0_1/include/c++/15/ostream:42,
                 from /home/linuxbrew/.linuxbrew/Cellar/gcc/15.2.0_1/include/c++/15/iostream:43,
                 from main.cpp:3:
In member function 'std::basic_ostream<_CharT, _Traits>::__ostream_type& std::basic_ostream<_CharT, _Traits>::operator<<(long long int) [with _CharT = char; _Traits = std::char_traits<char>]',
    inlined from 'int main()' at main.cpp:84:8:
/home/linuxbrew/.linuxbrew/Cellar/gcc/15.2.0_1/include/c++/15/bits/ostream.h:212:25: warning: 'ans' may be used uninitialized [-Wmaybe-uninitialized]
  212 |       { return _M_insert(__n); }
      |                ~~~~~~~~~^~~~~
main.cpp: In function 'int main()':
main.cpp:76:18: note: 'ans' was declared here
   76 |         ll mx=0, ans;
      |                  ^~~

ソースコード

diff #
raw source code 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
using namespace std;
typedef long long int ll;
typedef pair<ll, int> P;

const ll MAX=110001;
 
vector<ll> prime;
bool isprime[MAX];
 
void sieve(){
	for(ll i=3; i<MAX; i+=2){
		isprime[i]=1;
	}
	isprime[2]=1;
	prime.push_back(2);
	for(ll i=3; i<MAX; i++){
		if(isprime[i]){
			prime.push_back(i);
			for(ll j=2*i; j<MAX; j+=i){
				isprime[j]=0;
			}
		}
	}
	return;
}

int main()
{
	ll l, h;
	cin>>l>>h;
	sieve();
	int i1=0, i2=prime.size()-1;
	while(i1!=i2){
		int i0=(i1+i2+1)/2;
		if(prime[i0]*prime[i0]<=h){
			i1=i0;
		}else{
			i2=i0-1;
		}
	}
	if(l<=prime[i1]*prime[i1]){
		cout<<(*(upper_bound(prime.begin(), prime.end(), h/prime[i1])-1))*prime[i1]<<endl;
		return 0;
	}
	vector<P> fac[5000000];
	ll x[5000000];
	for(ll i=l; i<=h; i++) x[i-l]=i;
	for(auto p:prime){
		ll x0=(l+p-1)/p*p;
		for(ll x1=x0; x1<=h; x1+=p){
			int e=0;
			while(x[x1-l]%p==0){
				x[x1-l]/=p;
				e++;
			}
			fac[x1-l].push_back(P(p, e));
		}
	}
	for(ll i=l; i<=h; i++){
		if(x[i-l]>1) fac[i-l].push_back(P(x[i-l], 1));
	}
	ll mx=0, ans;
	for(ll i=h; i>=l; i--){
		if(fac[i-l].size()==1 && fac[i-l][0].second==1) continue;
		if(fac[i-l][0].first>mx){
			mx=fac[i-l][0].first;
			ans=i;
		}
	}
	cout<<ans<<endl;
	return 0;
}
0