結果

問題 No.371 ぼく悪いプライムじゃないよ
ユーザー どらら
提出日時 2016-05-14 02:14:41
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 5 ms / 1,000 ms
コード長 1,525 bytes
コンパイル時間 962 ms
コンパイル使用メモリ 82,732 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-12-30 18:08:52
合計ジャッジ時間 2,365 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 42
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cctype>
#include <cmath>
#include <iostream>
#include <queue>
#include <list>
#include <stack>
#include <map>
#include <numeric>
#include <set>
#include <sstream>
#include <string>
#include <vector>
using namespace std;
#define REP(i,a,n) for(int i=(a); i<(int)(n); i++)
#define rep(i,n) REP(i,0,n)
#define FOR(it,c) for(__typeof((c).begin()) it=(c).begin(); it!=(c).end(); ++it)
#define ALLOF(c) (c).begin(), (c).end()
typedef long long ll;

#define MAX_PRIME 105000
bool isprime[MAX_PRIME+1]; 
vector<int> primes; 
void init_prime(){
  fill(isprime, isprime+MAX_PRIME+1, true); 
  isprime[0] = isprime[1] = false; 
  for(int i = 2; i <= MAX_PRIME; i++) { 
    if (isprime[i]) { 
      primes.push_back(i); 
      for(int j = i*2; j <= MAX_PRIME; j += i) 
      isprime[j] = false; 
    } 
  } 
}

ll findmin(ll x){
  ll d, q;
  while(x>=4 && x%2==0){ return 2; x/=2;}
  d = 3; q = x/d;
  while(q>=d){
    if(x%d==0){ return d; x=q; }else{ d+=2; }
    q=x/d;
  }
  return x;
}

int main(){
  ios::sync_with_stdio(false);
  init_prime();

  ll L, H;
  cin >> L >> H;

  for(int i=primes.size()-1; i>=0; i--){
    ll p = primes[i];
    if(H<p) continue;
    ll q = H/p;

    ll ret = -1;
    for(ll j=q+1; j>=p; j--){
      if(p*j<L) break;
      if(L<=p*j && p*j<=H){
        ll res = findmin(p*j);
        if(res == p){
          ret = max(ret, p*j);
        }
      }
    }
    if(ret > 0){
      cout << ret << endl;
      break;
    }
  }
  return 0;
}
0