BZOJ2693jzptab
Sol
增加了数据组数T<=10000
推到
\(ans=\sum_{d=1}^{N}d*\sum_{i=1}^{\lfloor\frac{N}{d}\rfloor}\mu(i)*i^2*\frac{(\lfloor\frac{N}{d*i}\rfloor + 1) * \lfloor\frac{N}{d*i}\rfloor}{2}*\frac{(\lfloor\frac{M}{d*i}\rfloor + 1) * \lfloor\frac{M}{d*i}\rfloor}{2}\)
\(设S(i)=\frac{(i+1)*i}{2},将d*i换成k\)
\(原式=\sum_{k=1}^{N}S(\lfloor\frac{N}{k}\rfloor)*S(\lfloor\frac{M}{k}\rfloor)*k*\sum_{i|k}i*\mu(i)\)
\(设f(n)=\sum_{i|n}i*\mu(i),它是个积性函数,可以线性筛\)
# include <bits/stdc++.h>
# define RG register
# define IL inline
# define Fill(a, b) memset(a, b, sizeof(a))
using namespace std;
typedef long long ll;
const int _(1e7 + 1), MOD(1e8 + 9);
IL ll Read(){
char c = '%'; ll x = 0, z = 1;
for(; c > '9' || c < '0'; c = getchar()) if(c == '-') z = -1;
for(; c >= '0' && c <= '9'; c = getchar()) x = x * 10 + c - '0';
return x * z;
}
int prime[_], f[_], num, s[_];
bool isprime[_];
IL void Prepare(){
isprime[1] = 1; s[1] = f[1] = 1;
for(RG int i = 2; i < _; ++i){
if(!isprime[i]) prime[++num] = i, f[i] = 1 - i;
for(RG int j = 1; j <= num && i * prime[j] < _; ++j){
isprime[i * prime[j]] = 1;
if(i % prime[j]) f[i * prime[j]] = 1LL * f[i] * f[prime[j]] % MOD;
else{ f[i * prime[j]] = f[i]; break; }
}
s[i] = (s[i - 1] + 1LL * i * f[i] % MOD) % MOD;
}
}
IL ll S(RG ll x){ return x * (x + 1) / 2 % MOD; }
int main(RG int argc, RG char *argv[]){
Prepare();
for(RG ll T = Read(), n, m; T; --T){
RG ll ans = 0; n = Read(); m = Read();
if(n > m) swap(n, m);
for(RG ll i = 1, j; i <= n; i = j + 1){
j = min(n / (n / i), m / (m / i));
(ans += 1LL * (s[j] - s[i - 1] + MOD) % MOD * S(n / i) % MOD * S(m / i) % MOD) %= MOD;
}
printf("%lld\n", (ans + MOD) % MOD);
}
return 0;
}
BZOJ2693jzptab的更多相关文章
随机推荐
- zabbix 网络模板自动发现端口时,过滤掉某些特定规则的端口,减少item的方法
1.需求描述 默认情况下Zabbix 模板 中网络接口自动发现接口时,会产生很多item,有时候会有我们不需要的一些接口,这时候需要过滤掉他们. 比如我有一台运行kvm的服 ...
- Ubuntu配置Nginx虚拟主机和支持ThinkPHP
[Nginx配置虚拟主机] 每一个 server { listen 80; server_name www.a.com; ..... } 就表示一台虚拟域名, 然后对应的 ...
- Spring Boot - Font Awesome OTS parsing error: Failed to convert 字体加载失败
字体文件,加载不出来 解决方案 一 问题是Maven正在过滤字体文件并破坏它们. <resource> <directory>${project.basedir}/src/m ...
- Oracle批量操作数据库
1:批量插入 <insert id="insertBatch" parameterType="Java.util.List" > insert in ...
- Python自动化--语言基础6--模块操作之re、MySQL、Excel
1.Python自有模块正则 import re # re.match只匹配字符串的开始,如果字符串开始不符合正则表达式,则匹配失败,函数返回None print(re.match("aaa ...
- solr6.6教程-从mysql数据库中导入数据(三)
整理多半天的solr6.6,终于算是把solr6.6中配置数据库弄出来了,网上的文章千篇一律,各说个的,没有一篇统一覆盖solr6.6版本配置的 帖子,本章节我会把配置的一些注意事项一一列举,由于时间 ...
- thinkpad E480 用户初体验
梦寐以求的E480终于到手了,经过几天使用之后,也该写一下用户体验了,也算是对这些天的调研的一个总结吧. 首先是外形,从外形上来看E480的外形的确是相对以前的塑料外壳,改进了不少,全金属的机身,磨砂 ...
- uva10003 区间DP
很清晰的区间dp问题.d(i,j)表示断点i到断点j的最小费用,由于开头和结尾也是断点,所以应该加入断点数组,即 cut[0]=0; cut[n+1]=len; 边界就是d(i,i+1)=0; 转移方 ...
- poj1011 && uva307 DFS + 剪枝
将木棒从大到小排列,保证每次的选择都是最长可选的木棒. 剪枝: 1 . 如果第 i 根木棒被选择却无法成功拼接,那么后面与其长度相同的也不能选择. 2 . 如果第 cnt + 1 根木棒无法成功拼接, ...
- Linux sftp 安全文件传输命令
sftp 是一个交互式文件传输程式.它类似于 ftp, 但它进行加密传输,比FTP有更高的安全性. 1.常用登陆方式: 格式:sftp <user>@<host> 通过sftp ...