結果

問題 No.937 Ultra Sword
ユーザー otamay6otamay6
提出日時 2019-11-30 02:45:07
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 41 ms / 3,000 ms
コード長 1,043 bytes
コンパイル時間 1,819 ms
コンパイル使用メモリ 172,056 KB
実行使用メモリ 6,824 KB
最終ジャッジ日時 2024-11-21 03:20:13
合計ジャッジ時間 3,928 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 47
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include"bits/stdc++.h"
#define REP(i,n) for(int i=0,i##_len=(n);i<i##_len;++i)
#define rep(i,a,b) for(int i=int(a);i<int(b);++i)
using namespace std;
using ll = long long;

int GCD(int x,int y){
    if(x>y) swap(x,y);
    int Lx=32-__builtin_clz(x),Ly=32-__builtin_clz(y);
    return x?GCD(x,y^(x<<(Ly-Lx))):y;
}

bool Usable(int x,int g){
    int Lx=32-__builtin_clz(x),Lg=32-__builtin_clz(g);
    while(Lx>=Lg&&x!=0){
        x^=g<<(Lx-Lg);
        Lx=32-__builtin_clz(x);
    }
    return x==0;
}

int main(){
    int N;cin>>N;
    assert(2<=N&&N<=100000);
    int M=0, G=0;
    ll S=0;
    vector<int> A(N);
    vector<ll> f(100001);
    REP(i,N){
        cin>>A[i];
        assert(1<=A[i]&&A[i]<=1000000);
        M=max(M, A[i]);
        G=GCD(G, A[i]);
        S += A[i];
        f[A[i]]++;
    }
    vector<ll> mi(M+1, 0);
    for(int i=1;i<=M;++i){
        for(int j=i;j<=M;j+=i){
            mi[i]+=f[j]*(j-j/i);
        }
    }
    ll X=S;
    rep(i,2,M+1) if(Usable(i,G)){
        X = min(X, S-mi[i]);
    }
    cout<<X<<endl;
}
0