結果
問題 | No.641 Team Contest Estimation |
ユーザー |
|
提出日時 | 2023-12-17 14:46:04 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 160 ms / 2,000 ms |
コード長 | 772 bytes |
コンパイル時間 | 358 ms |
コンパイル使用メモリ | 82,420 KB |
実行使用メモリ | 90,404 KB |
最終ジャッジ日時 | 2024-09-27 08:01:08 |
合計ジャッジ時間 | 2,244 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 9 |
ソースコード
from collections import *from itertools import *from functools import *from heapq import *import sys,mathinput = sys.stdin.readlineN,K = map(int,input().split())A = list(map(int,input().split()))mod = 10**9 + 9B = [0]*KS = 0I = pow(2,K-1,mod)for i in range(K):S = (S + N*(1<<i)*I)%modb1 = 0for a in A:if (a>>i)&1:b1 += 1B[i] = b1SS = 0for i in range(K):SS = (SS + I*(1<<(2*i))*(B[i]**2))%modSS = (SS + I*(1<<(2*i))*((N-B[i])**2))%mod# print(SS)for i in range(K-1):bi = B[i]Ii = (1<<i)for j in range(i+1,K):bj = B[j]Ij = (1<<j)SS = (SS + I*(Ii*Ij*N*N))%mod# SS = (SS + I*(Ii*Ij*(N-bi)*(N-bj)))%modprint(S)print((SS*(1<<K) - S**2)%mod)