[bzoj] 1257 余数之和sum || 数论

原题

给出正整数n和k,计算j(n, k)=k mod 1 + k mod 2 + k mod 3 + … + k mod n的值,其中k mod i表示k除以i的余数。


\(\sum^n_{i=1}k\%i\)
\(=\sum^n_{i=1}k-\lfloor k/i \rfloor*i\)
\(=n*k-\sum^n_{i=1}\lfloor k/i \rfloor*i\)

\(\lfloor k/i \rfloor\)只有\(\sqrt k\)个取值
证明:
对于所有\(>\sqrt k\)的数,\(\lfloor k/i \rfloor\)一定是一个对应的\(<\sqrt k\)的值,所以最多只有\(2\sqrt k\)个值
也就是说\(\lfloor k/i \rfloor\)的取值是这样的:
1279575-20180105085554424-1400124160.png

所以每次i为左端点,k/(k/i)为右端点,这一段就可以直接处理。复杂度为\(O(\sqrt n)\)

#include<cstdio>
typedef long long ll;
using namespace std;
int n,k;
ll ans;

int main()
{
    scanf("%d%d",&n,&k);
    if (n>k) ans=(ll)(n-k)*k,n=k;
    int r;
    for (int i=1;i<=n;i=r+1)
    {
    int t=k/i;r=k/t;
    if (r>=n) r=n;
    ans+=(ll)(r-i+1)*k-(ll)(r-i+1)*(i+r)/2*t;
    }
    printf("%lld\n",ans);
    return 0;
}

转载于:https://www.cnblogs.com/mrha/p/8202475.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值