結果

問題 No.377 背景パターン
ユーザー cielciel
提出日時 2017-12-31 22:42:23
言語 Crystal
(1.11.2)
結果
CE  
(最新)
AC  
(最初)
実行時間 -
コード長 901 bytes
コンパイル時間 1,979 ms
コンパイル使用メモリ 218,476 KB
最終ジャッジ日時 2024-04-27 02:30:42
合計ジャッジ時間 2,386 ms
ジャッジサーバーID
(参考情報)
judge2 / judge4
このコードへのチャレンジ
(要ログイン)
コンパイルエラー時のメッセージ・ソースコードは、提出者また管理者しか表示できないようにしております。(リジャッジ後のコンパイルエラーは公開されます)
ただし、clay言語の場合は開発者のデバッグのため、公開されます。

コンパイルメッセージ
Showing last frame. Use --error-trace for full trace.

In Main.cr:39:20

 39 | cache[key]=pow(k,key,M) if !cache.has_key?(key)
                       ^--
Error: expected argument #2 to 'pow' to be Int64, not Float64

Overloads are:
 - pow(x : Int64, y : Int64, m : Int64)

ソースコード

diff #

M=1000000007_i64
def pow(x : Int64,y : Int64,m : Int64)
	z=1_i64
	while y>0
		z=z*x%m if y%2>0
		x=x*x%m
		y/=2
	end
	z
end
def prime_division(n)
	Process.run("factor",["factor",n.to_s]){|io|
		h=Hash(Int64,Int64).new(0_i64)
		io.output.gets.not_nil!.split()[1..-1].each{|e|h[e.to_i64]+=1}
		h.map{|k,v|[k,v]}
	}
end
def divisor_totient(a,d,n,t,&blk : (Int64,Int64)->_)
	if d==a.size
		blk.call(n,t)
	else
		(0..a[d][1]).each{|i|
			divisor_totient(
				a,d+1,
				n*pow(a[d][0],i.to_i64,M),
				i==0 ? t : t*(a[d][0]-1)*pow(a[d][0],(i-1).to_i64,M),&blk)
		}
	end
end

cache=Hash(Int64,Int64).new
h,w,k=ARGF.gets_to_end.split.map(&.to_i64)
a=prime_division(h)
b=prime_division(w)
r=0_i64
divisor_totient(a,0,1_i64,1_i64){|ad,at|
	divisor_totient(b,0,1_i64,1_i64){|bd,bt|
		key=w*h/ad.lcm(bd)
		cache[key]=pow(k,key,M) if !cache.has_key?(key)
		r=(r+at*bt%M*cache[key])%M
	}
}
puts r*pow(w*h%M,M-2,M)%M
0