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;
}