結果
| 問題 |
No.1084 積の積
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2022-08-02 01:38:08 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 1,017 ms / 2,000 ms |
| コード長 | 3,504 bytes |
| コンパイル時間 | 1,054 ms |
| コンパイル使用メモリ | 82,176 KB |
| 実行使用メモリ | 98,048 KB |
| 最終ジャッジ日時 | 2024-07-22 17:12:13 |
| 合計ジャッジ時間 | 10,968 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 5 |
| other | AC * 27 |
ソースコード
class LazySegTree():
#区間和のセグメント木、要素ls
def __init__(self,ls):
self.default = 0
self.func = (lambda x, y: x + y)
self.N = len(ls)
self.K = (self.N-1).bit_length()
self.N0 = 1 << self.K
self.dat = [self.default]*(2**(self.K+1))
self.lazy = [0]*(2**(self.K+1)) #遅延評価
for i in range(self.N): #葉の構築
self.dat[2**self.K+i] = ls[i]
self.build()
def build(self):
for j in range(self.N0-1,0,-1):
self.dat[j] = self.func(self.dat[j<<1],self.dat[j<<1|1]) #親が持つ条件
def leafvalue(self,x): #x番目の値を出力
return self.query(x,x+1)
def update_add(self,x,y): #x番目にyを足す
return self.updatel_add(x,x+1,y)
def gindex(self,l, r): #伝播するインデックス列挙
L = l + self.N0; R = r + self.N0
lm = (L // (L & -L)) >> 1
rm = (R // (R & -R)) >> 1
lsindex = []
while L < R:
if R <= rm:
lsindex.append(R)
if L <= lm:
lsindex.append(L)
L >>= 1; R >>= 1
while L:
lsindex.append(L)
L >>= 1
return lsindex
def propagates(self,ids):
for i in reversed(ids):
v = self.lazy[i]
if not v:
continue
self.lazy[i<<1] += v//2; self.lazy[i<<1|1] += v//2
self.dat[i<<1] += v//2; self.dat[i<<1|1] += v//2
self.lazy[i] = 0
def updatel_add(self,l, r, x):#区間にyを足す
self.propagates(self.gindex(l, r))
L = self.N0 + l
R = self.N0 + r
ii = 0
while L < R:
if L & 1:
self.lazy[L] += x*(2**ii)
self.dat[L] += x*(2**ii)
L += 1
if R & 1:
R -= 1
self.lazy[R] += x*(2**ii)
self.dat[R] += x*(2**ii)
L >>= 1
R >>= 1
ii += 1
for i in self.gindex(l, r):
self.dat[i] = self.func(self.dat[i<<1],self.dat[i<<1|1])+self.lazy[i]
def query(self,l, r):
self.propagates(self.gindex(l, r))
L = self.N0 + l
R = self.N0 + r
vL = self.default
vR = self.default
while L < R:
if L & 1:
vL = self.func(vL,self.dat[L])
L += 1
if R & 1:
R -= 1
vR = self.func(self.dat[R],vR)
L >>= 1
R >>= 1
return self.func(vL,vR)
def modPow(a,n,mod):#繰り返し二乗法 a**n % mod
if n==0:
return 1
if n==1:
return a%mod
if n & 1:
return (a*modPow(a,n-1,mod)) % mod
t = modPow(a,n>>1,mod)
return (t*t)%mod
N = int(input())
lsA = list(map(int,input().split()))+[1]
mod = 10**9+7
if 0 in lsA:
print(0)
exit()
rh = [0]*(N)
l = 0
r = 0
val = 1
while l < N:
while val < 10**9 and r < N:
val *= lsA[r]
if val >= 10**9:
val //= lsA[r]
break
if r == N:
break
r += 1
rh[l] = r
val //= lsA[l]
l += 1
LSG = LazySegTree([0]*(N+1))
for i in range(N):
LSG.update_add(i, rh[i]-i)
LSG.updatel_add(i+1, rh[i]+1, -1)
ll = [LSG.leafvalue(i) for i in range(N+1)]
for i in range(1,N+1):
ll[i] += ll[i-1]
ans = 1
for i in range(N):
ans *= modPow(lsA[i],ll[i],mod)
ans %= mod
print(ans)