結果

問題 No.184 たのしい排他的論理和(HARD)
ユーザー None
提出日時 2021-04-04 16:03:25
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 253 ms / 5,000 ms
コード長 3,400 bytes
コンパイル時間 552 ms
コンパイル使用メモリ 82,432 KB
実行使用メモリ 92,928 KB
最終ジャッジ日時 2024-12-27 20:18:40
合計ジャッジ時間 8,248 ms
ジャッジサーバーID
(参考情報)
judge4 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 34
権限があれば一括ダウンロードができます

ソースコード

diff #

class GE():
    """
    mod2の場合、
    A=[int,int,int,...]
    M=max(A[i].bit_length() for i in range(len(A)))
    とする
    """
    def __init__(self,A,mod,M=0):
        self.N=len(A)
        if M:
            self.M=M
        else:
            self.M=len(A[0])
        self.mod=mod
        self.A=A[:]
        self.uptri=None
        self.A2=None
        self.pivot=[]
        self.R=None

    def ut(self):
        if self.mod!=2:
            if self.uptri is not None:
                return self.uptri
            self.uptri=self.A[:]
            c=0
            for i in range(self.N):
                if i+c>=self.M:
                    break
                while self.uptri[i][i+c]==0:
                    for j in range(i+1,self.N):
                        if self.uptri[j][i+c]:
                            self.uptri[i],self.uptri[j]=self.uptri[j],self.uptri[i]
                            break
                    else:
                        c+=1
                        if i+c==self.M:
                            break
                else:
                    self.pivot.append((i,i+c))
                    t=pow(self.uptri[i][i+c],self.mod-2,self.mod)
                    for j in range(i+1,self.N):
                        tj=t*self.uptri[j][i+c]
                        self.uptri[j][i+c:]=[(aj-tj*ai)%self.mod for ai,aj in
                                             zip(self.uptri[i][i+c:],self.uptri[j][i+c:])]
                    continue
                break
            for pi,pj in self.pivot[::-1]:
                t=pow(self.uptri[pi][pj],self.mod-2,self.mod)
                self.uptri[pi][pj:]=[(ai*t)%self.mod for ai in self.uptri[pi][pj:]]
                for i in range(pi-1,-1,-1):
                    ti=self.uptri[i][pj]
                    self.uptri[i][pj:]=[(ai-api*ti)%self.mod for ai,api in zip(self.uptri[i][pj:],self.uptri[pi][pj:])]

            return self.uptri
        else:
            if self.A2 is not None:
                return self.A2
            self.A2=self.A[:]
            c=0
            for i in range(self.N):
                if i+c>=self.M:
                    break
                while not self.A2[i]&2**(i+c):
                    for j in range(i+1,self.N):
                        if self.A2[j]&2**(i+c):
                            self.A2[i],self.A2[j]=self.A2[j],self.A2[i]
                            break
                    else:
                        c+=1
                        if i+c==self.M:
                            break
                else:
                    self.pivot.append((i,i+c))
                    for j in range(i+1,self.N):
                        if self.A2[j]&2**(i+c):
                            self.A2[j]^=self.A2[i]
                    continue
                break
            for pi,pj in self.pivot[::-1]:
                for i in range(pi-1,-1,-1):
                    if 2**pj&self.A2[i]:
                        self.A2[i]^=self.A2[pi]

            return self.A2

    def rank(self):
        if self.R is not None:
            return self.R
        self.ut()
        self.R=len(self.pivot)
        return self.R

    def solve(self):
        self.ut()

############################################################################
import sys
input = sys.stdin.readline


N=int(input())
A=list(map(int, input().split()))
mat=GE(A,2,M=62)
rank=mat.rank()
print(pow(2,rank))
0