結果

問題 No.937 Ultra Sword
コンテスト
ユーザー otamay6
提出日時 2019-11-21 06:42:38
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 164 ms / 3,000 ms
コード長 1,070 bytes
コンパイル時間 648 ms
コンパイル使用メモリ 72,008 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-11-21 01:47:33
合計ジャッジ時間 5,000 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 47
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<vector>
#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);
    REP(i,N){
        for(int j=1;j*j<=A[i];++j){
            if(A[i]%j==0){
                mi[j] += A[i] - A[i]/j;
                if(j*j != A[i]){
                    mi[A[i]/j] += A[i]-j;
                }
            }
        }
    }
    ll X=S;
    rep(i,2,M+1) if(Usable(i,G)){
        X = min(X, S-mi[i]);
    }
    cout<<X<<endl;
}
0