結果

問題 No.16 累乗の加算
ユーザー jj
提出日時 2016-09-12 21:56:04
言語 Fortran
(gFortran 14.2.0)
結果
AC  
実行時間 7 ms / 5,000 ms
コード長 796 bytes
コンパイル時間 1,877 ms
コンパイル使用メモリ 32,384 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-26 05:16:35
合計ジャッジ時間 2,567 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 14
権限があれば一括ダウンロードができます

ソースコード

diff # 

program main
  implicit none
  integer*8::x,N,total,i
  integer*8,allocatable::a(:)
  integer::memo(0:100000000)
  integer,parameter::modulo=1000003
  data memo/100000001*0/,total/0/
  read *,x,N
  allocate(a(N))
  read *,a
  memo(0) = 1
  memo(1) = x
  do i=1,N
     total = total + power(x,a(i),memo)
  end do
  print '(i0)',MOD(total,modulo)
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(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 function power
end program
0