結果

問題 No.937 Ultra Sword
ユーザー mfbgjsczmfbgjscz
提出日時 2019-11-25 19:48:20
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
RE  
実行時間 -
コード長 2,159 bytes
コンパイル時間 2,003 ms
コンパイル使用メモリ 180,768 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-05-01 02:25:25
合計ジャッジ時間 8,095 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
testcase_01 RE -
testcase_02 RE -
testcase_03 RE -
testcase_04 RE -
testcase_05 RE -
testcase_06 RE -
testcase_07 RE -
testcase_08 RE -
testcase_09 RE -
testcase_10 RE -
testcase_11 RE -
testcase_12 RE -
testcase_13 RE -
testcase_14 RE -
testcase_15 RE -
testcase_16 RE -
testcase_17 RE -
testcase_18 RE -
testcase_19 RE -
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 RE -
testcase_28 RE -
testcase_29 RE -
testcase_30 RE -
testcase_31 RE -
testcase_32 RE -
testcase_33 RE -
testcase_34 RE -
testcase_35 RE -
testcase_36 RE -
testcase_37 RE -
testcase_38 RE -
testcase_39 RE -
testcase_40 RE -
testcase_41 RE -
testcase_42 RE -
testcase_43 RE -
testcase_44 RE -
testcase_45 RE -
testcase_46 RE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define lli long long int
#define uli unsigned long long int
#define INF 999999999999999999
#define rep(i,m,n) for(lli i = m;i < n;i++)
#define rrep(i,m,n) for(lli i=m-1;i>=n;i--)
#define pb(n) push_back(n)
#define UE(N) N.erase(unique(N.begin(),N.end()),N.end());
#define Sort(n) sort(n.begin(), n.end())
#define Rev(n) reverse(n.begin(),n.end())
#define Out(S) cout << S << endl
#define NeOut(S) cout << S
#define HpOut(S) cout << setprecision(30) << S << endl
#define Vec(K,L,N,S) vector<L> K(N,S)
#define DV(K,L,N,M,S) vector<vector<L>> K(N,vector<L>(M,S))
#define TV(K,L,N,M,R,S) vector<vector<vector<L>>> K(N,vector<vector<L>>(M,vector<L>(R,S)))
#define pint pair<lli,lli>
#define paf(L,R) pair<L,R>
#define mod 1000000007
#define MAX 5100000
#define chmax(a, b) a = (((a)<(b)) ? (b) : (a))
#define chmin(a, b) a = (((a)>(b)) ? (b) : (a))
map<lli,lli> divisor(vector<lli>HP){
  map<lli,lli>res;
  for(auto N:HP){
    for(lli i=1;i*i<=N;i++){
      if(N%i==0){
        res[i]+=N-N/i;
        if(i!=N/i)res[N/i]+=N-N/(N/i);
      }
    }
  }
  return res;
}
lli f2_gcd(bitset<20> A,bitset<20> B) {
	int i;
	for(i=20;i>=0;i--) {
		if(A[i]==1 && B[i]==1) {
			B^=A;
			break;
		}
		if(A[i]==1) break;
		if(B[i]==1) {
			swap(A,B);
			break;
		}
	}
	
	// A>B
	while(1){
		int x;
		for(x=i-1;x>=0;x--)if(B[x]) break;
		if(x<0) break;
		while(i>=x) {
			if(A[i]) A^=B<<(i-x);
			i--;
		}
		i=x;
		swap(A,B);
	}
  lli x=0;
  rep(i,0,20)if(A[i]==1)x+=pow(2,i);
  return x;
}
bool is_multiple(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(){
  lli A,B,C,D,E,F,N,M,K,L,R,X,Y,H,W,sum=0,num=0,flag=0;string S,T;
  assert(cin >> N);
  assert(N>100000||2>N);
  Y=0;
  Vec(HP,lli,N,0);
  rep(i,0,N){
      assert(cin >> HP[i]);
      assert(HP[i]>400000);
  }
  X=f2_gcd(HP[0],HP[1]);
  rep(i,2,N)X=f2_gcd(HP[i],X);
  rep(i,0,N){
    chmax(Y,HP[i]);
    sum+=HP[i];
  }
  num=sum;
  auto St=divisor(HP);
  for(auto v:St)if(is_multiple(v.first,X))num=min(num,sum-v.second);
  Out(num);
}
0