hdu 4790 Just Random

hdu 4790 Just Random

题意：
给定两个区间，求在两个区间中选出两个数的和模p为m的概率。

 

这个题目是个数学问题，模拟下就会发现选出的数呈有规律的矩阵分布，假设a+d大于b+c，x与y的和从a+c到b+c其数量依次递增，直接利用等差数列的求和公式计算就好了，从b+c到a+d其数量都为b-a+1，直接统计个数再相乘就行了，从a+d到b+d其数量依次递减，也是等差数列的求和。注意这是个区间，需要注意每个部分判断时的细节问题。

 

比赛做这个题的时候，方向本来已经对了，知道了数量的规律分布，但是没有多想就放弃了，后面比赛完了才知道真是这么做的，为何我不能深入思考呢，很多题目看似复杂，但我总是没有静下心来一步一步分析题目….必须有深搜终点的觉悟才行！

 

代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
#define INF 0x3f3f3f3f
#define LL long long

using namespace std;

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

int main()
{
    LL a, b, c, d, p, m, sum, n, i, tot, k;
    int t, j = 0;
    scanf("%d", &t);
    while (t--){
        sum = 0;
        scanf("%I64d %I64d %I64d %I64d %I64d %I64d", &a, &b, &c, &d, &p, &m);
        if (b+c > a+d) swap(a, c), swap(b, d);
        n = (b+c-1-m)/p;
        if (b+c-1 < m) n = -1;
        i = (a+c-m)/p;
        if ((a+c-m) %p != 0) i++;
        if (a+c <= m) i = 0;
        if (n>=0 &&  n>=i){
            k = i*p+m-a-c+1;
            sum = sum+(n-i+1)*(2*k+(n-i)*p)/2;
        }
        n++; i = (a+d-m)/p;
        if (a+d < m) i = -1;
        if (i >= n) sum = sum+(b-a+1)*(i-n+1);
        n = (b+d-m)/p; i++;
        if (b+d < m) n = -1;
        if (n >= i){
            k = b-a-(i*p+m-a-d-1);
            sum = sum+(n-i+1)*(2*k-(n-i)*p)/2;
        }
        tot = (b-a+1)*(d-c+1);
        LL tmp = gcd(tot, sum);
        printf("Case #%d: %I64d/%I64d\n", ++j, sum/tmp, tot/tmp);
    }
    return 0;
}