BZOJ1257 洛谷2261 [CQOI2007]余数求和

最新推荐文章于 2019-12-29 21:30:52 发布

WT_cnyali

最新推荐文章于 2019-12-29 21:30:52 发布

阅读量497

点赞数

分类专栏：数学分块文章标签：数学分块

本文链接：https://blog.csdn.net/WT_cnyali/article/details/78369526

版权

数学同时被 2 个专栏收录

9 篇文章 0 订阅

订阅专栏

分块

2 篇文章 0 订阅

订阅专栏

Description
给定 $n$ , $x$ ，求

$\sum i = 1 n x m o d i$ $\sum_{i=1}^{n}x \ mod \ i$
Solution

$= = \sum i = 1 n x m o d i \sum i = 1 n x - ⌊ x i ⌋ \times i n x - \sum i = 1 n ⌊ x i ⌋ \times i$ $\begin{eqnarray} &&\sum_{i=1}^{n}x \ mod \ i \\ &=& \sum_{i=1}^{n}{x - \lfloor \dfrac{x}{i} \rfloor \times i} \\ &=& nx - \sum_{i=1}^{n}{\lfloor \dfrac{x}{i} \rfloor \times i} \end{eqnarray}$
后面那一堆显然有大段大段都是一样的，最多不同的个数是 $n$ 的约数个数。时间复杂度 $O(\sqrt{n})$ 。

#include <bits/stdc++.h>
using namespace std ;
#define LL long long
LL ans, n, x ;
int main() {
#ifndef ONLINE_JUDGE
    freopen ( "BZOJ1257.in", "r", stdin ) ;
    freopen ( "BZOJ1257.out", "w", stdout ) ;
#endif
    LL i, j, k, lim, l, r ;
    scanf ( "%lld%lld", &n, &x ) ;
    ans = n*x ;
    lim = sqrt(x) ;
    for ( i = 1 ; i <= n ; ) {
        k = x/i ;
        if (!k) break ;
        l = i ;
        if ( (x+1)%k == 0 ) r = (x+1)/k-1 ;
        else r = floor( (x+1)*1.0/k ) ;
        r = min(r, n) ;
        ans -= (l+r)*(r-l+1)/2*k ;
        i = r+1 ;
    }
    printf ( "%lld\n", ans ) ;  
    return 0 ;
}