1584 加权约数和

题意:求\(\sum_{i=1}^{N} \sum_{j=1}^{N} {\max(i,j)\cdot \sigma(i\cdot j)}\)

多组数据\(n \le 10^6, T \le 50000\)


这道题有两步我感到非常神奇。tls好强啊。


首先,怎么处理\(max(i,j)\)

\[max(i,j) = \sum_{k=1}^n[k\le i \ or\ k \le j] = n-\sum_{k=1}^n[k>i][k>j]
\]

这样转化之后再代入,可以得到

\[ans = n*\sum_{i=1}^n\sum_{j=1}^n \sigma(ij)-\sum_{k=1}^n\sum_{i=1}^{k-1} \sum_{j=1}^{k-1}\sigma(ij)= n*f(n) - s(n) \\
\]

之前的题目已经推倒过\(f(n)\)

\[f(n) = \sum_{d=1}^n \mu(d) \cdot d \cdot (\sum_{i=1}^{\frac{n}{d}} \sigma_1(i))^2 \\
令g(n)= \mu(d)\cdot d,\ h(n) = (\sum_{i=1}^{n} \sigma_1(i))^2 \\
f(n) = \sum_{i=1}^n g(i) h(\frac{n}{i})
\]

因为还需要f的前缀和,整除分块计算的话复杂度\(O(n\sqrt{n})\)承受不了...


第二步很神的做法

对于一个\(i\),\(\frac{x}{i}\)的取值对于一段\(x\)也是相同的!就是每段\([k*i,k*i+i-1]\)

我们枚举\(i\),然后枚举取值相同的\(x\)的段,给这些\(f(x)\)加上\(g(i) h(\frac{x}{i})\)

对f进行差分可以做到\(O(1)\)区间加,那么根据调和级数求和复杂度\(O(nlogn)\)


对于单个函数需要整除分块的计算的,求他们的前缀和,都可以使用这个技巧优化到\(O(nlogn)\)


我的常数有点大,改了几个long long就好了,看来类型转化很慢啊

#include <iostream>
#include <stdio.h>
#include <cstring>
#include <algorithm>
#include <cmath>
//using namespace std;
typedef long long ll;
const int N=1e6+5, mo=1e9+7;
int U=1e6;
inline int read(){
char c=getchar(); int x=0,f=1;
while(c<'0' || c>'9') {if(c=='-')f=-1; c=getchar();}
while(c>='0' && c<='9') {x=x*10+c-'0'; c=getchar();}
return x*f;
} inline void mod(int &x) {if(x>=mo) x-=mo; else if(x<0) x+=mo;}
inline void mod(ll &x) {if(x>=mo) x-=mo; else if(x<0) x+=mo;}
bool notp[N]; int p[N/10], lp[N], mu[N]; ll si[N], g[N], h[N];
void sieve(int n) {
mu[1] = 1; si[1] = 1;
for(int i=2; i<=n; i++) {
if(!notp[i]) p[++p[0]] = i, mu[i] = -1, lp[i] = i, si[i] = 1+i;
for(int j=1; j <= p[0] && i*p[j] <= n; j++) {
int t = i*p[j]; notp[t] = 1;
if(i % p[j] == 0) {
mu[t] = 0;
lp[t] = lp[i] * p[j];
if(t == lp[t]) mod(si[t] = si[i] + lp[t]);
else si[t] = si[t / lp[t]] * si[lp[t]] %mo;
break;
}
lp[t] = p[j];
mu[t] = -mu[i];
si[t] = si[i] * (p[j] + 1) %mo;
}
}
for(int i=1; i<=n; i++) mod(g[i] = mu[i] * i), mod(si[i] += si[i-1]), h[i] = si[i] * si[i] %mo;
} int f[N], s[N];
void init(int n) {
for(int i=1; i<=n; i++) {
for(int x=i, k=1; x<=n; x+=i, k++) {
int l = x, r = x+i-1, d = g[i] * h[k] %mo;
if(r > n) r = n;
mod(f[l] += d); mod(f[r+1] -= d);
}
}
for(int i=1; i<=n; i++) mod(f[i] += f[i-1]), mod(s[i] = s[i-1] + f[i]);
} int n;
int main() {
freopen("in", "r", stdin);
sieve(U);
init(U);
int T=read();
for(int i=1; i<=T; i++) {
n=read();
int ans = (ll) f[n] * n %mo - s[n-1]; mod(ans);
printf("Case #%d: %d\n", i, ans);
}
}

51NOD 1584 加权约数和 [莫比乌斯反演 转化 Trick]的更多相关文章

  1. 51nod 1584 加权约数和 约数和函数小trick 莫比乌斯反演

    LINK:加权约数和 我曾经一度认为莫比乌斯反演都是板子题. 做过这道题我认输了 不是什么东西都是板子. 一个trick 设\(s(x)\)为x的约数和函数. 有 \(s(i\cdot j)=\sum ...

  2. [51Nod 1584] 加权约数和

    Description 在整理以前的试题时,他发现了这样一道题目:"求 \(\sum\sigma(i)\),其中 \(1≤i≤N\),\(σ(i)\) 表示 \(i\) 的约数之和.&quo ...

  3. 51nod 1584加权约数和

    学到了好多东西啊这题... https://blog.csdn.net/sdfzyhx/article/details/72968468 #include<bits/stdc++.h> u ...

  4. 【51nod】1222 最小公倍数计数 莫比乌斯反演+组合计数

    [题意]给定a和b,求满足a<=lcm(x,y)<=b && x<y的数对(x,y)个数.a,b<=10^11. [算法]莫比乌斯反演+组合计数 [题解]★具体 ...

  5. BZOJ 3994: [SDOI2015]约数个数和 [莫比乌斯反演 转化]

    2015 题意:\(d(i)\)为i的约数个数,求\(\sum\limits_{i=1}^n \sum\limits_{j=1}^m d(ij)\) \(ij\)都爆int了.... 一开始想容斥一下 ...

  6. 51nod 1222 最小公倍数计数【莫比乌斯反演】

    参考:https://www.cnblogs.com/SilverNebula/p/7045199.html 所是反演其实反演作用不大,又是一道做起来感觉诡异的题 转成前缀和相减的形式 \[ \sum ...

  7. 51nod 1594 Gcd and Phi(莫比乌斯反演)

    题目链接 传送门 思路 如果这题是这样的: \[ F(n)=\sum\limits_{i=1}^{n}\sum\limits_{j=1}^{n}\phi(gcd(i,j)) \] 那么我们可能会想到下 ...

  8. [51 Nod 1584] 加权约数和

    题意 求∑i=1N∑j=1Nmax(i,j)⋅σ1(ij)\large \sum_{i=1}^N\sum_{j=1}^Nmax(i,j)\cdot\sigma_1(ij)i=1∑N​j=1∑N​max ...

  9. YbtOJ#943-平方约数【莫比乌斯反演,平衡规划】

    正题 题目链接:http://www.ybtoj.com.cn/contest/122/problem/3 题目大意 \(S(i)\)表示\(i\)的约数个数,\(Q\)次询问给出\(n,m\)求 \ ...

随机推荐

  1. Educational Codeforces Round 2_B. Queries about less or equal elements

    B. Queries about less or equal elements time limit per test 2 seconds memory limit per test 256 mega ...

  2. Generator函数异步应用

    转载请注明出处: Generator函数异步应用 上一篇文章详细的介绍了Generator函数的语法,这篇文章来说一下如何使用Generator函数来实现异步编程. 或许用Generator函数来实现 ...

  3. java如何获取一个对象的大小

    When---什么时候需要知道对象的内存大小 在内存足够用的情况下我们是不需要考虑java中一个对象所占内存大小的.但当一个系统的内存有限,或者某块程序代码允许使用的内存大小有限制,又或者设计一个缓存 ...

  4. Ajax beforeSend和complete 方法

    http://blog.csdn.net/chenjianandiyi/article/details/52274591 .ajax({ beforeSend: function(){ // Hand ...

  5. LAMP和LNMP,你更愿意选择谁,为什么?

    https://www.zhihu.com/question/19697826 http://www.simongong.net lamp 的全称是linux + apache + mysql +ph ...

  6. js立体旋转展示效果

    <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/ ...

  7. IOS学习笔记27—使用GDataXML解析XML文档

    http://blog.csdn.net/ryantang03/article/details/7868246

  8. MySQL 查询最近几天的记录 最近7天的记录 本周内的记录

    本周内:select * from wap_content where week(created_at) = week(now) 查询一天:select * from table where to_d ...

  9. Jmeter下载时Binaries和Source两类包的区别

    在下载Jmeter时,存在两种类型的下载包,分别为Binaries和Source: 一般开放原代码软件都会有两个版本发布: Source Distribution 和 Binary Distribut ...

  10. linux_samba服务安装

    什么是samba服务? 用于Windows和linux系统之间实现共享文件的目的服务 如何配置其服务? Linux端: 搭建服务 1. 安装samba yum install -y samba 2. ...