ZOJ - 3483 - Gaussian Prime

先上题目：

Gaussian Prime

Time Limit: 3 Seconds Memory Limit: 65536 KB

In number theory, a Gaussian integer is a complex number whose real and imaginary part are both integers. The Gaussian integers, with ordinary addition and multiplication of complex numbers, form an integral domain, usually written as Z[i]. The prime elements of Z[i] are also known as Gaussian primes. Gaussian integers can be uniquely factored in terms of Gaussian primes up to powers of i and rearrangements.

A Gaussian integer a + bi is a Gaussian prime if and only if either:

One of a, b is zero and the other is a prime number of the form 4n + 3 (with n a nonnegative integer) or its negative -(4n + 3), or
Both are nonzero and a² + b² is a prime number (which will not be of the form 4n + 3).

0 is not Gaussian prime. 1, -1, i, and -i are the units of Z[i], but not Gaussian primes. 3, 7, 11, ... are both primes and Gaussian primes. 2 is prime, but is not Gaussian prime, as 2 =i(1-i)².

Your task is to calculate the density of Gaussian primes in the complex plane [x1, x2] × [y1, y2]. The density is defined as the number of Gaussian primes divided by the number of Gaussian integers.

Input

There are multiple test cases. The first line of input is an integer T ≈ 100 indicating the number of test cases.

Each test case consists of a line containing 4 integers -100 ≤ x1 ≤ x2 ≤ 100, -100 ≤ y1 ≤ y2 ≤ 100.

Output

For each test case, output the answer as an irreducible fraction.

Sample Input

Sample Output

References

http://en.wikipedia.org/wiki/Gaussian_integer
Weisstein, Eric W. "Gaussian Prime." From MathWorld--A Wolfram Web Resource. http://mathworld.wolfram.com/GaussianPrime.html

　　题意：告诉你一种数的定义，这种数是一个复数，告诉你题目的整个区间的范围，然后给你一个区间，问这个区间里面这种数的密度是多少(这种数比上区间里面的数的总个数)？

　　题目给的范围比较小，所以可以先预处理把区间里面的这种数都标记出来，然后对于每一次询问就搜一次。

　　这题需要注意的地方是这种数的定义。这里现需要把0~20000的素数都筛出来，然后根据定义把区间的这种数都找出来就行了。

　　还有一个需要注意的地方是如果分子是零的时候需要输出的是0/1，而不是0/大于1的分母。

上代码：

 #include <cstdio>

 #include <cstring>

 #define MAX 20002

 #define LL long long

 using namespace std;

 bool f[MAX];

 void deal(){

     LL n=MAX-;

     f[]=f[]=;

     memset(f,,sizeof(f));

     for(LL i=;i<=n;i++){

         if(!f[i]){

             for(LL j=i*i;j<=n;j+=i){

                 f[j]=;

             }

         }

     }

 }

 bool s[][];

 bool check(int y,int x){

     if(y== && x==) return ;

     else if((y== && x!= )|| (y!= && x==)){

         int k=x+y;

         if(k<) k=-k;

         if(k%==%){

             return !f[k];

         }

         return ;

     }else{

         LL sum=y*y+x*x;

         if(!f[sum] && (sum-+)%!=) return ;

         //if(!f[sum]) return 1;

     }

     return ;

 }

 void work(){

     int y,x;

     for(int i=;i<=;i++){

         for(int j=;j<=;j++){

             y=i-;

             x=j-;

             if(check(y,x)) s[i][j]=;

         }

     }

 }

 int gcd(int a,int b){

     return b== ? a : gcd(b,a%b);

 }

 void ask(){

     int a,b,c,d,g;

     int pr,num;

     scanf("%d %d %d %d",&a,&b,&c,&d);

     a+=;

     b+=;

     c+=;

     d+=;

     pr=;

     num=;

     for(int i=a;i<=b;i++){

         for(int j=c;j<=d;j++){

             if(s[i][j]) pr++;

             num++;

         }

     }

     g=gcd(num,pr);

     printf("%d/%d\n",pr/g,num/g);

 }

 int main()

 {

     int t;

     //freopen("data.txt","r",stdin);

     deal();

     work();

     scanf("%d",&t);

     while(t--){

         ask();

     }

     return ;

 }

3483

巴特西