1220 约数之和 
                                题目来源: Project Euler
基准时间限制:3 秒 空间限制:131072 KB 分值: 640 难度:8级算法题
 
Discription
d(k)表示k的所有约数的和。d(6) = 1 + 2 + 3 + 6 = 12。
定义S(N) = ∑1<=i<=N ∑1<=j<=N d(i*j)。
例如:S(3) = d(1) + d(2) + d(3) + d(2) + d(4) + d(6) + d(3) + d(6) + d(9) = 59,S(1000) = 563576517282。
给出正整数N,求S(N),由于结果可能会很大,输出Mod 1000000007(10^9 + 7)的结果。
 
Input
输入一个数N(2 <= N <= 10^9)。
Output
输出S(N) Mod 1000000007(10^9 + 7)的结果。
Input示例
1000
Output示例
576513341

我们知道的是,当d(x)表示x的约数的时候,d(i*j)=Σ(p|i)Σ(q|j) [gcd(p,q)==1]
但是当d(x)表示x的约数之和的时候,d(i*j)=Σ(p|i)Σ(q|j) p*q [gcd(p,j/q)==1]
这两者都可以通过质因子分解来证明。
然后式子就好推了,我这里就不推了23333
#include<bits/stdc++.h>
#define ll long long
#define ha 1000000000
using namespace std;
const int maxn=10000000;
ll zs[maxn/5],miu[maxn+5];
ll t=0,low[maxn+5];
bool v[maxn+5];
ll d[maxn+5];
map<ll,ll> mmpd;
map<ll,ll> mmpmiu;
ll n; inline ll add(ll x,ll y){
x+=y;
return x>=ha?x-ha:x;
} inline ll c(ll x){
if(x>ha) x%=ha;
return (x*(x+1)>>1)%ha;
} inline void init(){
d[1]=1,low[1]=1,miu[1]=1;
for(int i=2;i<=maxn;i++){
if(!v[i]) zs[++t]=i,miu[i]=-1,d[i]=i+1,low[i]=i;
for(int j=1,u;j<=t&&(u=zs[j]*i)<=maxn;j++){
v[u]=1;
if(!(i%zs[j])){
low[u]=low[i]*zs[j];
if(low[i]==i) d[u]=d[i]*zs[j]+1ll;
else d[u]=d[low[u]]*d[i/low[i]];
break;
}
low[u]=zs[j];
miu[u]=-miu[i];
d[u]=d[i]*(zs[j]+1);
}
} for(int i=1;i<=maxn;i++){
d[i]=add(d[i-1],d[i]);
miu[i]=add(add(miu[i]*i,ha),miu[i-1]);
}
} inline ll getmiu(ll x){
if(x<=maxn) return miu[x];
if(mmpmiu.count(x)) return mmpmiu[x]; ll an=ha-1;
for(ll i=2,j,now;i<=x;i=j+1){
now=x/i,j=x/now;
an=add(an,add(c(j),ha-c(i-1))*(ll)getmiu(now)%ha);
}
an=ha-an;
mmpmiu[x]=an;
return an;
} inline ll getd(ll x){
if(x<=maxn) return d[x];
if(mmpd.count(x)) return mmpd[x]; ll an=0;
for(ll i=1,j,now;i<=x;i=j+1){
now=x/i,j=x/now;
an=add(an,add(c(j),ha-c(i-1))*(now%ha)%ha);
} mmpd[x]=an;
return an;
} inline void solve(){
ll pre=0,an=0,oops,val;
for(ll i=1,j,now;i<=n;i=j+1){
now=n/i,j=n/now,oops=getmiu(j);
val=getd(now),val=val*(ll)val%ha;
an=add(an,add(oops,ha-pre)*(ll)val%ha);
pre=oops;
} printf("%lld\n",an);
} int main(){
init();
scanf("%lld",&n);
solve(); return 0;
}

  

 

  

 

51Nod 约数之和的更多相关文章

  1. 51NOD 1220 约数之和 [杜教筛]

    1220 约数之和 题意:求\(\sum_{i=1}^n \sum_{j=1}^n \sigma_1(ij)​\) \[ \sigma_0(ij) = \sum_{x\mid i}\sum_{y\mi ...

  2. [51Nod 1220] - 约数之和 (杜教筛)

    题面 令d(n)d(n)d(n)表示nnn的约数之和求 ∑i=1n∑j=1nd(ij)\large\sum_{i=1}^n\sum_{j=1}^nd(ij)i=1∑n​j=1∑n​d(ij) 题目分析 ...

  3. 【动态规划】mr359-最大公约数之和

    [题目大意] 选取和不超过S的若干个不同的正整数,使得所有数的约数(不含它本身)之和最大. 输入一个正整数S. 输出最大的约数之和. 样例输入 Sample Input 11 样例输出 Sample ...

  4. 约数之和(POJ1845 Sumdiv)

    最近应老延的要求再刷<算法进阶指南>(不得不说这本书不错)...这道题花费了较长时间~(当然也因为我太弱了)所以就写个比较易懂的题解啦~ 原题链接:POJ1845 翻译版题目(其实是AcW ...

  5. POJ1845Sumdiv题解--约数之和

    题目链接 https://cn.vjudge.net/problem/POJ-1845 分析 \(POJ\)里的数学题总是这么妙啊 首先有一个结论就是\(A=\prod{ \ {p_i}^{c_i} ...

  6. 51nod 1220 约数之和【莫比乌斯反演+杜教筛】

    首先由这样一个式子:\( d(ij)=\sum_{p|i}\sum_{q|j}[gcd(p,q)==1]\frac{pj}{q} \)大概感性证明一下吧我不会证 然后开始推: \[ \sum_{i=1 ...

  7. 【51nod1220】约数之和

    题目 d(k)表示k的所有约数的和.d(6) = 1 + 2 + 3 + 6 = 12. 定义S(N) = ∑1<=i<=N ∑1<=j<=N d(i*j). 例如:S(3) ...

  8. 51nod 约数和(数论)

    题目链接: 约数和 基准时间限制:2 秒 空间限制:131072 KB 分值: 80 有三个下标从1到n的数组a.b.c. a数组初始全为0. b[i]=∑j|ia[j] c[i]=∑j|ib[j] ...

  9. 51nod1220 约数之和

    http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1220 $G(n)=\sum\limits_{i=1}^n\sum\lim ...

随机推荐

  1. 使用python3调用MyQR库生成动态二维码(附源代码)

    可生成普通二维码.带图片的艺术二维码(黑白与彩色).动态二维码(黑白与彩色). GitHub:https://github.com/sylnsfar/qrcode 中文版:https://github ...

  2. STM8S与IAR程序常用错误

    一.IAR中的重复定义问题 在自己写头文件时,要记得将常量定义在.c文件中,如果将常量定义在.h文件中,当在main.c或者其他地方包含该头文件时,会将头文件中的常量定义包 含到main.c中,同时, ...

  3. body标签相关

    03-body标签中相关标签   今日主要内容: 列表标签 <ul>.<ol>.<dl> 表格标签 <table> 表单标签 <fom> 一 ...

  4. LiveScript 流程控制、循环以及列表推导式

    The LiveScript Book     The LiveScript Book Generators and Yield 你可以在你的 LiveScript 代码中使用 Ecmascript ...

  5. Python3异常-AttributeError: module 'sys' has no attribute 'setdefaultencoding'

    基于python3.6.1版本,在一个.py文件中,加入这3行: import requests, re, sys reload(sys) sys.setdefaultencoding("u ...

  6. 聊聊、Nginx 安装启动

    首先说下安装 Nginx 的步骤: (1)window 下安装 进入 http://nginx.org/en/download.html 下载版本 Mainline version 或者 Stable ...

  7. [办公软件][2]screenToGif

    https://github.com/NickeManarin/ScreenToGif/wiki/help 下载: .Net Framework 4.6.1   https://www.microso ...

  8. C语言总结(2)

    1.函数printf(" ")可以输出双引号中任何固定不变的内容. 2.必须在程序前面加:预处理命令. 3.#include<stdio.h>后面不需要“:”,“:”. ...

  9. Thanks for your encourage!

    将近三个月的学习,我的努力换回了代表荣誉的小黄衫,这令我很开心啊...我想是不是要写点什么来表达自己的心情呢=,=  于是就有了以下文字ahhhhhh... 学习心得: (1)学习中总会有失败和成功, ...

  10. spring AOP详解一

    AOP是Aspect Oriented Programing的简称,面向切面编程.AOP适合于那些具有横切逻辑的应用:如性能监测,访问控制,事务管理以及日志记录.AOP将这些分散在各个业务逻辑中的代码 ...