【BZOJ1257】余数之和(数论分块,暴力)

【BZOJ1257】余数之和(数论分块,暴力)

题解

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

7

HINT

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

题解

很Interesting的一道题目
所求为
\[\sum_{i=1}^n{k}\ \%\ {i}\]
而这个式子可以再改一下
\[\sum_{i=1}^nk-i*(k\ div \ i)\]
其中\(div\)是整除
前面的\(n*k\)直接先算出来
后面的东西,很容易观察到,\(i\)是单增的
而\(k/i\)在一段范围内是不会变化的
因此,每次求出这一段范围,然后直接计算等差数列即可

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
inline int read()
{
    int x=0,t=1;char ch=getchar();
    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
    if(ch=='-')t=-1,ch=getchar();
    while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
    return x*t;
}
long long ans,n,k;
int main()
{
    n=read();k=read();
    ans=n*k;
    for(long long l=1,r=0;l<=n;l=r+1)
    {
        if(k/l)r=min(n,k/(k/l));
        else r=n;
        ans-=(k/l)*(r-l+1)*(l+r)>>1;
    }

    printf("%lld\n",ans);
    return 0;
}
时间: 2024-10-13 20:07:54

【BZOJ1257】余数之和(数论分块,暴力)的相关文章

bzoj 1257 [CQOI2007]余数之和——数论分块

题目:https://www.lydsy.com/JudgeOnline/problem.php?id=1257 \( n\%i = n - \left \lfloor n/i \right \rfloor * i \) 注意 n<k 时当前块的右端点可能超过 n ! #include<cstdio> #include<cstring> #include<algorithm> #define ll long long using namespace std; in

余数(数论分块)

题目描述: 题解: 首先容易想到:当 i > n 时,n mod i = n 所以如果 m > n , ans+=((m-n)%mod)*(n%mod)%mod; m=n; 接下来考虑 i<=n 的情况: 后面这个东西很明显可以用数论分块算,时间复杂度是O(sqrt(N)) 但我当时做的时候还没有学数论分块,所以代码十分丑陋 #include<iostream> #include<cstdio> #include<cmath> using namesp

bzoj 1257 : [CQOI2007]余数之和 (数学+分块)

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. 1<=n ,k<=10^9 Output 输出仅一行,即j(n, k). Sample Input 5 3 Sample

ACM学习历程—SNNUOJ1132 余数之和(数论)

Description F(n) = (n % 1) + (n % 2) + (n % 3) + ...... (n % n).其中%表示Mod,也就是余数.例如F(6) = 6 % 1 + 6 % 2 + 6 % 3 + 6 % 4 + 6 % 5 + 6 % 6 = 0 + 0 + 0 + 2 + 1 + 0 = 3.给出n,计算F(n). Input 输入1个数N(2 <= N <= 10^12). Output 输出F(n). Sample Input 6 Sample Output

[日常摸鱼]bzoj1257余数之和

题意:输入$k,n$,求$\sum_{i=1}^n k \mod i$ $k \mod i=k-i*\lfloor \frac{k}{i} \rfloor $,$n$个$k$直接求和,后面那个东西像比较套路的分段求和 算k/(k/i)这种东西的时候还要注意判一下分母为0什么的- #include<cstdio> typedef long long lint; lint n,k,ans; inline lint min(lint a,lint b){return a<b?a:b;} inl

BZOJ 1257: [CQOI2007]余数之和sum( 数论 )

n >= k 部分对答案的贡献为 k * (n - k) n < k 部分贡献为 ∑ (k - ⌊k / i⌋ * i)  = ∑  , ⌊k / i⌋ 相等的数是连续的一段, 此时这段连续的数对答案的贡献成等差数列, 可以O(1)求出..然后就分⌊k / i⌋ 相等的一块一块来就行了. 分出来大概是sqrt(k)块.这个sqrt(k)我并不会证Orz...写了个程序验证了一下, 分出来的块数和2 * sqrt(k)非常接近. 所以时间复杂度为O(sqrt(k)) --------------

bzoj千题计划173:bzoj1257: [CQOI2007]余数之和sum

http://www.lydsy.com/JudgeOnline/problem.php?id=1257 k%i=k-int(k/i)*i 除法分块,对于相同的k/i用等差序列求和来做 #include<cstdio> #include<iostream> using namespace std; int main() { int n,k; scanf("%d%d",&n,&k); long long ans=0; if(n>k) { an

数论分块之整除分块

前言 最近在学莫比乌斯反演,然而只看懂了莫比乌斯函数,然后反演看着一脸懵逼,最后只看懂了数论分块里面的一个分支内容(也是莫比乌斯反演的前置姿势),整除分块 于是写一篇博文记录一下整除分块 整除分块 整除分块是用于快速处理形似 \[ \sum_{i=1}^{n}{\lfloor \frac{n}{i} \rfloor} \] 的式子的方法 很显然,这个可以\(O(n)\)得到答案.但是,在某些题目中,毒瘤出题人将数据加强到了\(10^{10}\)以上,这个时候我们就无法通过\(O(n)\)的解法来

余数之和sum

余数之和sum Time Limit:5000MS   Memory Limit:165888KB   64bit IO Format:%lld & %llu SubmitStatus 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