結果
| 問題 | No.16 累乗の加算 |
| コンテスト | |
| ユーザー |
jj
|
| 提出日時 | 2016-09-12 21:53:33 |
| 言語 | Fortran (gFortran 15.2.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 876 bytes |
| 記録 | |
| コンパイル時間 | 299 ms |
| コンパイル使用メモリ | 32,384 KB |
| 実行使用メモリ | 6,016 KB |
| 最終ジャッジ日時 | 2024-11-17 04:06:03 |
| 合計ジャッジ時間 | 1,020 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 4 WA * 10 |
ソースコード
program main
implicit none
integer*8::x,N,total,i
integer*8,allocatable::a(:)
integer::memo(0:100000000)
data memo/100000001*0/,total/0/
read *,x,N
allocate(a(N))
read *,a
do i=1,N
total = total + power(x,a(i),memo)
end do
print '(i0)',total
contains
recursive function power(x,p,memo) result(y)
integer*8::x,y
integer*8::p
integer,parameter::modulo=1000003
integer::memo(0:100000000)
if(memo(p).ne.0) then
y = memo(p)
else if(p.eq.0) then
y = 1
memo(p) = y
else if(p.eq.1) then
y = x
memo(p) = y
else
if(MOD(p,2).eq.0) then
y = mod(power(x,p/2,memo)*power(x,p/2,memo),modulo)
memo (p) = y
else
y = mod(power(x,p/2+1,memo)*power(x,p/2,memo),modulo)
memo (p) = y
end if
end if
end function power
end program main
jj