Description

  



  

  

  

Solution

  

  一开始的时候我是这么推的(\(f(n)\)表示斐波那契数列的第\(n\)项)

\[\begin{aligned}
Ans&=\prod_{x=1}^{\min(n,m)}f(x)^{(\sum_{i=1}^n\sum_{j=1}^m[\gcd(i,j)=x])}\\
&=\prod_{x=1}^{\min(n,m)}f(x)^{\sum_{e=1}^{\min(\lfloor\frac nx\rfloor,\lfloor \frac mx\rfloor)}\mu(e)\lfloor\frac n {ex}\rfloor\lfloor\frac m {ex}\rfloor}\\
\end{aligned}
\]

​  然后我想,根号分段套根号分段,\(\mathcal O(\sqrt n(\log+\sqrt n))\)解决!嗯不错,一看数据组数1000........不过还是能拿70的分。

  

  

  

  题解的思路非常神。

  

  假设我们能构造一个函数\(g\),使得

\[f(n)=\prod_{d|n}g(d)
\]

​  那么答案就变成

\[\begin{aligned}
Ans&=\prod_{i=1}^n\prod_{j=1}^mf({\gcd(i,j))}\\
&=\prod_{i=1}^n\prod_{j=1}^m\prod_{d|i,d|j}g(d)\\
&=\prod_{d=1}^{\min(n,m)}g(d)^{\lfloor\frac nd\rfloor\lfloor\frac md\rfloor}
\end{aligned}
\]

  这样就可以在\(\mathcal O(2\sqrt n)\)的时间内处理每一个询问了。前提是我们知道\(g\)及其前缀积。

  

​  考虑式子\(f(n)=\prod_{d|n}g(d)\)十分像莫比乌斯反演,能否用类似的形式反演出\(g\)呢?

  

​  在\(\sum\)的意义下

\[g_n=\sum_{d|n}\mu(d)f({\frac nd})
\]

​  反演的本质是通过加减来容斥出所需要的组合。而在乘法的意义下,不就是通过乘除来容斥出所需要的组合吗?所以有:

\[g_n=\prod_{d|n}f({\frac nd})^{\mu(d)}=\prod_{d|n}f(d)^{\mu(\frac nd)}
\]

​  因此我们可以在\(\mathcal O(n \lg n)\)的时间内处理出\(g\)的取值和前缀积。那么上面的根号分段也就迎刃而解了。

  

  总体思路是仿造莫比乌斯反演构造一个可求函数,利用该函数化简式子使得答案求和式变得简明且复杂度较低,再通过传统根号分段求解。

  

  

  

Code

  

#include <cstdio>
using namespace std;
const int N=1e6+10,MOD=1e9+7,PMOD=MOD-1;
int n,m;
int fib[N],ifib[N],g[N],ig[N];
bool vis[N];
int p[N],pcnt,mu[N];
inline int min(int x,int y){return x<y?x:y;}
inline int fmi(int x,int y){
int res=1;
for(;y;x=1LL*x*x%MOD,y>>=1)
if(y&1) res=1LL*res*x%MOD;
return res;
}
void sieve(){
mu[1]=1;
for(int i=2;i<=1e6;i++){
if(!vis[i]) p[++pcnt]=i,mu[i]=-1;
for(int j=1;j<=pcnt&&i*p[j]<=1e6;j++){
int x=i*p[j];
vis[x]=true;
if(i%p[j]==0){
mu[x]=0;
break;
}
mu[x]=-mu[i];
}
}
}
void prework(){
sieve();
fib[0]=0; fib[1]=1; ifib[1]=1;
for(int i=2;i<=1e6;i++){
fib[i]=(fib[i-2]+fib[i-1])%MOD;
ifib[i]=fmi(fib[i],MOD-2);
}
for(int i=1;i<=1e6;i++) g[i]=1;
for(int d=1;d<=1e6;d++)
for(int n=d;n<=1e6;n+=d)
g[n]=1LL*g[n]*(mu[n/d]==1?fib[d]:(mu[n/d]==-1?ifib[d]:1))%MOD;
g[0]=ig[0]=1;
for(int i=1;i<=1e6;i++){
g[i]=1LL*g[i]*g[i-1]%MOD;
ig[i]=fmi(g[i],MOD-2);
}
}
int main(){
prework();
int T,up,ans;
scanf("%d",&T);
while(T--){
scanf("%d%d",&n,&m);
up=min(n,m);
ans=1;
for(int i=1,j;i<=up;i=j+1){
j=min(n/(n/i),m/(m/i));
ans=1LL*ans*fmi(1LL*g[j]*ig[i-1]%MOD,1LL*(n/i)*(m/i)%PMOD)%MOD;
}
printf("%d\n",ans<0?ans+MOD:ans);
}
return 0;
}

数字表格(product)的更多相关文章

  1. 数字表格(product)

    Portal -->broken qwq Description ​  求\(\prod\limits_{i=1}^n\prod\limits_{j=1}^m f[gcd(i,j)](mod\ ...

  2. [Sdoi2017]数字表格 [莫比乌斯反演]

    [Sdoi2017]数字表格 题意:求 \[ \prod_{i=1}^n \prod_{j=1}^m f[(i,j)] \] 考场60分 其实多推一步就推倒了... 因为是乘,我们可以放到幂上 \[ ...

  3. BZOJ 2154: Crash的数字表格 [莫比乌斯反演]

    2154: Crash的数字表格 Time Limit: 20 Sec  Memory Limit: 259 MBSubmit: 2924  Solved: 1091[Submit][Status][ ...

  4. 【BZOJ】【2154】Crash的数字表格

    莫比乌斯反演 PoPoQQQ讲义第4题 题解:http://www.cnblogs.com/jianglangcaijin/archive/2013/11/27/3446169.html 感觉两次sq ...

  5. 【莫比乌斯反演】关于Mobius反演与lcm的一些关系与问题简化(BZOJ 2154 crash的数字表格&&BZOJ 2693 jzptab)

    BZOJ 2154 crash的数字表格 Description 今天的数学课上,Crash小朋友学习了最小公倍数(Least Common Multiple).对于两个正整数a和b,LCM(a, b ...

  6. 【BZOJ 2154】Crash的数字表格 (莫比乌斯+分块)

    2154: Crash的数字表格 Description 今天的数学课上,Crash小朋友学习了最小公倍数(Least Common Multiple).对于两个正整数a和b,LCM(a, b)表示能 ...

  7. BZOJ 4816 数字表格

    首先是惯例的吐槽.SDOI题目名称是一个循环,题目内容也是一个循环,基本上过几年就把之前的题目换成另一个名字出出来,喜大普奔亦可赛艇.学长说考SDOI可以考出联赛分数,%%%. 下面放解题报告.并不喜 ...

  8. BZOJ:4816: [Sdoi2017]数字表格

    4816: [Sdoi2017]数字表格 Time Limit: 50 Sec  Memory Limit: 128 MBSubmit: 501  Solved: 222[Submit][Status ...

  9. 【BZOJ4816】数字表格(莫比乌斯反演)

    [BZOJ4816]数字表格(莫比乌斯反演) 题面 BZOJ 求 \[\prod_{i=1}^n\prod_{j=1}^mf[gcd(i,j)]\] 题解 忽然不知道这个要怎么表示... 就写成这样吧 ...

  10. 【BZOJ2154】Crash的数字表格(莫比乌斯反演)

    [BZOJ2154]Crash的数字表格(莫比乌斯反演) 题面 BZOJ 简化题意: 给定\(n,m\) 求\[\sum_{i=1}^n\sum_{j=1}^mlcm(i,j)\] 题解 以下的一切都 ...

随机推荐

  1. appium 元素定位方法汇总

    以上图为例,要定位到右下角的 我的 ,并点击 # appium的webdriver提供了11种元素定位方法,在selenium的基础上扩展了三个,可以在pycharm里面输入driver.find_e ...

  2. 3.编写sub过程及开发函数——《Excel VBA 程序开发自学宝典》

    3.1 编写sub过程 实例: Sub 建立10个表() If sheets.count>=10 then exit sub Sheets.add , sheets(sheets.count) ...

  3. Centos7 安装与破解 Confluence 6.7.1

    1.1硬件需求建议: CPU:32/64 bit 2.27GHz双核心以上之CPU: 内存:8GB以上: 硬盘:300GB,7200转以上: 建议数据库.Confluence等各自独立一台服务器(本测 ...

  4. Hyperledger Fabric CA User’s Guide——配置设置(四)

    配置设置 Fabric CA提供了三种方案去配置Fabric CA服务端和客户端,优先顺序是: CLI flags(标识) 环境变量 配置文件 在本文档的其余部分中,我们将对配置文件进行更改.但是,可 ...

  5. 【神经网络】自编码聚类算法--DEC (Deep Embedded Clustering)

    1.算法描述 最近在做AutoEncoder的一些探索,看到2016年的一篇论文,虽然不是最新的,但是思路和方法值得学习.论文原文链接 http://proceedings.mlr.press/v48 ...

  6. xocde missing file 解决方法

    创建一个警告里丢失的文件 托到警报的目录下 重启xcode 警告消失 在xcode中选中这个文件右键delete再次删掉文件,就OK了.

  7. ms cms

    很多开源CMS,都是半开源或假开源. 看看微软的CMS怎么样吧.这么著名的CMS,下载网址不好找. https://orchard.codeplex.com/SourceControl/latest

  8. Linux 做网关

    首先创建两张路由表,只需要添加到相应的文件中即可,Linux一共支持255个路由表,rt_tables文件中默认已经存在了三张路由表,分别是:   255  local   254  main   2 ...

  9. 第33次Scrum会议(11/21)【欢迎来怼】

    一.小组信息 队名:欢迎来怼小组成员队长:田继平成员:李圆圆,葛美义,王伟东,姜珊,邵朔,阚博文 小组照片 二.开会信息 时间:2017/11/21 11:35~11:57,总计22min.地点:东北 ...

  10. 冲刺One之站立会议3 /2015-5-16

    2015-5-16 今天我们主要完成一部分服务器端的内容,因为只有服务器端完成了主要功能其他的部分才可以测试有没有成功实现目标.具体包括服务器登陆时需要的端口号.启动时间.服务器状态的显示.在线人数等 ...