結果
問題 | No.377 背景パターン |
ユーザー |
|
提出日時 | 2021-06-15 00:16:15 |
言語 | Ruby (3.4.1) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,109 bytes |
コンパイル時間 | 49 ms |
コンパイル使用メモリ | 7,552 KB |
実行使用メモリ | 24,960 KB |
最終ジャッジ日時 | 2024-12-25 11:51:40 |
合計ジャッジ時間 | 15,285 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 5 |
other | AC * 12 TLE * 2 |
コンパイルメッセージ
Main.rb:31: warning: `*' interpreted as argument prefix Syntax OK
ソースコード
#!/usr/bin/rubydef yield_from(*funcsyms)funcsyms.each{|funcsym|func = method(funcsym)define_method(funcsym){|*args,&blk|return to_enum(funcsym,*args) if !blkfunc.(*args){|*a|a.each{|e|blk.(e)}}}}endif falserequire 'prime'def prime_division(n) n.prime_division endelsedef prime_division(n)IO.popen('factor '+n.to_s){|io|a=[]line=io.getsline.split[1..-1].group_by{|e|e}.each{|k,v|a<<[k.to_i,v.size]}a}endenddef divisor_totient(a,d,n,t)return to_enum(:divisor,a,d,n,t) if !block_given?if d==a.sizeyield [n,t]else(0..a[d][1]).each{|i|yield *divisor_totient(a,d+1,n*a[d][0]**i,i==0 ? t : t*(a[d][0]-1)*a[d][0]**(i-1))}endendyield_from :divisor_totientdef pow(x,y,m)z=1while y>0z=z*x%m if y%2>0x=x*x%my/=2endzendcache={}M=1000000007H,W,K=gets.split.map(&:to_i)A=prime_division(H)B=prime_division(W)r=0divisor_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