[cf1485c]Floor and Mod(数论)

题意很简单,就是要求满足a=kb+k(0<k<b)的(a,b)的数量,且1≤a≤x,1≤b≤y。目测是个很可搞的数学题,那么我们来胡乱分析一波。

假如没有k<b的限制。对于处于1到y范围内的每个b,它的贡献就是x/(b+1)。这个是可以用整除分块做的。

以防有人不知道(没有

什么是整除分块呢?就是因为从1到y枚举b肯定是不现实的,但我们可以注意到,因为这里的除是要向下取整,所以其实对于一个区间内的b,x/(b+1)的值应该都是一样的。那么我们只需要对于每个区间左端点计算出x/(b+1),再算出区间长度,相乘即为该区间的贡献,然后直接去计算下一个区间。而关于如何找区间端点,举几个数试一试也可以发现,每个区间的右端点显然是可以被整除的点。实现可以看下面的代码,由于此题有y的限制,而且是x/(b+1),而非x/b,所以跟纯整除分块板子略有区别。

回过头考虑本题还有k<b的限制。所以每个b的贡献其实是min(x/(b+1),b-1)。显然,b大到一定程度之后就都是取x/(b+1)了,只要找到分界点,分界点以前用求和公式,以后整除分块即可。分界点大概在b=sqrt(x+1)附近,记得不能大于y的限制。

#include<cstdio>
#include<cmath>
#include<algorithm>
#define ll long long
using namespace std;
const int N=100010;
int t,x,y,l,r;
ll ans;
inline int read(){
	int x=0,f=0;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=1;ch=getchar();}
	while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
	return f?-x:x;
}
inline int min1(int i,int j){return i<j?i:j;}
int main(){
	t=read();
	while(t--){
		x=read(),y=read();
		l=min1((int)sqrt(x-1),y);
		while(l<y&&l<x/(l+2))++l;
		ans=1ll*l*(l-1)/2;
		for(l=l+2,r;l<=x&&l-1<=y;l=r+1){
			r=min1(x/(x/l),y+1);
			ans+=1ll*(x/l)*(r-l+1);
		}
		printf("%lld\n",ans);
	}
	return 0;
}

 

  • 2
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值