結果

問題 No.371 ぼく悪いプライムじゃないよ
コンテスト
ユーザー rickytheta
提出日時 2016-05-14 00:11:00
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 3 ms / 1,000 ms
コード長 2,278 bytes
コンパイル時間 2,269 ms
コンパイル使用メモリ 161,008 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-12-30 18:05:33
合計ジャッジ時間 3,523 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 42
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:37:8: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   37 |   scanf("%lld%lld",&L,&H);
      |   ~~~~~^~~~~~~~~~~~~~~~~~

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef vector<int> vi;
typedef vector<ll> vl;
typedef complex<double> P;
typedef pair<int,int> pii;
#define REP(i,n) for(ll i=0;i<n;++i)
#define REPR(i,n) for(ll i=1;i<n;++i)
#define FOR(i,a,b) for(ll i=a;i<b;++i)

#define DEBUG(x) cout<<#x<<": "<<x<<endl
#define DEBUG_VEC(v) cout<<#v<<":";REP(i,v.size())cout<<" "<<v[i];cout<<endl
#define ALL(a) (a).begin(),(a).end()

#define MOD (ll)(1e9+7)
#define ADD(a,b) a=((a)+(b))%MOD
#define FIX(a) ((a)%MOD+MOD)%MOD

ll L,H;
bool isp[100000];
vl ps;

ll pow(ll a,ll b){
  ll r = 1;
  while(b){
    if(b&1)r=r*a;
    a=a*a;
    b>>=1;
  }
  return r;
}

int main(){
  scanf("%lld%lld",&L,&H);
  REP(i,100000)isp[i]=true;
  isp[0]=isp[1]=false;
  FOR(x,2,100000){
    if(!isp[x])continue;
    ps.push_back(x);
    {
      int it = x+x;
      while(it<100000){
        isp[it]=false;
        it+=x;
      }
    }
  }
  if(H-L<0){
    // brute-force
    ll mx = L-1;
    ll mfct = 0;
    FOR(x,L,H+1){
      REP(i,ps.size()){
        ll p = ps[i];
        if(x<=p)break;
        if(x%p==0){
          if(p>=mfct){
            mfct = p;
            mx = x;
          }
          break;
        }
      }
    }
    printf("%lld\n",mx);
  }else{
    ll ans = 1;
    ll mxmn = 1;
    REP(i,ps.size()){
      ll x = ps[ps.size()-1-i];
      if(x<mxmn)break;
      ll ub = H/x;
      while(ub*x>=L && ub>=x){
        ll mfct=x;
        if(ub<100000 && isp[ub]){
          mfct = min(x,ub);
        }else{
          REP(j,ps.size()){
            if(ub%ps[j]==0){
              mfct=min(mfct,ps[j]);
              break;
            }
          }
        }
        ll fct = min(mfct,x);
        if(fct>mxmn || (fct==mxmn&&ub*x>ans)){
          ans = ub*x;
          mxmn = fct;
          if(fct==x)break;
        }
        --ub;
      }
    }
    printf("%lld\n",ans);
    // set<ll> dp;
    // dp.insert(1);
    // FOR(i,0,ps.size()){
    //   ll x = ps[ps.size()-1-i];
    //   ll ans = -1;
    //   for(ll p : dp){
    //     if(p*x>H)break;
    //     dp.insert(p*x);
    //     if(p!=1 && L<=p*x){
    //       ans = max(ans,p*x);
    //     }
    //   }
    //   if(ans!=-1){
    //     printf("%lld\n",ans);
    //     return 0;
    //   }
    // }
  }
  return 0;
}
0