結果
問題 |
No.377 背景パターン
|
ユーザー |
|
提出日時 | 2021-06-15 00:17:32 |
言語 | Ruby (3.4.1) |
結果 |
AC
|
実行時間 | 2,529 ms / 5,000 ms |
コード長 | 1,093 bytes |
コンパイル時間 | 135 ms |
コンパイル使用メモリ | 7,424 KB |
実行使用メモリ | 12,544 KB |
最終ジャッジ日時 | 2024-12-25 11:54:34 |
合計ジャッジ時間 | 8,182 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
外部呼び出し有り |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 5 |
other | AC * 14 |
コンパイルメッセージ
Syntax OK
ソースコード
#!/usr/bin/ruby def yield_from(*funcsyms) funcsyms.each{|funcsym| func = method(funcsym) define_method(funcsym){|*args,&blk| return to_enum(funcsym,*args) if !blk func.(*args){|*a|a.each{|e|blk.(e)}} } } end if false require 'prime' def prime_division(n) n.prime_division end else def prime_division(n) IO.popen('factor '+n.to_s){|io| a=[] line=io.gets line.split[1..-1].group_by{|e|e}.each{|k,v|a<<[k.to_i,v.size]} a } end end def divisor_totient(a,d,n,t,&blk) return to_enum(:divisor,a,d,n,t) if !block_given? if d==a.size blk.call [n,t] else (0..a[d][1]).each{|i| divisor_totient( a,d+1, n*a[d][0]**i, i==0 ? t : t*(a[d][0]-1)*a[d][0]**(i-1), &blk ) } end end def pow(x,y,m) z=1 while y>0 z=z*x%m if y%2>0 x=x*x%m y/=2 end z end cache={} M=1000000007 H,W,K=gets.split.map(&:to_i) A=prime_division(H) B=prime_division(W) r=0 divisor_totient(A,0,1,1){|a,at| divisor_totient(B,0,1,1){|b,bt| key=W*H/a.lcm(b) cache[key]=pow(K,key,M) if !cache.has_key?(key) r=(r+at*bt*cache[key])%M } } p r*pow(W*H,M-2,M)%M