結果

問題 No.377 背景パターン
ユーザー ciel
提出日時 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

ソースコード

diff #
プレゼンテーションモードにする

#!/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
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0