結果

問題 No.1084 積の積
ユーザー tyawanmusityawanmusi
提出日時 2020-06-13 15:42:27
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 363 ms / 2,000 ms
コード長 1,218 bytes
コンパイル時間 999 ms
コンパイル使用メモリ 10,684 KB
実行使用メモリ 19,224 KB
最終ジャッジ日時 2023-09-16 10:06:57
合計ジャッジ時間 8,023 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 16 ms
7,840 KB
testcase_01 AC 15 ms
7,832 KB
testcase_02 AC 16 ms
7,808 KB
testcase_03 AC 15 ms
7,872 KB
testcase_04 AC 363 ms
17,156 KB
testcase_05 AC 33 ms
10,052 KB
testcase_06 AC 360 ms
17,156 KB
testcase_07 AC 15 ms
8,192 KB
testcase_08 AC 16 ms
7,736 KB
testcase_09 AC 46 ms
18,140 KB
testcase_10 AC 15 ms
8,152 KB
testcase_11 AC 23 ms
8,528 KB
testcase_12 AC 110 ms
12,648 KB
testcase_13 AC 216 ms
17,712 KB
testcase_14 AC 89 ms
11,584 KB
testcase_15 AC 82 ms
11,308 KB
testcase_16 AC 242 ms
18,840 KB
testcase_17 AC 102 ms
12,456 KB
testcase_18 AC 170 ms
15,044 KB
testcase_19 AC 223 ms
17,644 KB
testcase_20 AC 63 ms
10,252 KB
testcase_21 AC 88 ms
11,744 KB
testcase_22 AC 73 ms
10,696 KB
testcase_23 AC 137 ms
13,868 KB
testcase_24 AC 125 ms
13,452 KB
testcase_25 AC 219 ms
17,640 KB
testcase_26 AC 262 ms
19,104 KB
testcase_27 AC 263 ms
18,500 KB
testcase_28 AC 270 ms
18,608 KB
testcase_29 AC 258 ms
19,040 KB
testcase_30 AC 265 ms
18,476 KB
testcase_31 AC 269 ms
19,224 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

# 問題を、変更しました！笑
#
# 解法メモ
# p[i] = (l=i)としたときのf(l,r)<10^9である(rの最大値)-i+1、とする
# つまりp[i]はA_iを左端とする総乗が10^9を超えない区間の個数
# これはしゃくとりで常勝
# q[i] = 答えに含まれるA_iの数 
# つまり、A_1^q[1]*A_2^q[2]*...*A_N^q[N]=ans
# pが求まると、少し実装は複雑だけどqを求められる
# p[i]=jとしたとき、q[k]+=j-(k-i)(i <= k < i+j)
# 例えば、p[i]=5のときq[i]+=5,q[i+1]+=4,q[i+2]+=3,...q[i+4]+=1
# これはいい感じに実装するとO(N)で操作できる
# 私は天才なのでimosを使った天才実装をします
# 末尾とかが怖いけれど、両端に10^9を加えるなどで対処可能

n=int(input())
a=list(map(int,input().split()))+[10**9]
mod=10**9+7
if 0 in a:exit(print(0))
p=[0]*n
s=1
j=0
for i in range(n):
  while s*a[j]<10**9:
    s*=a[j]
    j+=1
  p[i]=j-i
  s//=a[i]
imos=[0]*(n+1)
for i in range(n):
  imos[i]+=1
  imos[i+p[i]]-=1
for i in range(n):
  imos[i+1]+=imos[i]
for i in range(n):
  if i!=0:imos[i-1]-=p[i]
for i in range(n,0,-1):
  imos[i-1]+=imos[i]
ans=1
for i in range(n):
  ans*=pow(a[i],imos[i],mod)
  ans%=mod
print(ans)
0