結果

問題 No.184 たのしい排他的論理和(HARD)
ユーザー koyumeishikoyumeishi
提出日時 2015-04-17 03:41:00
言語 C++11
(gcc 13.3.0)
結果
AC  
実行時間 1,359 ms / 5,000 ms
コード長 2,464 bytes
コンパイル時間 722 ms
コンパイル使用メモリ 81,672 KB
実行使用メモリ 62,208 KB
最終ジャッジ日時 2024-07-04 11:22:35
合計ジャッジ時間 24,049 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,816 KB
testcase_01 AC 1 ms
6,944 KB
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 2 ms
6,940 KB
testcase_04 AC 2 ms
6,940 KB
testcase_05 AC 2 ms
6,944 KB
testcase_06 AC 1 ms
6,940 KB
testcase_07 AC 2 ms
6,944 KB
testcase_08 AC 965 ms
44,928 KB
testcase_09 AC 190 ms
11,776 KB
testcase_10 AC 722 ms
34,944 KB
testcase_11 AC 510 ms
25,856 KB
testcase_12 AC 1,128 ms
51,840 KB
testcase_13 AC 1,254 ms
55,808 KB
testcase_14 AC 678 ms
33,152 KB
testcase_15 AC 1,322 ms
60,160 KB
testcase_16 AC 1,105 ms
50,688 KB
testcase_17 AC 1,177 ms
54,400 KB
testcase_18 AC 2 ms
6,940 KB
testcase_19 AC 2 ms
6,940 KB
testcase_20 AC 157 ms
62,080 KB
testcase_21 AC 1,359 ms
62,016 KB
testcase_22 AC 1,348 ms
62,208 KB
testcase_23 AC 2 ms
6,948 KB
testcase_24 AC 2 ms
6,940 KB
testcase_25 AC 2 ms
6,940 KB
testcase_26 AC 2 ms
6,944 KB
testcase_27 AC 2 ms
6,944 KB
testcase_28 AC 796 ms
38,656 KB
testcase_29 AC 1,105 ms
53,504 KB
testcase_30 AC 954 ms
47,488 KB
testcase_31 AC 821 ms
41,472 KB
testcase_32 AC 1,024 ms
51,072 KB
testcase_33 AC 1,358 ms
60,548 KB
testcase_34 AC 1,346 ms
59,776 KB
testcase_35 AC 1,326 ms
61,564 KB
testcase_36 AC 1,347 ms
61,440 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <cstdio>
#include <sstream>
#include <map>
#include <string>
#include <algorithm>
#include <queue>
#include <cmath>
#include <set>
using namespace std;

long long extgcd(long long a, long long b, long long &x, long long &y){
	long long d=a;
	if(b!=0){
		d = extgcd(b, a%b, y, x);
		y -= (a/b) * x;
	}else{
		x = 1;
		y = 0;
	}
	return d;
}

long long mod_inverse(long long a, long long m){
	long long x,y;
	extgcd(a,m,x,y);
	return (m+x%m)%m;
}


// A[n*p] * X[p*m] = B[n*m]
template<class T = int>
int gaussian_elimination_with_mod(vector<vector<T>>& A, vector<vector<T>>& B, int n, int p, int m, const T mod){
	int R = max(n,p); int C = p+m;
	vector<vector<T>> V( R, vector<T>(C, 0) );
	for(int i=0; i<n; i++)
		for(int j=0; j<p; j++) V[i][j] = A[i][j];
	for(int i=0; i<n; i++)
		for(int j=0; j<m; j++) V[i][j+p] = B[i][j];


	int rank = 0;
	int row = 0;

	vector<int> left(R, -1);

	//foward
	for(int col=0; col<C && row<R; col++){
		//pivot
		T val = abs( V[row][col] );
		int pivot = row;
		for(int j=row+1; j<R; j++){
			if(val < abs( V[j][col] )){
				val = abs( V[j][col] );
				pivot = j;
			}
		}
		if(pivot != row) swap(V[row], V[pivot]);

		if(val == 0) continue;

		T inv = mod_inverse(val, mod);
		for(int j=row+1; j<R; j++){
			T c = (V[j][col] * inv + mod) % mod;
			for(int k=col; k<p+m; k++) V[j][k] = ((V[j][k] - V[row][k] * c) % mod + mod) % mod;
		}

		left[row] = col;
		row++;
		rank++;
	}

	return rank;

	//backward
	for(int i=R-1; i>=0; i--){
		bool zero = true;
		bool valid = true;
		for(int col=0; col<p; col++) if(V[i][col] != 0) zero = false;
		for(int col = p; zero && col<C; col++) if(V[i][col] != 0) valid = false;

		if(valid == false) return -1;	//no solution

		if(left[i] == -1) continue;

		T inv = mod_inverse(V[i][ left[i] ], mod);
		for(int j=left[i]; j<C; j++) V[i][j] = (V[i][j] * inv) % mod;

		for(int j=i-1; j>=0; j--){
			for(int k=0; k<C; k++){
				V[j][k] = ( (V[j][k] - V[j][i] * V[i][k]) % mod + mod ) % mod;
			}
		}
	}


	//return V;
	return rank;

}

int main(){
	int n;
	cin >> n;
	vector<long long> a(n);
	for(int i=0; i<n; i++){
		cin >> a[i];
	}

	vector<vector<int>> A(n, vector<int>(64, 0));
	for(int i=0; i<n; i++){
		for(long long j=0; j<64; j++){
			A[i][j] = (a[i] >> j) & 1LL;
		}
	}

	vector<vector<int>> B(64);

	long long rank = gaussian_elimination_with_mod<int>(A,B, n,64,0, 2);
	cerr << rank << endl;
	cout << (1LL<<rank) << endl;

	return 0;
}
0