#302953 (C++14) No.719 Coprime

提出ソース

結果

問題	No.719 Coprime
ユーザー	jell
提出日時	2018-12-09 21:59:58
言語	C++14 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	139 ms / 3,000 ms
コード長	6,207 bytes
コンパイル時間	1,464 ms
コンパイル使用メモリ	121,728 KB
実行使用メモリ	5,888 KB
最終ジャッジ日時	2024-09-16 15:37:43
合計ジャッジ時間	3,430 ms
ジャッジサーバーID （参考情報）	judge1 / judge4

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 5,248 KB
testcase_01	AC	2 ms 5,376 KB
testcase_02	AC	1 ms 5,376 KB
testcase_03	AC	1 ms 5,376 KB
testcase_04	AC	1 ms 5,376 KB
testcase_05	AC	1 ms 5,376 KB
testcase_06	AC	1 ms 5,376 KB
testcase_07	AC	1 ms 5,376 KB
testcase_08	AC	1 ms 5,376 KB
testcase_09	AC	2 ms 5,376 KB
testcase_10	AC	2 ms 5,376 KB
testcase_11	AC	2 ms 5,376 KB
testcase_12	AC	1 ms 5,376 KB
testcase_13	AC	2 ms 5,376 KB
testcase_14	AC	2 ms 5,376 KB
testcase_15	AC	2 ms 5,376 KB
testcase_16	AC	2 ms 5,376 KB
testcase_17	AC	2 ms 5,376 KB
testcase_18	AC	2 ms 5,376 KB
testcase_19	AC	2 ms 5,376 KB
testcase_20	AC	2 ms 5,376 KB
testcase_21	AC	2 ms 5,376 KB
testcase_22	AC	2 ms 5,376 KB
testcase_23	AC	2 ms 5,376 KB
testcase_24	AC	2 ms 5,376 KB
testcase_25	AC	2 ms 5,376 KB
testcase_26	AC	2 ms 5,376 KB
testcase_27	AC	2 ms 5,376 KB
testcase_28	AC	2 ms 5,376 KB
testcase_29	AC	2 ms 5,376 KB
testcase_30	AC	2 ms 5,376 KB
testcase_31	AC	2 ms 5,376 KB
testcase_32	AC	2 ms 5,376 KB
testcase_33	AC	2 ms 5,376 KB
testcase_34	AC	2 ms 5,376 KB
testcase_35	AC	2 ms 5,376 KB
testcase_36	AC	2 ms 5,376 KB
testcase_37	AC	2 ms 5,376 KB
testcase_38	AC	1 ms 5,376 KB
testcase_39	AC	2 ms 5,376 KB
testcase_40	AC	2 ms 5,376 KB
testcase_41	AC	1 ms 5,376 KB
testcase_42	AC	2 ms 5,376 KB
testcase_43	AC	2 ms 5,376 KB
testcase_44	AC	2 ms 5,376 KB
testcase_45	AC	2 ms 5,376 KB
testcase_46	AC	2 ms 5,376 KB
testcase_47	AC	2 ms 5,376 KB
testcase_48	AC	2 ms 5,376 KB
testcase_49	AC	2 ms 5,376 KB
testcase_50	AC	2 ms 5,376 KB
testcase_51	AC	3 ms 5,376 KB
testcase_52	AC	4 ms 5,376 KB
testcase_53	AC	9 ms 5,376 KB
testcase_54	AC	21 ms 5,376 KB
testcase_55	AC	43 ms 5,376 KB
testcase_56	AC	100 ms 5,504 KB
testcase_57	AC	104 ms 5,376 KB
testcase_58	AC	131 ms 5,888 KB
testcase_59	AC	139 ms 5,888 KB
testcase_60	AC	131 ms 5,888 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

// #include <boost/functional/hash.hpp>
// #include <boost/multiprecision/cpp_int.hpp>
// using mlint = boost::multiprecision::cpp_int;

#include <cassert>
#include <deque>
#include <queue>
#include <stack>
#include <algorithm>
#include <iostream>
#include <set>
#include <map>
#include <list>
#include <cmath>
#include <complex>
#include <functional>
#include <random>
#include <unordered_map>
#include <unordered_set>
#include <bitset>
#include <iomanip>

using namespace std;
using i64 = int_fast64_t;
using db = double;
using pii = pair<int,int>;
using pli = pair<int_fast64_t,int>;
using pdi = pair<db,int>;
using plpii = pair<int_fast64_t,pair<int,int>>;
template <class T> using mat = vector<vector<T>>;
template <class T> using rprque = priority_queue<T,vector<T>,greater<T>>;

#define esc(x) cout << (x) << endl,exit(0)
#define inf(T) (numeric_limits<T>::max() / 2 - 1)
#define each(i,v) for(auto i = begin(v); i != end(v); ++i)
#define reach(i,v) for(auto i = rbegin(v); i != rend(v); ++i)
#define urep(i,s,t) for(int_fast64_t i = (s); i <= (t); ++i)
#define drep(i,s,t) for(int_fast64_t i = (s); i >= (t); --i)
#define rep(i,n) urep(i,0,(n) - 1)
#define rep1(i,n) urep(i,1,(n))
#define rrep(i,n) drep(i,(n) - 1,0)
#define all(v) begin(v),end(v)
#define rall(v) rbegin(v),rend(v)
#define vct vector
#define u_map unordered_map
#define u_set unordered_set
#define l_bnd lower_bound
#define u_bnd upper_bound
#define rsz resize
#define ers erase
#define emp emplace
#define emf emplace_front
#define emb emplace_back
#define pof pop_front
#define pob pop_back
#define mkp make_pair
#define mkt make_tuple
#define fir first
#define sec second

struct setupper {
    setupper(uint_fast8_t prec) {
        cin.tie(0);
        ios::sync_with_stdio(false);
        cout << fixed << setprecision(prec);
        #ifdef Local
        #define debug 1
        #define print(x) cout << #x << " = " << (x) << " (L" << __LINE__ << ")" << endl
        cout << "-- Execution At Local ---\n" << "\n<Standard Output>\n" << "\n----------------\n\n";
        auto print_atexit = []() {
            time_t end_time = time(NULL);
            struct tm *ret = localtime(&end_time);
            cout << "\n----------------\n";
            cout << "\nSuccessfully Executed At " << (ret -> tm_hour) << ":" << (ret -> tm_min) << ":" << (ret -> tm_sec) << "\n\n";
        };
        atexit(print_atexit);
        #endif
    }
} setupper_obj(10);

template <class T, class U> ostream& operator << (ostream &s, pair<T,U> p) { return s << p.fir << " " << p.sec; }
template <class T> ostream& operator << (ostream &s, vector<T> &v) { for(auto i = v.begin(); i != v.end(); ++i) { if(i != begin(v)) s << " "; s << *i; } return s; }

template <class T> void hash_combine(size_t &seed, T const &key) {
    hash<T> hasher;
    seed ^= hasher(key) + 0x9e3779b9 + (seed << 6) + (seed >> 2);
}
namespace std {
    template <class T, class U> struct hash<pair<T,U>> {
        size_t operator()(pair<T,U> const &p) const {
            size_t seed(0);
            hash_combine(seed,p.first);
            hash_combine(seed,p.second);
            return seed;
        }
    };
    template <class T, class U, class V> struct hash<tuple<T,U,V>> {
        size_t operator()(tuple<T,U,V> const &t) const {
            size_t seed(0);
            hash_combine(seed,get<0>(t));
            hash_combine(seed,get<1>(t));
            hash_combine(seed,get<2>(t));
            return seed;
        }
    };
}

const auto add = [](auto &x, auto y) { x += y; };
const auto mul = [](auto &x, auto y) { x *= y; };
const auto lam_min = [](auto x, auto y) { return min(x,y); };
const auto lam_max = [](auto x, auto y) { return max(x,y); };
const auto chmax = [](auto &m, auto x) { if(m < x){ m = x; return true; } return false; };
const auto chmin = [](auto &m, auto x) { if(m > x){ m = x; return true; } return false; };
bool odd(i64 x) { return x & 1; }
bool even(i64 x) { return !odd(x); }
bool parity(i64 x, i64 y) { return odd(x) ^ even(y); }
bool bit(i64 n, uint8_t e) { return n >> e & 1; }
i64 mask(i64 n, uint8_t e) { return n & (1 << e) - 1; }
int ilog(i64 x, uint64_t b = 2) { if(x) return ilog(x / b, b) + 1; return -1; }
i64 qceil(i64 x, i64 y) { return x > 0 ? (x - 1) / y + 1 :  x / y; }
i64 gcd(i64 a, i64 b) { if(!b) return abs(a); return gcd(b, a % b); }
i64 lcm(i64 a, i64 b) { if(a | b) return abs(a / gcd(a, b) * b); return 0; }
i64 extgcd(i64 a, i64 b, i64 &x, i64 &y) {
    i64 d = a;
    if(b) d = extgcd(b,a % b,y,x),y -= (a / b) * x;
    else x = 1,y = 0;
    return d;
}
template <class F> i64 binsr(i64 ok, i64 ng, const F &fn) { while(abs(ok - ng) > 1){ i64 mid = ok + ng >> 1; (fn(mid) ? ok : ng) = mid; } return ok; }
template <class T> T cmprs(T &v) {
    T tmp,ret = v; sort(all(tmp = v));
    tmp.erase(unique(all(tmp)),end(tmp));
    each(i,ret) *i = l_bnd(all(tmp),*i) - begin(tmp) + 1;
    return ret;
}

constexpr int dx[9] = {1,0,-1,0,1,-1,-1,1,0};
constexpr int dy[9] = {0,1,0,-1,1,1,-1,-1,0};
constexpr i64 mod = 1e9 + 7;
constexpr db gold = 1.618033988;
constexpr db eps = 1e-15;

bool isp[1300];
vct<int> smp;
vct<vct<int>> grp; 
int n;
int dp[220][1<<20];

int main() {
    cin>>n;
    fill(isp+2,isp+1300,true);
    for(int d=2;d*d<=n;d++) {
        if(isp[d]) {
            smp.emb(d);
            int k=d;
            while((k+=d)<=n) {
                isp[k]=false;
            }
        }
    }
    bool in[1300]={};
    int idx=0;
    drep(i,n,2) {
        if(isp[i]) {
            vct<int> v;
            grp.emb(v);
            urep(j,2,n) {
                if(!in[j]&&j%i==0) in[j]=true,grp.back().emb(j);
            }
        }
    }
    int sz=smp.size();
    reach(i,grp) {
        auto nowv=*i;
        int *crt=dp[idx];
        int *nxt=dp[++idx];
        rep(s,1<<sz) {
            chmax(nxt[s],crt[s]);
            for(int k:nowv) {
                bool psbl=1;
                rep(j,sz) if(bit(s,j)&&k%smp[j]==0) { psbl=0; break; }
                if(psbl) {
                    int msk=0;
                    rep(j,sz) if(k%smp[j]==0) msk|=1<<j;
                    chmax(nxt[s|msk],crt[s]+k);
                } 
            }
        }
        if(idx==grp.size()) esc(nxt[(1<<sz)-1]);
    }
}

yukicoder

結果

テストケース

ソースコード