非常好题目,使我思考良多。
经典的数学问题,题目给出了两道数学式子,当然就是要让我们推啊。
首先纠正一个错误,当我发现了这两个式子能凑出来平方式的时候我就想着去用两个式子来互相简化,但其实这样存在以下一个错误:
当我们有一个对于两个正整数的式子:a2<x2a^2 < x^2a2<x2 的时候,假如我们对等式两边进行平方,就会得到 a2<x2a^2 < x^2a2<x2 ,这个式子当然是没问题的,但是会导致取到的范围精度变化了,也就是说中间范围不准了,对我们这道题是不适用的。
那么我们就直接从这两个式子上下手。
如果我们要直接枚举a,可以知道a的范围是 1 ~ n/2 ,因为假如使得b和c都取到最小值,即1,那么就能够通过式子 a∗b−a∗c−b∗c≤na*b - a*c - b*c \le na∗b−a∗c−b∗c≤n得到a的范围。(由于a,b,c都是正整数,所以n一定是大于等于x的,所以不必拿x的石子来取范围)。
知道了a的区间,我们来思考b的。
我们至少知道:a∗b≤na*b \le na∗b≤n,所以b的取值范围就是na\frac{n}{a}an,这里不妨思考一下取值范围,当a从1取到n/2时,b就从n1\frac{n}{1}1n取到了nn2\frac{n}{\frac{n}{2}}2nn是一个调和级数,大约是nlognnlognnlogn级别的。
这时候重点来了,我们怎么确定c?
如果你想要直接通过循环来找c,那你就寄了,因为指定会爆。
这时候仔细思考一下题目,我们要求的是整个小于等于一个范围内的数,所以只要是我们找到了一个最大值,那c的取值范围可以是1 ~ 最大值
,那么我们直接加上最大值就好了,也就省去了一层循环,这样一定能过。
c的最大值当然是通过题目给出的题干来求。
通过第一个式子,我们可以得出:c≤n−a∗ba+bc \le \frac{n - a*b}{a+b}c≤a+bn−a∗b
通过第二个式子,我们可以得出:c≤x−a−bc \le x - a - bc≤x−a−b
那么直接从二者中取最小值就可以了,值得注意的是,当其中有最小值取成了负数的时候,就代表了没有可用的c,直接跳过就可以了。
CODE:
#include<bits/stdc++.h>
using namespace std;
const int N = 1e6 + 10;
void solve(){
int n,x;
cin >> n >> x;
int res = 0;
for(int a = 1;a <= n/2;a++){
for(int b = 1;b <= n/a;b++){
if(a*b > n)continue;
if(x - a - b < 0)continue;
res += min((n - a*b)/(a+b),x - a - b);
}
}
cout << res << endl;
}
int main(){
int T;cin >> T;
while(T--){
solve();
}
return 0;
}