hdu 4790 Just Random

思路：對於a<=x<=b,c<=y<=d，知足條件的結果爲ans=f(b,d)-f(b,c-1)-f(a-1,d)+f(a-1,c-1)。

而函數f(a,b)是計算0<=x<=a,0<=y<=b知足條件的結果。這樣計算就很方便了。

例如：求f(16,7),p=6,m=2.

對於x有：0 1 2 3 4 5 0 1 2 3 4 5 0 1 2 3 4

對於y有：0 1 2 3 4 5 0 1

很容易知道對於xy中的（0 1 2 3 4 5）對知足條件的數目爲p。

這樣取A集合爲（0 1 2 3 4 5 0 1 2 3 4 5），B集合爲（0 1 2 3 4）。

C集合爲（0 1 2 3 4 5），D集合爲（0 1）。

這樣就能夠分紅4部分來計算了。

f(16,7)=A和C知足條件的數+A和D知足條件的數+B和C知足條件的數+B和D知足條件的數。

其中前3個很好求的，關鍵是B和D知足條件的怎麼求！

這個要根據m來分狀況。

#include<cstdio>
#include<iostream>
using namespace std;
#define ll __int64
ll p,m;
ll gcd(ll a,ll b)
{
    if(a<b) swap(a,b);
    while(b){
        ll t=a;
        a=b;
        b=t%b;
    }
    return a;
}
ll f(ll a,ll b)
{
    if(a<0||b<0) return 0;
    ll ma=a%p,mb=b%p,ans;
    ans=(a/p)*(b/p)*p;//1
    ans+=(ma+1)*(b/p)+(mb+1)*(a/p);//2+3
    if(ma>m){ //4
        ans+=min(m,mb)+1;
        ll t=(p+m-ma)%p;//根據ma求出知足最小的y來
        if(t<=mb) ans+=mb-t+1;
    }else{
        ll t=(p+m-ma)%p;//根據ma求出知足最小的y來
        if(t<=mb) ans+=min(m-t+1,mb-t+1);
    }
    return ans;
}
int main()
{
    int ca=0,t;
    ll a,b,c,d;
    scanf("%d",&t);
    while(t--){
        cin>>a>>b>>c>>d>>p>>m;
        ll ans=f(b,d)-f(b,c-1)-f(a-1,d)+f(a-1,c-1);
        ll tot=(b-a+1)*(d-c+1);
        ll g=gcd(ans,tot);
        printf("Case #%d: %I64d/%I64d\n",++ca,ans/g,tot/g);
    }
    return 0;
}

View Code