結果

問題 No.371 ぼく悪いプライムじゃないよ
ユーザー tskrextskrex
提出日時 2016-07-05 04:58:14
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 35 ms / 1,000 ms
コード長 2,799 bytes
コンパイル時間 851 ms
コンパイル使用メモリ 98,328 KB
実行使用メモリ 5,004 KB
最終ジャッジ日時 2023-08-29 01:08:48
合計ジャッジ時間 3,197 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 13 ms
4,688 KB
testcase_01 AC 14 ms
4,692 KB
testcase_02 AC 13 ms
4,788 KB
testcase_03 AC 13 ms
4,748 KB
testcase_04 AC 13 ms
4,688 KB
testcase_05 AC 13 ms
4,700 KB
testcase_06 AC 13 ms
4,784 KB
testcase_07 AC 13 ms
4,640 KB
testcase_08 AC 13 ms
4,692 KB
testcase_09 AC 14 ms
4,804 KB
testcase_10 AC 13 ms
4,804 KB
testcase_11 AC 14 ms
4,752 KB
testcase_12 AC 14 ms
4,804 KB
testcase_13 AC 13 ms
4,696 KB
testcase_14 AC 14 ms
4,804 KB
testcase_15 AC 14 ms
4,640 KB
testcase_16 AC 14 ms
4,700 KB
testcase_17 AC 13 ms
4,732 KB
testcase_18 AC 14 ms
4,692 KB
testcase_19 AC 14 ms
4,808 KB
testcase_20 AC 14 ms
4,988 KB
testcase_21 AC 14 ms
4,640 KB
testcase_22 AC 14 ms
4,632 KB
testcase_23 AC 14 ms
4,640 KB
testcase_24 AC 13 ms
4,640 KB
testcase_25 AC 35 ms
4,808 KB
testcase_26 AC 14 ms
4,748 KB
testcase_27 AC 14 ms
4,784 KB
testcase_28 AC 13 ms
4,792 KB
testcase_29 AC 14 ms
4,644 KB
testcase_30 AC 17 ms
4,720 KB
testcase_31 AC 15 ms
4,636 KB
testcase_32 AC 14 ms
4,688 KB
testcase_33 AC 24 ms
5,004 KB
testcase_34 AC 17 ms
4,636 KB
testcase_35 AC 22 ms
4,740 KB
testcase_36 AC 14 ms
4,692 KB
testcase_37 AC 13 ms
4,792 KB
testcase_38 AC 17 ms
4,752 KB
testcase_39 AC 16 ms
4,688 KB
testcase_40 AC 17 ms
4,696 KB
testcase_41 AC 24 ms
4,776 KB
testcase_42 AC 28 ms
4,696 KB
testcase_43 AC 14 ms
4,804 KB
testcase_44 AC 14 ms
4,636 KB
testcase_45 AC 13 ms
4,632 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <algorithm>
#include <cassert>
#include <climits>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <deque>
#include <iomanip>
#include <iostream>
#include <limits>
#include <map>
#include <queue>
#include <set>
#include <sstream>
#include <stack>
#include <string>
#include <vector>

#define FOR(i,k,n) for (int (i)=(k); (i)<(n); ++(i))
#define rep(i,n) FOR(i,0,n)
#define pb push_back
#define eb emplace_back
#define all(v) begin(v), end(v)
#define debug(x) cerr<< #x <<": "<<x<<endl
#define debug2(x,y) cerr<< #x <<": "<< x <<", "<< #y <<": "<< y <<endl

using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> i_i;
typedef pair<i_i, int> p_i;
typedef vector<int> vi;
typedef vector<vector<int> > vvi;
typedef vector<ll> vll;
typedef vector<vector<ll> > vvll;
typedef vector<char> vc;
typedef vector<vector<char> > vvc;
typedef vector<double> vd;
typedef vector<vector<double> > vvd;
template<class T> using vv=vector<vector< T > >;
typedef deque<int> di;
typedef deque<deque<int> > ddi;

// cout vector
template<typename T> ostream& operator<<(ostream& s, const vector<T>& v) {
    int len = v.size();
    for (int i = 0; i < len; ++i) {
        s << v[i]; if (i < len - 1) s << "\t";
    }
    return s;
}

// cout 2-dimentional vector
template<typename T> ostream& operator<<(ostream& s, const vector< vector<T> >& vv) {
    int len = vv.size();
    for (int i = 0; i < len; ++i) {
        s << vv[i] << endl;
    }
    return s;
}

int MAX_PRIME; // in this problem up to 10^6
deque<bool> isprime;
vector<int> primes;
void init_prime(int mp) {
    MAX_PRIME = mp;
    isprime.resize(MAX_PRIME+1, true);
    isprime[0] = isprime[1] = false;
    primes.clear();
    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 min_factor(ll n) {
    rep (i, upper_bound(all(primes), (int)sqrt(n)) - begin(primes)) {
        if ( n % primes[i] == 0 ) {
            return primes[i];
        }
    }
    return (ll)INT_MAX;
}

int main() {
    ll l, h;
    cin >> l >> h;
    init_prime(1000000);
    int maxminfact = 2;
    ll maxminproduct = h - (h % 2 == 1 ? 1 : 0);
    for (int i = upper_bound(all(primes), (int)sqrt(h)) - begin(primes) - 1; i >= 0; --i) {
        if ( primes[i] < maxminfact ) {
            break;
        }
        for (ll j = h / primes[i]; j * primes[i] >= l; --j) {
            int tmp_fact = min(primes[i], (int)min_factor(j));
            if ( tmp_fact > maxminfact ) {
                maxminfact = tmp_fact;
                maxminproduct = j * primes[i];
            }
        }
    }
    printf("%lld\n", maxminproduct);

    return 0;
}
0