結果

問題 No.886 Direct
ユーザー tempura_pp
提出日時 2019-08-20 18:44:25
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,355 ms / 4,000 ms
コード長 609 bytes
コンパイル時間 426 ms
コンパイル使用メモリ 81,972 KB
実行使用メモリ 319,276 KB
最終ジャッジ日時 2024-10-06 12:45:08
合計ジャッジ時間 14,827 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 32
権限があれば一括ダウンロードができます

ソースコード

diff # 

h, w = map(int,input().split())
n = max(h+1, w+1)
cnt1 = [max(0,h-i) for i in range(n)]
cnt2 = [max(0,w-i) for i in range(n)]
isprime = [1]*n
isprime[0]=0
isprime[1]=0
prime = []
for i,b in enumerate(isprime):
    if b == 0: continue
    prime.append(i)
    j = 2*i
    while j < n :
        isprime[j] = 0
        j += i
for p in prime :
    m = (n-1)//p
    for i in range(m,0,-1):
        cnt1[i]+=cnt1[i*p]
        cnt2[i]+=cnt2[i*p]
cnt = [cnt1[i]*cnt2[i] for i in range(n)]
for p in prime :
    m = (n-1)//p
    for i in range(1,m+1):
        cnt[i]-=cnt[i*p]
print((2*cnt[1]+h*(w-1)+w*(h-1))%(10**9+7))
0