【bzoj 1257】[CQOI2007]余数之和sum-CSDN博客

本文链接：https://blog.csdn.net/reverie_mjp/article/details/53202468

1257: [CQOI2007]余数之和sum

Time Limit: 5 Sec Memory Limit: 162 Mb Submit: 3708 Solved: 1704
[Submit][Status][Discuss]

Description

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

Input

输入仅一行，包含两个整数n, k。

Output

输出仅一行，即j(n, k)。

Sample Input

5 3

Sample Output

HINT

50%的数据满足：1<=n, k<=1000
100%的数据满足：1<=n ,k<=10^9

Source
[Submit][Status][Discuss]

【题解】

\sum i = 1 n k m o d i

$\sum_{i=1}^n k mod i$

= \sum i = 1 n k - ⌊ k i ⌋ * i

$=\sum_{i=1}^n k-\lfloor\frac{k}{i}\rfloor *i$

= n * k - \sum i = 1 n ⌊ k i ⌋ * i

$=n*k-\sum_{i=1}^n \lfloor\frac{k}{i}\rfloor*i$

当 i > k 时 ， ⌊ k i ⌋ * i = 0

$当i>k时，\lfloor\frac{k}{i}\rfloor*i =0$
so,只需考虑

i≤k $i≤k$ 的情况，即：

= n * k - \sum i = 1 m i n (n, k) ⌊ k i ⌋ * i

$=n*k-\sum_{i=1}^{min(n,k)}\lfloor\frac{k}{i}\rfloor *i$

根 据 性 质 可 知 ： ⌊ k i ⌋ 的 取 值 不 超 过 k \sqrt 个 连 续 区 间

$根据性质可知：\lfloor\frac{k}{i}\rfloor的取值不超过\sqrt{k}个连续区间$

用 i 从 1 到 m i n (n, k) 枚 举 ， 每 次 找 到 取 值 w ， 算 出 左 右 区 间 。 左 区 间 ： l = i ， 右 区 间 ： 当 ⌊ k i ⌋ = w 时 ， 根 据 高 斯 消 元 的 性 质 ， w = ⌊ k i ⌋ \leq k i

$用i从1到min(n,k)枚举，每次找到取值w，算出左右区间。左区间：l=i，右区间：当\lfloor\frac{k}{i}\rfloor = w 时，根据高斯消元的性质，w=\lfloor\frac{k}{i}\rfloor≤\frac{k}{i}$

s o, i \geq k w, r = k w

$so,i≥\frac{k}{w},r=\frac{k}{w}$

由 于 可 能 存 在 r \geq n ， r = m i n (r ， w)

$由于可能存在r≥n，r=min(r，w)$

对 于 每 个 区 间 ， a n s + = w * \sum q = l r q ， 然 后 到 下 一 区 间 ： i = r + 1

$对于每个区间，ans+=w*\sum_{q=l}^r q，然后到下一区间：i=r+1$

结 果 为 a n s

$结果为ans$

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
long long n,m,ans;
inline long long js(long long x)
{
    return x*(x+1)>>1;
}
int main()
{
    long long i,j;
    scanf("%I64d%I64d",&n,&m);
    for (i=1,j=0;i<=min(n,m);i=j+1)
     {
        j=m/(m/i);//计算一段数的结尾 
        if (j>n) j=n;//如果j跑到n之外，就把它拽回n 
        ans+=m/i*((js(j)-js(i-1)));//用前缀和来维护每一段的和 
     }
    ans=n*m-ans; 
    printf("%I64d\n",ans); 
    return 0;
 }