[BZOJ2154]Crash的数字表格莫比乌斯反演

orz PoPoQQQ 课件上的例题啊orzorz

话说这种根号划分的方法好像次次都有的样子orzorz

http://wenku.baidu.com/link?url=RRtdDApIUqzKmUDS4IOPU0MQnOLMJ6xXv0hFSajFA05YkoiGQdrN6koBLNMrmxnMs0MXXRL_ELxQfwSOvVdtSsNLaKao96PtBY63tJai903

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#include<queue>
#define SF scanf
#define PF printf
using namespace std;
typedef long long LL;
const int MAXN = 10000000;
const int MOD = 20101009;
int mu[MAXN+10], prime[MAXN+10], S[MAXN+10];
int n, m, tot;
bool vis[MAXN+10];
void init() {mu[1] = 1;for(int i = 2; i <= n; i++) {if(!vis[i]) prime[++tot] = i, mu[i] = -1;for(int j = 1; 1LL * prime[j] * i <= n; j++) {vis[prime[j]*i] = true;if(i % prime[j] == 0) {mu[prime[j]*i] = 0;break;}mu[prime[j]*i] = -mu[i];}}for(int i = 1; i <= n; i++) S[i] = (S[i-1] + 1LL * i * i * mu[i] % MOD) % MOD;
}
LL sum(LL x, LL y) {return (x * (x+1) / 2) % MOD * (y * (y+1) / 2 % MOD) % MOD;
}
LL calc(int x, int y) {LL ret = 0, Next_pos;if(x > y) swap(x, y);for(int i = 1; i <= x; i = Next_pos+1) {Next_pos = min(x / (x/i), y / (y/i));ret = (ret + (S[Next_pos] - S[i-1]) * sum(x/i, y/i) % MOD) % MOD;}return ret;
}
int main() {SF("%d%d", &n, &m);if(n > m) swap(n, m);init();LL Next_pos, ans = 0;for(int i = 1; i <= n; i = Next_pos+1) {Next_pos = min(n / (n/i), m / (m/i));ans = (ans + (i+Next_pos) * (Next_pos-i+1) / 2 % MOD * calc(n/i, m/i) % MOD) % MOD;}cout << (ans+MOD)%MOD;
}

[BZOJ2154]Crash的数字表格莫比乌斯反演相关推荐

【BZOJ2154】Crash的数字表格 [莫比乌斯反演]
Crash的数字表格 Time Limit: 20 Sec Memory Limit: 259 MB [Submit][Status][Discuss] Description 今天的数学课上,Cr ...
BZOJ 2154 Crash的数字表格 (莫比乌斯反演)
Crash的数字表格今天的数学课上,Crash小朋友学习了最小公倍数(Least Common Multiple).对于两个正整数a和b,LCM(a, b)表示能同时被a和b整除的最小正整数.例如, ...
【bzoj2154】Crash的数字表格莫比乌斯反演
题目描述今天的数学课上,Crash小朋友学习了最小公倍数(Least Common Multiple).对于两个正整数a和b,LCM(a, b)表示能同时被a和b整除的最小正整数.例如,LCM(6, ...
洛谷 - P1829 - Crash的数字表格 - 莫比乌斯反演
求: $S(n,m)=\sum\limits_{i=1}^{n}\sum\limits_{j=1}^{m}lcm(i,j)$ 显然: \(S(n,m)=\sum\limits_{i=1}^{n}\ ...
bzoj2154: Crash的数字表格
好神的莫比乌斯函数然后O(sqrt(n)*sqrt(n))好神的优化啊. #include<cstdio> #include<cstring> #include<ccty ...
bzoj2154 Crash的数字表格
2154: Crash的数字表格 Time Limit: 20 Sec Memory Limit: 259 MB Submit: 4549 Solved: 1643 [Submit][Status ...
BZOJ2154: Crash的数字表格 BZOJ2693: jzptab
[传送门:BZOJ2154&BZOJ2693] 简要题意: 给出n,m,求$\sum_{i=1}^{n}\sum_{j=1}^{m}LCM(i,j)$ 题解: 莫比乌斯反演(因为BZOJ269 ...
BZOJ 2154 [国家集训队]Crash的数字表格 / JZPTAB（莫比乌斯反演，经典好题）（Luogu P1829）
整理的算法模板合集: ACM模板点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 P1829 [国家集训队]Crash的数字表格 / JZPTAB(反演,经典好题) Problem S ...
P1829 [国家集训队]Crash的数字表格 / JZPTAB（莫比乌斯反演）
[国家集训队]Crash的数字表格 / JZPTAB 题目描述今天的数学课上,Crash 小朋友学习了最小公倍数(Least Common Multiple).对于两个正整数 a a a 和 b b ...

[BZOJ2154]Crash的数字表格莫比乌斯反演

[BZOJ2154]Crash的数字表格莫比乌斯反演相关推荐

最新文章

热门文章

[BZOJ2154]Crash的数字表格 莫比乌斯反演

[BZOJ2154]Crash的数字表格 莫比乌斯反演相关推荐

最新文章

热门文章

[BZOJ2154]Crash的数字表格莫比乌斯反演

[BZOJ2154]Crash的数字表格莫比乌斯反演相关推荐