結果

問題 No.937 Ultra Sword
ユーザー leaf_1415
提出日時 2019-11-29 23:34:02
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 1,702 bytes
コンパイル時間 504 ms
コンパイル使用メモリ 61,456 KB
実行使用メモリ 36,396 KB
最終ジャッジ日時 2024-11-21 03:04:09
合計ジャッジ時間 27,394 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 24 WA * 23
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <bitset>
#define llint long long

using namespace std;

llint n;
llint cnt[100005];
bool prime[100005];
vector<llint> pvec;
llint a[100005], b[4000005];

void GaussianElimination(llint a[], int n)
{
	llint r = 0;
	for(int i = 59; i >= 0 && r < n; i--){ //上位ビットから見ていく
		if((a[r]&(1LL<<i)) == 0){
			int p = -1;
			for(int j = r+1; j < n; j++){
				if(a[j] & (1LL<<i)){
					p = j;
					break;
				}
			}
			if(p == -1) goto end;
			swap(a[r], a[p]);
		}
		for(int j = 0; j < n; j++){
			if(j == r) continue;
			if(a[j]&(1LL<<i)) a[j] ^= a[r];
		}
		r++;
		end:;
	}
}

int main(void)
{
	cin >> n;
	for(int i = 0; i < n; i++){
		cin >> a[i];
		cnt[a[i]] += a[i];
	}
	
	for(int i = 2; i < 100005; i++){
		if(prime[i]) continue;
		for(int j = 2*i; j < 100005; j+=i) prime[j] = true;
	}
	
	for(int i = 2; i < 100005; i++){
		if(prime[i]) continue;
		for(int j = 100004/i; j >= 1; j--){
			cnt[j] += cnt[j*i];
		}
	}
	
	for(int i = 0; i < n; i++){
		for(int j = 0; j < 40; j++){
			b[j*n+i] = a[i] << j;
		}
	}
	GaussianElimination(b, n*40);
	
	//for(int i = 0; i < 50; i++) cout << bitset<30>(b[i]) << endl;
	
	llint sum = 0;
	for(int i = 0; i < n; i++) sum += a[i];
	
	llint ans = sum+1;
	for(int i = 1; i <= 100000; i++){
		llint tmp = 0, pre = -1;
		for(int j = 59; j >= 0; j--){
			bool found = false;
			for(int k = pre+1; k < 50; k++){
				if(b[k]&(1LL<<k)){
					found = true;
					tmp ^= b[k];
					pre = k;
					break;
				}
			}
			if(tmp&(1LL<<j) != (i&(1LL<<j)) && !found) goto end;
		}
		//cout << i << " " << sum << " " << cnt[i] << endl;
		ans = min(ans, sum-cnt[i]+cnt[i]/i);
		end:;
	}
	cout << ans << endl;
	
	return 0;
}
0