UVA 10200 Prime Time (打表)
Prime Time
Euler is a well-known matematician, and, among many other things, he discovered that the formula n 2 + n + 41 produces a prime for 0 ≤ n < 40. For n = 40, the formula produces 1681, which is 41 ∗ 41. Even though this formula doesn’t always produce a prime, it still produces a lot of primes. It’s known that for n ≤ 10000000, there are 47,5% of primes produced by the formula! So, you’ll write a program that will output how many primes does the formula output for a certain interval.
Input
Each line of input will be given two positive integer a and b such that 0 ≤ a ≤ b ≤ 10000. You must read until the end of the file.
Output
For each pair a, b read, you must output the percentage of prime numbers produced by the formula in this interval (a ≤ n ≤ b) rounded to two decimal digits.
Sample Input
0 39
0 40
39 40
Sample Output
100.00
97.56
50.00
n可根据公式 f(n) = n*n+n+41得到f(n),f(n)可能是素数也可能不是素数,求a到b之间的数可根据公式得到素数的百分比
#include<stdio.h>
#include<iostream>
#include<algorithm>
#include<string.h>
using namespace std;
int a[10010];
int judge(int x)
{for(int i=2;i*i<=x;i++)
if(x%i==0)
return 0;
return 1;
}
void init()
{a[0]=0;
for(int i=0;i<=10005;i++)
{if(judge(i*i+i+41)==1)
a[i+1]=a[i]+1;
else
a[i+1]=a[i];
}
}
int main()
{
int c,b;
init();
while(~scanf("%d%d",&c,&b))
{float p=(float)(a[b+1]-a[c])/(float)(b-c+1)*100.0;
printf("%.2f\n",p+0.0000001);//坑在这了
}
}