結果

問題 No.937 Ultra Sword
ユーザー otamay6
提出日時 2019-11-25 17:48:22
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 133 ms / 3,000 ms
コード長 1,053 bytes
コンパイル時間 2,188 ms
コンパイル使用メモリ 174,476 KB
実行使用メモリ 14,848 KB
最終ジャッジ日時 2024-11-21 02:10:06
合計ジャッジ時間 6,611 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
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;
    int M=0, G=0;
    ll S=0;
    vector<int> A(N);
    REP(i,N){
        cin>>A[i];
        M=max(M, A[i]);
        G=GCD(G, A[i]);
        S += A[i];
    }
    vector<ll> mi(M+1, 0);
    vector<vector<int>> div(M+1);
    for(int i=1;i<=M;++i){
        for(int j=i;j<=M;j+=i){
            div[j].push_back(i);
        }
    }
    REP(i,N){
        for(auto y:div[A[i]]){
            mi[y]+=A[i]-A[i]/y;
        }
    }
    ll X=S;
    rep(i,2,M+1) if(Usable(i,G)){
        X = min(X, S-mi[i]);
    }
    cout<<X<<endl;
}
0