結果

問題 No.36 素数が嫌い！
ユーザー jj
提出日時 2016-10-10 17:16:45
言語 Fortran
(gFortran 14.2.0)
結果
AC  
実行時間 163 ms / 5,000 ms
コード長 1,327 bytes
コンパイル時間 1,456 ms
コンパイル使用メモリ 34,288 KB
実行使用メモリ 77,952 KB
最終ジャッジ日時 2024-06-27 00:49:52
合計ジャッジ時間 3,529 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

program main
  implicit none
  integer*8::N,upper,i
  integer*8,allocatable::base(:),exponents(:)
  read *,N
  upper = DINT(log(DBLE(N)))+1_8
  allocate(base(upper),exponents(upper))
  base = 0
  exponents = 0
  call factorization(N, base, exponents)

contains
  subroutine factorization(n, base, exponents)
    integer*8::n,base(:),exponents(:)
    integer*8::i,j,itr,upper,total_exponents
    integer*8,allocatable::furui(:)
    data itr/0/, total_exponents/0/
    upper = INT(SQRT(DBLE(n)+1))
    allocate(furui(2:upper))
    furui = 0

    do i=2,upper
       if(furui(i).ne.0) cycle
       if(MOD(n, i).ne.0) cycle
       itr = itr + 1
       base(itr) = i
       do j=i+i,upper,i
          furui(j) = 1
       end do

       do j=1,upper
          n = n/i
          if(MOD(n,i).ne.0) exit
       end do
       exponents(itr) = j
       total_exponents = total_exponents + j
       if(n.ne.1) then
          if(total_exponents.ge.2) then
             print '(a)',"YES"
             return
          end if
       else
          if(total_exponents.ge.3) then
             print '(a)',"YES"
          else
             print '(a)',"NO"
          end if
          return
       end if
    end do
    itr = itr + 1
    base(itr) = n
    exponents(itr) = 1
    print '(a)', "NO"
  end subroutine factorization
end program main
0