51NOD1965:奇怪的式子
传送门
拆开变成
\]
考虑 \(\prod_{i=1}^{n}\sigma_0(i)^{\mu(i)}\)
运用 \(\mu\) 的性质,设 \(c(i)\) 表示 \(i\) 的质数因子个数
那么就是
\]
只需要设
\(f(x,j)\) 表示 \(1\) 到 \(x\) 最小质因子大于等于第 \(j\) 个质数的合数的 \(\mu(i)c(i)\) 的和
\(g(x,j)\) 表示 \(1\) 到 \(x\) 最小质因子大于等于第 \(j\) 个质数的合数的 \(\mu(i)\) 的和
预处理出 \(cntp(i)\) 表示 \(1\) 到 \(x\) 的质数个数就可以递推了
这一部分直接 \(min25\) 筛即可
考虑 \(\prod_{i=1}^{n}\sigma_0(i)^{i}\)
分开考虑每个质数的贡献
设 \(s(n,p,k)\) 表示 \(1\) 到 \(n\) 中只含有 \(p^k\) 这个因子不含 \(p^i,i\ne k\) 的数的和
设 \(P\) 为质数集合
那么
\]
设 \(Sum(x)=\sum_{i=1}^{x}i\)
那么
\]
当 \(p\le \sqrt{n}\) 的时候直接暴力计算
当 \(p > \sqrt{n}\) 的时候,显然 \(k\) 最多就是 \(1\)
那么
\]
\(min25\) 筛预处理出 \(S(x)\) 表示 \(1\) 到 \(x\) 的质数的和,然后数论分块即可
注意常数优化QwQ
# include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll mod(1e12 + 39);
const ll phi(1e12 + 38);
const int maxn(1e6 + 5);
inline ll Mul1(ll x, ll y) {
register ll ret = x * y - (ll)((long double)x / mod * y + 0.5) * mod;
return ret < 0 ? ret + mod : ret;
}
inline ll Mul2(ll x, ll y) {
register ll ret = x * y - (ll)((long double)x / phi * y + 0.5) * phi;
return ret < 0 ? ret + phi : ret;
}
inline ll Pow(ll x, ll y) {
register ll ret = 1;
for (; y; y >>= 1, x = Mul1(x, x))
if (y & 1) ret = Mul1(ret, x);
return ret;
}
inline void Inc1(ll &x, ll y) {
x = x + y >= mod ? x + y - mod : x + y;
}
inline void Inc2(ll &x, ll y) {
x = x + y >= phi ? x + y - phi : x + y;
}
inline ll Dec1(ll x, ll y) {
return x - y < 0 ? x - y + mod : x - y;
}
inline ll Dec2(ll x, ll y) {
return x - y < 0 ? x - y + phi : x - y;
}
int test, pr[maxn], tot, d, id1[maxn], id2[maxn], cnt;
ll n, val[maxn], cntp[maxn], f[maxn], g[maxn], s[maxn], sp[maxn];
bitset <maxn> ispr;
# define ID(x) ((x) <= d ? id1[x] : id2[n / (x)])
inline ll Sum1(ll x) {
return (x & 1) ? Mul1(x, (x + 1) / 2) : Mul1(x / 2, x + 1);
}
inline ll Sum2(ll x) {
return (x & 1) ? Mul2(x, (x + 1) / 2) : Mul2(x / 2, x + 1);
}
inline ll GetSum(ll x, ll v1, ll v2) {
return Dec2(Mul2(v1, Sum2(x / v1)), Mul2(v2, Sum2(x / v2)));
}
inline void Solve() {
while (cnt) f[cnt] = g[cnt] = 0, cnt--;
register ll i, j, id, ans1, ans2, v, ret, lst, cur;
for (d = sqrt(n), i = 1; i <= n; i = j + 1) {
j = n / (n / i), val[++cnt] = n / i;
val[cnt] <= d ? id1[val[cnt]] = cnt : id2[j] = cnt;
cntp[cnt] = val[cnt] - 1, s[cnt] = Sum2(val[cnt]) - 1;
}
for (i = 1; i <= tot && pr[i] <= n / pr[i]; ++i)
for (j = 1; j <= cnt && pr[i] <= val[j] / pr[i]; ++j) {
id = ID(val[j] / pr[i]), cntp[j] -= cntp[id] - i + 1;
Inc2(s[j], Mul2(pr[i], Dec2(sp[i - 1], s[id])));
}
for (--i; i; --i)
for (j = 1; j <= cnt && pr[i] <= val[j] / pr[i]; ++j) {
id = ID(val[j] / pr[i]);
v = Dec2(cntp[id] - i, g[id]), Inc2(g[j], v);
Inc2(v, cntp[id] - i), Inc2(f[j], Dec2(v, f[id]));
}
for (i = 1; i <= cnt; ++i) Inc2(f[i], phi - cntp[i]);
ans1 = Pow(2, f[ID(n)]), ans2 = 1;
for (i = 1; i <= tot && pr[i] <= n / pr[i]; ++i)
for (j = 1, v = pr[i]; v <= n; v = v * pr[i], ++j) ans2 = Mul1(ans2, Pow(j + 1, GetSum(n, v, v * pr[i])));
for (lst = sp[i - 1], ret = 0, i = pr[i]; i <= n; i = j + 1) {
j = n / (n / i);
Inc2(ret, Mul2(Dec2(cur = s[ID(j)], lst), Sum2(n / i))), lst = cur;
}
ans2 = Mul1(ans2, Pow(2, ret)), printf("%lld\n", Mul1(ans1, ans2));
}
int main() {
register int i, j;
ispr[1] = 1;
for (i = 1; i < maxn; ++i) {
if (!ispr[i]) pr[++tot] = i, sp[tot] = (sp[tot - 1] + i) % phi;
for (j = 1; j <= tot && i * pr[j] < maxn; ++j) {
ispr[i * pr[j]] = 1;
if (!(i % pr[j])) break;
}
}
scanf("%d", &test);
while (test) scanf("%lld", &n), Solve(), --test;
return 0;
}
51NOD1965:奇怪的式子的更多相关文章
- [51nod1965]奇怪的式子
noteskey 怎么说,魔性的题目...拿来练手 min_25 正好...吧 首先就是把式子拆开来算贡献嘛 \[ANS=\prod_{i=1}^n \sigma_0(i)^{\mu(i)} \pro ...
- 51nod1965. 奇怪的式子(min_25筛)
题目链接 http://www.51nod.com/Challenge/Problem.html#!#problemId=1965 题解 需要求的式子显然是个二合一形式,我们将其拆开,分别计算 \(\ ...
- 【51nod1965】奇怪的式子
Portal --> 51nod1965 Solution 怎么说呢..这题..做的有点痛苦.. 首先看这个式子长得..比较奇怪,指数里面那个加号有点烦人,而且这个函数不是一个积性函数也有点烦人 ...
- 【51NOD1965】奇怪的式子 min_25筛
题目描述 给你\(n\),求 \[ \prod_{i=1}^n{\sigma_0(i)}^{i+\mu(i)} \] 对\({10}^{12}+39\)取模. \(\sigma_0(i)\)表示约数个 ...
- 51nod 1965 奇怪的式子 —— min_25筛
题目:http://www.51nod.com/Challenge/Problem.html#!#problemId=1965 推式子就同这里:https://www.cnblogs.com/yoyo ...
- 51nod 1965 奇怪的式子——min_25筛
题目:http://www.51nod.com/Challenge/Problem.html#!#problemId=1965 考虑 \( \prod_{i=1}^{n}\sigma_0^i \) \ ...
- light oj 1236 分解质因数
题目链接:http://acm.hust.edu.cn/vjudge/contest/view.action?cid=70017#problem/H 题意:求满足1<=i<=j<=n ...
- BZOJ 1045 糖果传递
奇怪的式子.最后发现取中位数. #include<iostream> #include<cstdio> #include<cstring> #include< ...
- BZOJ3456城市规划
题目描述 刚刚解决完电力网络的问题, 阿狸又被领导的任务给难住了.刚才说过, 阿狸的国家有n个城市, 现在国家需要在某些城市对之间建立一些贸易路线, 使得整个国家的任意两个城市都直接或间接的连通.为了 ...
随机推荐
- 基于注解的Spring容器源码分析
从spring3.0版本引入注解容器类之后,Spring注解的使用就变得异常的广泛起来,到如今流行的SpringBoot中,几乎是全部使用了注解.Spring的常用注解有很多,有@Bean,@Comp ...
- linux C API连接并查询mysql5.7.9
开发环境: ubuntu16.04 mysql5.7.9 原生C API VIM 配置远程连接 配置mysql允许远程连接的方法默认情况下,mysql只允许本地登录,如果要开启远程连接,则需要修改/e ...
- Asp.net的生命周期之应用程序生命周期
参考:http://msdn.microsoft.com/zh-cn/library/ms178473(v=vs.100).aspx 参考:http://www.cnblogs.com/JimmyZh ...
- 可以修改类不用重启Tomcat加载整个项目
修改类不重启Tomcat(不用手动重启) 修改tomcat conf目录下的server.xml <Context path="/struts2" docBase=&quo ...
- 博客主题皮肤探索-主题的本地化和ChromeCacheView使用
还有前言 用了大佬的主题之后,有些资源是无法在暴露的接口处更改的,需要自己去css更改.但后台给的都是压缩过的,找起来比较困难,所以特意记录了这篇. 本地资源替换 侧边栏: .introduce-bo ...
- mongoengine在python中的使用
# /usr/bin/python # encoding: utf-8 # Author: masako from mongoengine import * host ='127.0.0.1' por ...
- 全面解析C#中参数传递
一.引言 对于一些初学者(包括工作几年的人在内)来说,有时候对于方法之间的参数传递的问题感觉比较困惑的,因为之前在面试的过程也经常遇到参数传递的基础面试题,这样的面试题主要考察的开发人员基础是否扎实, ...
- (转)AIX的SVMON命令详解
原文:http://czmmiao.iteye.com/blog/1153499 svmon概述 svmon 命令用于显示当前内存状态的信息,可通过 # lslpp bos.perf.tools 查看 ...
- Steps to install Docker on Manjaro 16.10--转
https://manjaro-tutorial.blogspot.com/2016/12/how-to-install-docker-on-manjaro-1610.html Open Termin ...
- WindowsServer2008安装IIS相关服务
控制面板->程序->打开或关闭Windows功能 添加角色,选择IIS服务器,选择以下角色服务,如果添加过就选择添加角色服务 如果出现500错误,这个跟程序没有多大关系,可以试一下以下操作 ...