题目描述
JerryC对数字痴迷到了一种非正常的境界。每天JerryC都有喜欢的一些数字。第 iii 天JerryC就喜欢Ai−BiA_i-B_iAi−Bi中的数字。但是他觉得这样并不是很有趣,于是他就定义了一个函数
f(i)=i∗∑j=1i(i mod j)f(i)=i*\sum^{i}_{j=1}(i\ mod\ j)f(i)=i∗j=1∑i(i mod j)
但是JerryC还是觉得不够有趣,于是她觉得只有满足
∑i=ABf(i)\sum^{B}_{i=A}f(i)i=A∑Bf(i)
为奇数的一天才是她特别高兴的一天。
现在告诉你最近TTT天JerryC喜欢的数的区间,请你猜猜她每一天会不会特别高兴。
解释:
第一个式子就是表示 f(i)f(i)f(i) 等于 iii 乘以 iii 的约数个数。
第二个式子就是求出 [A,B][A, B][A,B] 区间内所有的 f(i)f(i)f(i) 的和。
输入输出格式
输入格式:
第111行一个整数TTT表示有多少天。
第2−(T+1)2-(T+1)2−(T+1)行两个非负整数Ai, BiA_i,\ B_iAi, Bi表示每天JerryC喜欢的数的区间。
输出格式:
TTT行。每行为"YesYesYes“或”NoNoNo",表示第 iii 天是不是JerryC特别高兴的一天。
输入输出样例
输入样例:
333
1 101\ 101 10
11 1511\ 1511 15
5 125\ 125 12
#####输出样例:
NoNoNo
NoNoNo
YesYesYes
数据范围
0<=T<=106, 1<=Ai<=Bi<=10180<=T<=10^6,\ 1<=A_i<=B_i<=10^{18}0<=T<=106, 1<=Ai<=Bi<=1018
Solution\text{Solution}Solution
我们发现,只有f(c2), cf(c^2),\ cf(c2), c是奇数的情况才能对答案产生贡献。设c=⌈a⌉, d=bc=\lceil\sqrt a\rceil,\ d=\sqrt bc=⌈a⌉, d=b。
对于区间[a,b][a,b][a,b],只有[c,d][c,d][c,d]区间长度中的奇数对答案有贡献。我们求出[c,d][c,d][c,d]区间长度便可以求出答案。
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#define reg register
typedef long long ll;
int n,ans;
ll a,b;
ll c,d;
//读入优化
inline ll read(){
ll x=0;char c;
do c=getchar(); while(c<'0'||c>'9');
while(c>='0'&&c<='9')
x=x*10+c-48,c=getchar();
return x;
}
int main(){
scanf("%d",&n);
for(reg int i=1;i<=n;++i){
a=read(),b=read();
c=ceil(sqrt(a)),d=sqrt(b);
if(a==b)//特判
{
printf(c*c==a&&(c&1)?"Yes\n":"No\n");
continue;
}
//如果c是偶数,他对答案不产生贡献
if(!(c&1)) c++;
printf(((d-c)/2+1)&1?"Yes\n":"No\n");
}
}