洛谷题面传送门

hot tea.

首先注意到这个 \(\text{lcm}\) 特别棘手,并且这里的 \(k\) 大得离谱,我们也没办法直接枚举每个质因子的贡献来计算答案。不过考虑到如果我们把这里的 \(\text{lcm}\) 改为 \(\gcd\) 那么一遍莫比乌斯反演即可搞定,因此考虑将这里的 \(\text{lcm}\) 与 \(\gcd\) 联系在一起。那么什么能将这两个东西联系在一起呢?Min-Max 容斥,具体来说,考虑式子

\[\text{lcm}(S)=\prod\limits_{\varnothing\ne T\subseteq S}\gcd(T)^{(-1)^{|T|-1}}
\]

因此我们考虑开始推式子:

\[\begin{aligned}
res&=\prod\limits_{i_1,i_2,\cdots,i_k}\text{lcm}(i_1,i_2,\cdots,i_k)\\
&=\prod\limits_{i_1,i_2,\cdots,i_k}\prod\limits_{S\subseteq\{1,2,3,\cdots,k\}}\gcd(\{i_x|x\in S\})^{(-1)^{|S|-1}}\\
&=\prod\limits_{j=1}^k\prod\limits_{|S|=j,S\subseteq\{1,2,3,\cdots,k\}}\prod\limits_{i_1,i_2,\cdots,i_k}\gcd(\{i_x|x\in S\})^{(-1)^{j-1}}\\
&=\prod\limits_{j=1}^k\prod\limits_{i_1,i_2,\cdots,i_j}\gcd(i_1,i_2,\cdots,i_j)^{(-1)^{j-1}·\dbinom{k}{j}·n^{k-j}}\\
&=\prod\limits_{j=1}^k\prod\limits_{d=1}^nd^{(-1)^{j-1}·\dbinom{k}{j}·n^{k-j}·\sum\limits_{i_1,i_2,\cdots,i_j}[\gcd(i_1,i_2,\cdots,i_j)=d]}\\
&=\prod\limits_{j=1}^k\prod\limits_{d=1}^nd^{(-1)^{j-1}·\dbinom{k}{j}·n^{k-j}·\sum\limits_{p=1}^{n/d}\mu(p)\lfloor\dfrac{n}{dp}\rfloor^j}\\
&=\prod\limits_{d=1}^n\prod\limits_{p=1}^n(d^{\mu(p)})^{\sum\limits_{j=1}^k(-1)^{j-1}·n^{k-j}·\lfloor\dfrac{n}{dp}\rfloor^j}\\
&=\dfrac{1}{\prod\limits_{d=1}^n\prod\limits_{p=1}^n(d^{\mu(p)})^{\sum\limits_{j=1}^k(-1)^{j}·n^{k-j}·\lfloor\dfrac{n}{dp}\rfloor^j}}\\
&=\dfrac{1}{\prod\limits_{d=1}^n\prod\limits_{p=1}^n(d^{\mu(p)})^{\sum\limits_{j=0}^k(-1)^{j}·n^{k-j}·\lfloor\dfrac{n}{dp}\rfloor^j-n^k}}\\
&=\dfrac{1}{\prod\limits_{d=1}^n\prod\limits_{p=1}^n(d^{\mu(p)})^{(n-\lfloor\dfrac{n}{dp}\rfloor)^k-n^k}}\\
&=\dfrac{1}{\prod\limits_{T=dp}f(T)^{(n-\lfloor\dfrac{n}{T}\rfloor)^k-n^k}}
\end{aligned}
\]

其中 \(f(T)=\prod\limits_{dp=T}d^{\mu(p)}\)。

整除分块求解即可,注意,根据费马小定理,指数上那托东西要 \(\bmod 998244352\),而由于指数上面还套了一个快速幂,且这里的 \(k\) 很大,因此需要欧拉降幂,具体来说如果 \(k\le\varphi(998244352)\) 那么直接算即可,否则需令 \(k\leftarrow k\bmod\varphi(998244352)+\varphi(998244352)\)

时间复杂度 \(T\sqrt{n}\log n\)

const int MAXN=1e6;
const int MAXLEN=100;
const int MOD=998244353;
const int PHI=MOD-1;
const int PPHI=402653184;
int getmod(ll x){return (x>PPHI)?(x%PPHI+PPHI):x;}
int pr[MAXN+5],prcnt=0,mu[MAXN+5],prd[MAXN+5];
bool vis[MAXN+5];
int qpow(int x,int e){
if(e<0) e+=PHI;int ret=1;
for(;e;e>>=1,x=1ll*x*x%MOD) if(e&1) ret=1ll*ret*x%MOD;
return ret;
}
int _qpow(int x,int e){
int ret=1;
for(;e;e>>=1,x=1ll*x*x%PHI) if(e&1) ret=1ll*ret*x%PHI;
return ret;
}
void sieve(int n){
mu[1]=1;
for(int i=2;i<=n;i++){
if(!vis[i]) mu[i]=-1,pr[++prcnt]=i;
for(int j=1;j<=prcnt&&pr[j]*i<=n;j++){
vis[pr[j]*i]=1;if(i%pr[j]==0) break;
mu[pr[j]*i]=-mu[i];
}
}
for(int i=1;i<=n;i++) prd[i]=1;
for(int i=1;i<=n;i++) for(int j=1;j*i<=n;j++)
prd[i*j]=1ll*prd[i*j]*qpow(i,mu[j])%MOD;
// for(int i=1;i<=n;i++) printf("%d\n",prd[i]);
prd[0]=1;
for(int i=1;i<=n;i++) prd[i]=1ll*prd[i-1]*prd[i]%MOD;
}
int n;ll k;
char buf[MAXLEN+5];
void solve(){
scanf("%d%s",&n,buf+1);k=0;int len=strlen(buf+1);
for(int i=1;i<=len;i++) k=getmod(10ll*k+buf[i]-'0');
int res=1;
for(int l=1,r;l<=n;l=r+1){
r=(n/(n/l));
res=1ll*res*qpow(1ll*prd[r]*qpow(prd[l-1],MOD-2)%MOD,(_qpow(-n/l+n,k)-_qpow(n,k)+PHI)%PHI)%MOD;
} printf("%d\n",qpow(res,MOD-2));
}
int main(){
int qu;scanf("%d",&qu);sieve(MAXN);
while(qu--) solve();
return 0;
}

洛谷 P7360 -「JZOI-1」红包(Min-Max 容斥+推式子)的更多相关文章

  1. [洛谷P3701]「伪模板」主席树

    题目大意:太暴力了,就不写了,看这儿 题解:对于每个$byx$的人,从源点向人连边,容量为此人的寿命. 对于每个手气君的人,从人向汇点连边,容量为此人的寿命. 对于每个$byx$的人与手气君的人,如果 ...

  2. LOJ 3119: 洛谷 P5400: 「CTS2019 | CTSC2019」随机立方体

    题目传送门:LOJ #3119. 题意简述: 题目说的很清楚了. 题解: 记恰好有 \(i\) 个极大的数的方案数为 \(\mathrm{cnt}[i]\),则答案为 \(\displaystyle\ ...

  3. LOJ 3120: 洛谷 P5401: 「CTS2019 | CTSC2019」珍珠

    题目传送门:LOJ #3120. 题意简述: 称一个长度为 \(n\),元素取值为 \([1,D]\) 的整数序列是合法的,当且仅当其中能够选出至少 \(m\) 对相同元素(不能重复选出元素). 问合 ...

  4. 【洛谷5643】[PKUWC2018] 随机游走(Min-Max容斥+待定系数法+高维前缀和)

    点此看题面 大致题意: 从一个给定点出发,在一棵树上随机游走,对于相邻的每个点均有\(\frac 1{deg}\)的概率前往.多组询问,每次给出一个点集,求期望经过多少步能够访问过点集内所有点至少一次 ...

  5. LOJ #3119「CTS2019 | CTSC2019」随机立方体 (容斥)

    博客链接 里面有个下降幂应该是上升幂 还有个bk的式子省略了k^3 CODE 蛮短的 #include <bits/stdc++.h> using namespace std; const ...

  6. 洛谷 P2522 [HAOI2011]Problem b (莫比乌斯反演+简单容斥)

    题目描述 对于给出的n个询问,每次求有多少个数对(x,y),满足a≤x≤b,c≤y≤d,且gcd(x,y) = k,gcd(x,y)函数为x和y的最大公约数. 输入输出格式 输入格式: 第一行一个整数 ...

  7. 洛谷 1447 [NOI2010]能量采集——容斥/推式子

    题目:https://www.luogu.org/problemnew/show/P1447 1.容斥原理 求 f [ i ] 表示 gcd==i 的对数,先 f [ i ] = (n/i) * (m ...

  8. 洛谷 P4710 「物理」平抛运动

    洛谷 P4710 「物理」平抛运动 洛谷传送门 题目描述 小 F 回到班上,面对自己 28 / 110 的物理,感觉非常凉凉.他准备从最基础的力学学起. 如图,一个可以视为质点的小球在点 A(x_0, ...

  9. 洛谷比赛 「EZEC」 Round 4

    洛谷比赛 「EZEC」 Round 4 T1 zrmpaul Loves Array 题目描述 小 Z 有一个下标从 \(1\) 开始并且长度为 \(n\) 的序列,初始时下标为 \(i\) 位置的数 ...

随机推荐

  1. Node.js CMS——基于 NestJS/NuxtJS 的完整开源项目

    这是一款轻量级的基于 Node.js 的开源 CMS,采用前后端分离开发模式,集成了 API.后台管理.WEB 展示三个完整项目.开箱即是一套完整的企业网站,适合企业.个人直接使用或二次开发. API ...

  2. JavaScript05

    显示和隐藏 元素的显示和隐藏 元素display属性可控制元素的显示和隐藏,先获取元素对象,再通过点语法调用style对象中的display属性 语法格式: 元素.style.display='non ...

  3. Java:并发笔记-05

    Java:并发笔记-05 说明:这是看了 bilibili 上 黑马程序员 的课程 java并发编程 后做的笔记 4. 共享模型之内存 本章内容 上一章讲解的 Monitor 主要关注的是访问共享变量 ...

  4. Uniapp云打包生成apk下载链接

    使用uni[]()app云打包生成安装包下载链接 manifest.json 中配置自动获取appid manifest.json中配置app 图标 按教程生成.keystore证书 使用云打包生成安 ...

  5. AFO记

    希望永远也不要动笔写这个. 发以自勉

  6. 如何优雅的处理 accept 出现 EMFILE 的问题

    通常情况下,服务端调用 accept 函数会返回一个新的文件描述符,用于和客户端之间的数据传输 在服务器的开发中,有时会遇到这种情况:当调用 accept 函数接受客户端连接,函数返回失败,对应的错误 ...

  7. 矩形覆盖 牛客网 剑指Offer

    矩形覆盖 牛客网 剑指Offer 题目描述 我们可以用21的小矩形横着或者竖着去覆盖更大的矩形.请问用n个21的小矩形无重叠地覆盖一个2*n的大矩形,总共有多少种方法? class Solution: ...

  8. 我为啥开始用CSDN博客

    今晚开通CSDN博客,并且决定以后每天都使用这个不错的东西.与此同时,在博客园也开通了一个:http://www.cnblogs.com/fish7/ 我原本是把做过的题都用WPS整理的,然后每次打印 ...

  9. Centos 系统常用编译环境

    centos编译环境配置 yum install -y autoconf make automake gcc gcc-c++

  10. swoole、swoft环境配置

    一.服务器环境 1.lnmp wget http://soft.vpser.net/lnmp/lnmp1.5.tar.gz -cO lnmp1.5.tar.gz && tar zxf ...