結果

問題 No.719 Coprime
ユーザー はまやんはまやん
提出日時 2018-08-05 17:26:15
言語 C++14
(gcc 13.3.0 + boost 1.87.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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 61
権限があれば一括ダウンロードができます

ソースコード

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