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. debian右键添加在终端中打开

    sudo apt-get install nautilus-open-terminal -y 注销,重启

  2. Linux学习-Tarball 的管理与建议

    使用原始码管理软件所需要的基础软件 从原始码的说明我们晓得要制作一个 binary program 需要很多咚咚的呢!这包括底下这些基础的软件: gcc 或 cc 等 C 语言编译程序 (compil ...

  3. Python IO多路复用select模块

    多路复用的分析实例:服务端.客户端 #服务端配置 from socket import * import time import select server = socket(AF_INET, SOC ...

  4. uiautomator 一个简单脚本创建流程

    http://www.codeceo.com/article/android-ui-auto-test.html

  5. 2018 ACM 国际大学生程序设计竞赛上海大都会赛重现赛 F Color it

    链接:https://www.nowcoder.com/acm/contest/163/F 来源:牛客网 2018 ACM 国际大学生程序设计竞赛上海大都会赛重现赛 F Color it 时间限制:C ...

  6. PHP 格式化显示时间 date()函数

    Y 4位数字年,y为2位数字,如99即1999年 m 数字月份,前面有前导0,如01.n 为无前导0数字月份 F 月份,完整的文本格式,例如 January 或者 March M 三个字母缩写表示的月 ...

  7. Jeddict1:简述

    备注:本文是一篇翻译文,文章所有的内容(排除目录)均来自于:https://jeddict.github.io/ ,但为了自己的理解,进行了目录分类! 一.什么是Jeddict jeddict是一个开 ...

  8. UVa——540Team Queue(STL练习map、queue数组的综合使用)

    Team Queue Time Limit:                                                        3000MS                 ...

  9. el-select绑定值为对象时,报错[Vue warn]: <transition-group> children must be keyed: <ElTag>

    解决方法: <el-select v-model="syncParams.toSlaveList" multiple value-key="ip" pla ...

  10. docker 查询或获取私有仓库(registry)中的镜像

    docker 查询或获取私有仓库(registry)中的镜像,使用 docker search 192.168.1.8:5000 命令经测试不好使. 解决: 1.获取仓库类的镜像: [root@sha ...