cf962：D.Fun（三元）

题目

给定两个整数 𝑛n 和 𝑥x ，求 𝑎𝑏+𝑎𝑐+𝑏𝑐≤𝑛ab+ac+bc≤n 和 𝑎+𝑏+𝑐≤𝑥a+b+c≤x 的个正整数的三联数( 𝑎,𝑏,𝑐a,b,c )的个数。

注意顺序问题(例如 ( 1,1,21,1,2 ) 和 ( 1,2,11,2,1 ) 被视为不同)， 𝑎a ， 𝑏b ， 𝑐c 必须严格大于 00 。

输入

第一行包含一个整数 𝑡t ( 1≤𝑡≤1041≤t≤104 ) - 测试用例数。

每个测试用例包含两个整数 𝑛n 和 𝑥x ( 1≤𝑛,𝑥≤1061≤n,x≤106 )。( 1≤𝑛,𝑥≤1061≤n,x≤106 ).

保证所有测试用例中 𝑛n 的总和不超过 106106 ，所有测试用例中 𝑥x 的总和不超过 106106 。

输出

输出一个整数 - 𝑎𝑏+𝑎𝑐+𝑏𝑐≤𝑛ab+ac+bc≤n 和 𝑎+𝑏+𝑐≤𝑥a+b+c≤x 的正整数三元组( 𝑎,𝑏,𝑐a,b,c )的个数。

做法

这题就直接枚举a，b，然后求c的最大值中的最小值（n和x分别可以求一个最小值）就好了。

为什么两重循环1e6的范围不会爆呢，因为枚举第二重b的时间复杂度是log(n)的。官方解释是这样的：因为 ab+ac+bc≤n ，我们至少知道 ab≤n 。除以两边得到 b≤na 。当 a=1 时， b 有 n 个选择。当 a=2 时， b 有 n/2 个选择。因此， b 总共有 n+n/2+n/3+...+n/n 个选项。这仅仅是调和数列。因此总复杂度是n*log(n)

#include<bits/stdc++.h> 
using namespace std;
int t;
int n,x;
int main(){
	scanf("%d",&t);
	while(t--){
		scanf("%d%d",&n,&x);
		long long ans=0;
		for(int i=1;i<=x-2;i++){//枚举a
			for(int j=1;j<=x-2;j++) {//枚举b
				if(i*j>=n) break;
				if(i+j>=x) break;
				int k=(n-i*j)/(i+j);//算出c
				ans+=min(k,x-i-j);//c的最大值的最小值 
			}
		}
		cout<<ans<<endl; 
	}
}