結果

問題 No.1084 積の積
ユーザー kohei2019kohei2019
提出日時 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
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 42 ms
52,864 KB
testcase_01 AC 42 ms
52,608 KB
testcase_02 AC 42 ms
52,992 KB
testcase_03 AC 42 ms
52,480 KB
testcase_04 AC 959 ms
96,256 KB
testcase_05 AC 62 ms
78,720 KB
testcase_06 AC 1,017 ms
96,512 KB
testcase_07 AC 42 ms
52,480 KB
testcase_08 AC 41 ms
52,992 KB
testcase_09 AC 68 ms
83,328 KB
testcase_10 AC 41 ms
52,736 KB
testcase_11 AC 135 ms
77,696 KB
testcase_12 AC 217 ms
85,120 KB
testcase_13 AC 343 ms
96,000 KB
testcase_14 AC 196 ms
81,536 KB
testcase_15 AC 195 ms
81,536 KB
testcase_16 AC 379 ms
98,048 KB
testcase_17 AC 216 ms
84,352 KB
testcase_18 AC 277 ms
90,112 KB
testcase_19 AC 349 ms
95,872 KB
testcase_20 AC 172 ms
79,744 KB
testcase_21 AC 200 ms
81,536 KB
testcase_22 AC 188 ms
80,256 KB
testcase_23 AC 244 ms
87,040 KB
testcase_24 AC 237 ms
85,632 KB
testcase_25 AC 339 ms
95,616 KB
testcase_26 AC 457 ms
96,384 KB
testcase_27 AC 456 ms
96,640 KB
testcase_28 AC 456 ms
96,384 KB
testcase_29 AC 446 ms
96,384 KB
testcase_30 AC 458 ms
96,384 KB
testcase_31 AC 444 ms
96,384 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

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)
0