Submission #986856

---

Source Code Expand

program main
  implicit none
  integer(8) r,n,m,i
  real(8),allocatable::a(:)
  read*,r,n,m
  allocate(a((n/2+1)))
  a(:)=dsqrt(n*n*16d0-((/(i*i*16,i=n,0,-2)/)))
  if(mod(n+m+1,2).eq.1)a((n-m)/2+1)=a((n-m)/2+1)*1.5d0
    a((n-m)/2+2:)=a((n-m)/2+2:)*2d0
  print*,dble(sum(a)/dble(n)-4d0*mod(n+1,2))*r
end program main

Submission Info

Judge Result