結果
問題 | No.16 累乗の加算 |
ユーザー |
|
提出日時 | 2021-02-09 22:57:29 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 43 ms / 5,000 ms |
コード長 | 3,834 bytes |
コンパイル時間 | 163 ms |
コンパイル使用メモリ | 82,340 KB |
実行使用メモリ | 64,008 KB |
最終ジャッジ日時 | 2024-07-07 04:44:39 |
合計ジャッジ時間 | 1,181 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 14 |
ソースコード
import mathimport syssys.setrecursionlimit(10**7)#競技プログラミング対整数問題のライブラリーですclass integerlib():def __init__(self):passdef primeset(self,N): #N以下の素数をsetで求める.エラトステネスの篩O(√Nlog(N))lsx = [1]*(N+1)for i in range(2,int(-(-N**0.5//1))+1):if lsx[i] == 1:for j in range(i,N//i+1):lsx[j*i] = 0setprime = set()for i in range(2,N+1):if lsx[i] == 1:setprime.add(i)return setprimedef defprime(self,N):#素数かどうかの判定、エラトステネスの篩O(√Nlog(N))return N in self.primeset(N)def gcd(self,ls):#最大公約数ls = list(ls)ans = 0for i in ls:ans = math.gcd(ans,i)return ansdef lmc(self,ls):#最小公倍数ls = list(ls)ans = self.gcd(ls)for i in ls:ans = self.lmcsub(ans,i)return ansdef lmcsub(self,a,b):gcd = math.gcd(a,b)lmc = (a*b)//gcdreturn lmcdef factorization(self,N):#素因数分解√Narr = []temp = Nfor i in range(2, int(-(-N**0.5//1))+1):if temp%i==0:cnt=0while temp%i==0:cnt+=1temp //= iarr.append([i, cnt])if temp!=1:arr.append([temp, 1])if arr==[]:arr.append([N, 1])return arr #[素因数、個数]def factorizationset(self,N):#素因数分解√N,含まれている素因数の種類if N == 1:return set()ls = self.factorization(N)setf = set()for j in ls:setf.add(j[0])return setfdef divisorsnum(self,N):#約数の個数ls = []for i in self.factorization(N):ls.append(i[1])d = 1for i in ls:d *= i+1return ddef Eulerfunc(self,N):#オイラー関数正の整数Nが与えられる。1,2,…,Nのうち、Nと互いに素であるものの個数を求めよ。ls = list(self.factorizationset(N))ls2 = [N]for i in ls:ls2.append(ls2[-1]-ls2[-1]//i)return ls2[-1]def make_divisors(self,N):#約数列挙O(√N)lower_divisors , upper_divisors = [], []i = 1while i*i <= N:if N % i == 0:lower_divisors.append(i)if i != N // i:upper_divisors.append(N//i)i += 1return lower_divisors + upper_divisors[::-1]def invmod(self,a,mod):#mod逆元if a == 0:return 0if a == 1:return 1return (-self.invmod(mod % a, mod) * (mod // a)) % moddef cmbmod(self,n, r, mod):#nCr % modinv = [0,1]for i in range(2, n + 1):inv.append((-inv[mod % i] * (mod // i)) % mod)cmd = 1for i in range(1,min(r,n-r)+1):cmd = (cmd*(n-i+1)*inv[i])%modreturn cmddef permmod(self,n, r, mod):#nPr % modperm = 1for i in range(n,r-1,-1):perm = (perm*i)%modreturn permdef modPow(self,a,n,mod):#繰り返し二乗法 a**n % modif n==0:return 1if n==1:return a%modif n % 2 == 1:return (a*self.modPow(a,n-1,mod)) % modt = self.modPow(a,n//2,mod)return (t*t)%modx,N = map(int,input().split())lsN = list(map(int,input().split()))mod = 1000003IT = integerlib()ans = 0for i in range(N):ans += IT.modPow(x,lsN[i],mod)ans %= modprint(ans)