結果
問題 | No.377 背景パターン |
ユーザー | ciel |
提出日時 | 2016-06-05 16:41:03 |
言語 | PyPy2 (7.3.15) |
結果 |
AC
|
実行時間 | 3,159 ms / 5,000 ms |
コード長 | 868 bytes |
コンパイル時間 | 326 ms |
コンパイル使用メモリ | 76,448 KB |
実行使用メモリ | 98,176 KB |
最終ジャッジ日時 | 2024-04-17 04:57:39 |
合計ジャッジ時間 | 10,909 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
外部呼び出し有り |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 181 ms
90,496 KB |
testcase_01 | AC | 185 ms
90,776 KB |
testcase_02 | AC | 185 ms
90,496 KB |
testcase_03 | AC | 186 ms
90,496 KB |
testcase_04 | AC | 180 ms
90,624 KB |
testcase_05 | AC | 179 ms
90,624 KB |
testcase_06 | AC | 175 ms
90,624 KB |
testcase_07 | AC | 179 ms
90,624 KB |
testcase_08 | AC | 175 ms
90,880 KB |
testcase_09 | AC | 179 ms
90,880 KB |
testcase_10 | AC | 181 ms
90,496 KB |
testcase_11 | AC | 178 ms
90,880 KB |
testcase_12 | AC | 176 ms
90,624 KB |
testcase_13 | AC | 242 ms
92,160 KB |
testcase_14 | AC | 231 ms
91,904 KB |
testcase_15 | AC | 180 ms
90,632 KB |
testcase_16 | AC | 3,026 ms
97,040 KB |
testcase_17 | AC | 3,159 ms
98,176 KB |
testcase_18 | AC | 179 ms
90,880 KB |
ソースコード
#!/usr/bin/python from collections import defaultdict from fractions import gcd import subprocess import sys if sys.version_info[0]>=3: raw_input=input xrange=range def prime_division(n): h=defaultdict(int) line=subprocess.check_output(['factor',str(n)]) for e in line.split()[1:]: h[int(e)]+=1 return [(k,v) for k,v in h.items()] def divisor_totient(a,d,n,t): if d==len(a): yield (n,t) else: for i in xrange(a[d][1]+1): for x,xt in divisor_totient( a,d+1, n*a[d][0]**i, (t if i==0 else t*(a[d][0]-1)*a[d][0]**(i-1)) ): yield (x,xt) cache={} M=1000000007 H,W,K=map(int,raw_input().split()) A=prime_division(H) B=prime_division(W) r=0 for a,at in divisor_totient(A,0,1,1): for b,bt in divisor_totient(B,0,1,1): key=W*H*gcd(a,b)//a//b if key not in cache: cache[key]=pow(K,key,M) r=(r+at*bt*cache[key])%M print(r*pow(W*H,M-2,M)%M)