2020牛客多校第七场H-Dividing

Dividing

题意

(1,k)为一个元组
如果(n,k)是元组的话,那么(n+k,k)和(nk,k)也是元组

思路

由于元组性质,如果(n,k)为元组,那么一定能够把n除k或减k到1。
故n=n’k+1或者n=n’k
即n%k=0或者n%k=1。
即为计算
∑ i = 1 k ⌊ N i ⌋ + ∑ i = 2 k ⌊ N − 1 i ⌋ + k − 1 \sum_{i=1}^{k}\left \lfloor \frac{N}{i} \right \rfloor+\sum_{i=2}^{k}\left \lfloor \frac{N-1}{i} \right \rfloor+k-1 i=1kiN+i=2kiN1+k1
注意:
因为k=1时N%k都为0,所以会被计算在第一项中,第二项就不需要再枚举k=1了。
被第一项和第二项计算到的元组的n等于n’k+1或者n’k,n’>=1,所以并未计算n=1的元组,由于k=1特殊,已经被第一项包括了,所以答案需要再加上k-1。

关于分块O( n \sqrt{n} n )计算 ∑ i = 1 k ⌊ N i ⌋ \sum_{i=1}^{k}\left \lfloor \frac{N}{i} \right \rfloor i=1kiN,参考大佬博客
里面有结论。

代码

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=1e9+7;

int main(){
    ll n,k,ans=0;
    cin>>n>>k;
    for(ll l=1,r,p;l<=min(n,k);l=r+1){
        p=n/l;
        r=min(k,n/p);
        ans=(ans+(r-l+1)%mod*p%mod)%mod;
    }
    n--;
    for(ll l=2,r,p;l<=min(n,k);l=r+1){
        p=n/l;
        r=min(k,n/p);
        ans=(ans+(r-l+1)%mod*p%mod)%mod;
    }
    ans=(ans+k-1)%mod;
    printf("%lld\n",ans);
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值