結果

問題 No.719 Coprime
ユーザー はまやんはまやんはまやんはまやん
提出日時 2018-08-05 17:26:15
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 66 ms / 3,000 ms
コード長 2,640 bytes
コンパイル時間 1,724 ms
コンパイル使用メモリ 173,760 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-19 18:06:40
合計ジャッジ時間 3,661 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,812 KB
testcase_01 AC 1 ms
6,944 KB
testcase_02 AC 1 ms
6,940 KB
testcase_03 AC 2 ms
6,940 KB
testcase_04 AC 2 ms
6,940 KB
testcase_05 AC 2 ms
6,944 KB
testcase_06 AC 2 ms
6,940 KB
testcase_07 AC 1 ms
6,940 KB
testcase_08 AC 2 ms
6,940 KB
testcase_09 AC 2 ms
6,940 KB
testcase_10 AC 2 ms
6,940 KB
testcase_11 AC 2 ms
6,940 KB
testcase_12 AC 2 ms
6,944 KB
testcase_13 AC 2 ms
6,940 KB
testcase_14 AC 1 ms
6,944 KB
testcase_15 AC 2 ms
6,940 KB
testcase_16 AC 2 ms
6,940 KB
testcase_17 AC 1 ms
6,940 KB
testcase_18 AC 2 ms
6,940 KB
testcase_19 AC 1 ms
6,944 KB
testcase_20 AC 2 ms
6,940 KB
testcase_21 AC 1 ms
6,944 KB
testcase_22 AC 1 ms
6,944 KB
testcase_23 AC 2 ms
6,944 KB
testcase_24 AC 1 ms
6,940 KB
testcase_25 AC 2 ms
6,944 KB
testcase_26 AC 1 ms
6,944 KB
testcase_27 AC 1 ms
6,940 KB
testcase_28 AC 2 ms
6,940 KB
testcase_29 AC 2 ms
6,944 KB
testcase_30 AC 1 ms
6,940 KB
testcase_31 AC 2 ms
6,944 KB
testcase_32 AC 1 ms
6,940 KB
testcase_33 AC 1 ms
6,944 KB
testcase_34 AC 2 ms
6,944 KB
testcase_35 AC 2 ms
6,940 KB
testcase_36 AC 2 ms
6,940 KB
testcase_37 AC 2 ms
6,940 KB
testcase_38 AC 2 ms
6,940 KB
testcase_39 AC 2 ms
6,944 KB
testcase_40 AC 1 ms
6,944 KB
testcase_41 AC 2 ms
6,940 KB
testcase_42 AC 2 ms
6,944 KB
testcase_43 AC 2 ms
6,940 KB
testcase_44 AC 1 ms
6,944 KB
testcase_45 AC 2 ms
6,940 KB
testcase_46 AC 2 ms
6,944 KB
testcase_47 AC 2 ms
6,944 KB
testcase_48 AC 2 ms
6,944 KB
testcase_49 AC 2 ms
6,940 KB
testcase_50 AC 2 ms
6,940 KB
testcase_51 AC 2 ms
6,944 KB
testcase_52 AC 4 ms
6,940 KB
testcase_53 AC 6 ms
6,940 KB
testcase_54 AC 20 ms
6,944 KB
testcase_55 AC 21 ms
6,940 KB
testcase_56 AC 48 ms
6,940 KB
testcase_57 AC 49 ms
6,940 KB
testcase_58 AC 65 ms
6,940 KB
testcase_59 AC 66 ms
6,940 KB
testcase_60 AC 66 ms
6,944 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<b;i++)
#define rrep(i,a,b) for(int i=a;i>=b;i--)
#define fore(i,a) for(auto &i:a)
#define all(x) (x).begin(),(x).end()
#pragma GCC optimize ("-O3")
using namespace std; void _main(); int main() { cin.tie(0); ios::sync_with_stdio(false); _main(); }
typedef long long ll; const int inf = INT_MAX / 2; const ll infl = 1LL << 60;
template<class T>bool chmax(T &a, const T &b) { if (a<b) { a = b; return 1; } return 0; }
template<class T>bool chmin(T &a, const T &b) { if (b<a) { a = b; return 1; } return 0; }
//---------------------------------------------------------------------------------------------------
vector<int> makePrimes(int n) { // [2,n]
    vector<int> res, primes(n + 1, 1);
    primes[0] = primes[1] = 0;
    rep(i, 2, sqrt(n)) if (primes[i]) for (int j = 0; i * (j + 2) <= n; j++) primes[i * (j + 2)] = 0;
    rep(i, 2, n + 1) if (primes[i]) res.push_back(i);
    return res;
}
/*---------------------------------------------------------------------------------------------------
　　　　　　　　　　　 ∧＿∧  
　　　　　 ∧＿∧ 　（´<_｀ ）　 Welcome to My Coding Space!
　　　　 （ ´_ゝ`）　/　 ⌒i     
　　　　／　　　＼　 　  |　|     
　　　 /　　 /￣￣￣￣/　　|  
　 ＿_(__ﾆつ/　    ＿/ .| .|＿＿＿＿  
　 　　　＼/＿＿＿＿/　（u　⊃  
---------------------------------------------------------------------------------------------------*/






int N;
int dp[206][1<<11];
int f[1 << 11];
//---------------------------------------------------------------------------------------------------
void _main() {
    cin >> N;

    auto allpr = makePrimes(N);
    auto halfpr = makePrimes(sqrt(N) + 1);

    int n = allpr.size();
    int m = halfpr.size();

    rep(i, m, n) rep(msk, 0, 1 << m) {
        for (int j = 1; j * allpr[i] <= N; j++) {
            int num = j * allpr[i];
            int ms = 0;
            rep(k, 0, m) if (num % halfpr[k] == 0) ms += 1 << k;
            if ((msk & ms) == 0) chmax(dp[i + 1][msk + ms], dp[i][msk] + num);
        }
    }

    rep(i, 2, N + 1) {
        int ms = 0;
        int ok = 1;
        rep(j, 0, n) {
            if (i % allpr[j] == 0) {
                if (j < m) ms += 1 << j;
                else ok = 0;
            }
        }
        if (ok) chmax(f[ms], i);
    }

    rep(msk, 0, 1 << m) for (int msk2 = msk; msk2 >= 0; msk2--) {
        msk2 &= msk;
        chmax(f[msk], f[msk2] + f[msk - msk2]);
    }

    int ans = 0;
    rep(msk, 0, 1 << m) chmax(ans, dp[n][msk] + f[(1 << m) - 1 - msk]);
    cout << ans << endl;
}
0