hdu 4790 Just Random （2013成都J题） 数学思路题 容斥

题意：在[a,b]  [c,d] 之间，和模p等于m的对数

详见代码

 #include <stdio.h>
 #include <algorithm>
 #include <string.h>
 #include<cmath>
 #define LL long long
 using namespace std;
 int T;
 LL a,b,c,d,p,m;

 LL gcd(LL a, LL b) {
     return b ? gcd(b, a % b) : a;
 }

 LL fun(LL x,LL y) {//表示0到x区间，0到y区间的组合对数
     LL ret;
     LL ra,rb;
     ra=x%p,rb=y%p;
     ret=(x/p)*(y/p)*p;
     ret+=(ra+)*(y/p)+(rb+)*(x/p);
     if(ra>m) {
         ret+=min(m+,rb+);
         LL tmp=(m+p-ra)%p;
         if(tmp<=rb) {
             ret+=rb-tmp+;
         }
     } else {
         LL tmp=(m+p-ra)%p;
         if(tmp<=rb) {
             ret+=min(m-tmp+,rb-tmp+);
         }
     }
     return ret;
 }
 int main() {
     scanf("%d",&T);
     int cas=;
     while(T--) {
         cas++;
         scanf("%I64d%I64d%I64d%I64d%I64d%I64d",&a,&b,&c,&d,&p,&m);
         LL ans=fun(b,d)-fun(a-,d)-fun(b,c-)+fun(a-,c-);//容斥原理求和
         LL tot=(b-a+)*(d-c+);
         if(ans<)
             ans=;
 //        printf("ans:%I64d tot:%I64d\n",ans,tot);
         LL l=gcd(ans,tot);
         ans/=l,tot/=l;
         printf("Case #%d: %I64d/%I64d\n", cas, ans, tot);
     }
     return ;
 }