Evanyou Blog 彩带

洛谷P2261 [CQOI2007] 余数求和 [数论分块]

  题目传送门

余数求和

题目背景

数学题,无背景

题目描述

给出正整数n和k,计算G(n, k)=k mod 1 + k mod 2 + k mod 3 + … + k mod n的值,其中k mod i表示k除以i的余数。例如G(10, 5)=5 mod 1 + 5 mod 2 + 5 mod 3 + 5 mod 4 + 5 mod 5 …… + 5 mod 10=0+1+2+1+0+5+5+5+5+5=29

输入输出格式

输入格式:

 

两个整数n k

 

输出格式:

 

答案

 

输入输出样例

输入样例#1: 
10 5
输出样例#1: 
29

说明

30%: n,k <= 1000

60%: n,k <= 10^6

100% n,k <= 10^9


  分析:

  之前没怎么写过数论分块(蒟蒻并不会莫比乌斯反演),于是做道题练下手。

  因为$a\%b=a-b*\lfloor \frac{a}{b}\rfloor$,所以我们所求的式子$\sum^n_{i=1}k\mod i$可以转化为$n*k-\sum^n_i{i*\lfloor\frac{k}{i}\rfloor}$。

  和式的部分就可以用整除分块来做,复杂度就是$O(\sqrt{n})$的。

  Code:

 

//It is made by HolseLee on 7th Nov 2018
//Luogu.org P2261
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;

typedef long long ll;
ll n,k,ans;

int main()
{
    cin>>n>>k;
    ans=n*k;
    for(ll l=1,r; l<=n; l=r+1) {
        if( l<=k ) r=min(k/(k/l),n);
        else r=n;
        ans-=(r-l+1)*(r+l)*(k/l)/2;
    }
    cout<<ans<<'\n';
    return 0;
}

 

posted @ 2018-11-07 20:21  HolseLee  阅读(210)  评论(0编辑  收藏  举报