UVA-10200 Prime Time 素数（注意除法精度）

题目链接：https://cn.vjudge.net/problem/UVA-10200

题意

给出一个公式$ m=n^2+n+41, n \in Z^+ $
现在$ a,b\in[0, 10000] $的范围内取n
问有几个m是素数

思路

不说了
关键是注意除法精度问题
当出现不得不使用除法时，一定要在除法的结果上加1e6，保证精度

代码

#include <cstdio>
#include <cmath>
const double eps=1e-6;
int sum[int(1e4)+15];
bool isprime(const long long &i){
long long n=i*i+i+41, size=sqrt(n)+1;
for (int mod=2; mod!=n && mod<=size; mod++)
if (n%mod==0) return false;
return true;
}

void init(void){
sum[0]=1;
for (int i=0; i<=int(1e4); i++)
if (isprime(i)) sum[i+1]=sum[i]+1;
else sum[i+1]=sum[i];
}

int main(void){
int a, b;
init();

while (scanf("%d%d", &a, &b)==2)
printf("%.2f\n", 100*(sum[b+1]-sum[a])/(double)(b-a+1)+eps);

return 0;
}