luogu P3704 [SDOI2017]数字表格
我是真的弱,推式子只能推一半
下面假设\(n<m\)
考虑题目要求的东西,可以考虑每个gcd的贡献,即$$\prod_{d=1}{n}f[d]{\sum_{i=1}{\lfloor\frac{n}{d}\rfloor}\sum_{j=1}{\lfloor\frac{m}{d}\rfloor}[gcd(i,j)=1]}$$
由\(n=\sum_{d|n} \mu[d]\),得$$\prod_{d=1}{n}f[d]{\sum_{i=1}{\lfloor\frac{n}{d}\rfloor}\sum_{j=1}{\lfloor\frac{m}{d}\rfloor}\sum_{k|i,k|j}\mu[k]}$$$$\prod_{d=1}{n}f[d]{\sum_{k=1}^{\lfloor\frac{n}{d}\rfloor}\mu[k]\lfloor\frac{n}{kd}\rfloor\lfloor\frac{m}{kd}\rfloor}$$
大力数论分块即可获得60'好成绩
如果我们令\(t=kd\),然后把t提出来,即则$$\prod_{t=1}{n}\prod_{d|t}f[d]{\lfloor\frac{n}{t}\rfloor\lfloor\frac{m}{t}\rfloor\mu[\lfloor\frac{n}{d}\rfloor]}$$$$\prod_{t=1}{n}(\prod_{d|t}f[d]{\mu[\lfloor\frac{n}{d}\rfloor]})^{\lfloor\frac{n}{t}\rfloor\lfloor\frac{m}{t}\rfloor}$$
里面的可以枚举每个数倍数预处理,然后就是数论分块
// luogu-judger-enable-o2
#include<bits/stdc++.h>
#define LL long long
#define ldb long double
#define il inline
#define re register
using namespace std;
const int N=1e6+10,mod=1e9+7;
il int rd()
{
int x=0,w=1;char ch=0;
while(ch<'0'||ch>'9') {if(ch=='-') w=-1;ch=getchar();}
while(ch>='0'&&ch<='9') {x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
return x*w;
}
il int fpow(int a,int b){int an=1;while(b){if(b&1) an=1ll*an*a%mod;a=1ll*a*a%mod,b>>=1;}return an;}
bool v[N];
int prm[N],mu[N],tt,fb[N],ifb[N],f[N],ivf[N];
int main()
{
fb[0]=0,fb[1]=1,ifb[0]=ifb[1]=1;
for(int i=2;i<=N-10;++i) fb[i]=(fb[i-1]+fb[i-2])%mod,ifb[i]=fpow(fb[i],mod-2);
mu[1]=1;
for(int i=2;i<=N-10;++i)
{
if(!v[i]) v[i]=1,prm[++tt]=i,mu[i]=-1;
for(int j=1;j<=tt&&i*prm[j]<=N-10;++j)
{
v[i*prm[j]]=1,mu[i*prm[j]]=-mu[i];
if(i%prm[j]==0) {mu[i*prm[j]]=0;break;}
}
}
for(int i=0;i<=N-10;++i) f[i]=1;
for(int i=1;i<=N-10;++i)
for(int j=1;i*j<=N-10;++j)
f[i*j]=1ll*f[i*j]*((~mu[j])?fpow(fb[i],mu[j]):ifb[i])%mod;
ivf[0]=1;
for(int i=1;i<=N-10;++i) f[i]=1ll*f[i]*f[i-1]%mod,ivf[i]=fpow(f[i],mod-2);
int T=rd();
while(T--)
{
int n=rd(),m=rd(),ans=1;
if(n>m) swap(n,m);
for(int i=1,j;i<=n;i=j+1)
{
j=min(n/(n/i),m/(m/i));
ans=1ll*ans*fpow(1ll*f[j]*ivf[i-1]%mod,1ll*(m/i)*(n/i)%(mod-1))%mod;
}
printf("%d\n",ans);
}
return 0;
}
luogu P3704 [SDOI2017]数字表格的更多相关文章
- P3704 [SDOI2017]数字表格
P3704 [SDOI2017]数字表格 链接 分析: $\ \ \ \prod\limits_{i = 1}^{n} \prod\limits_{j = 1}^{m} f[gcd(i, j)]$ $ ...
- P3704 [SDOI2017]数字表格 (莫比乌斯反演)
[题目链接] https://www.luogu.org/problemnew/show/P3704 [题解] https://www.luogu.org/blog/cjyyb/solution-p3 ...
- bzoj 4816: 洛谷 P3704: [SDOI2017]数字表格
洛谷很早以前就写过了,今天交到bzoj发现TLE了. 检查了一下发现自己复杂度是错的. 题目传送门:洛谷P3704. 题意简述: 求 \(\prod_{i=1}^{N}\prod_{j=1}^{M}F ...
- 洛谷P3704 [SDOI2017]数字表格
题目描述 Doris刚刚学习了fibonacci数列.用f[i]f[i] 表示数列的第ii 项,那么 f[0]=0f[0]=0 ,f[1]=1f[1]=1 , f[n]=f[n-1]+f[n-2],n ...
- 洛谷 P3704 [SDOI2017]数字表格(莫比乌斯函数)
题面传送门 题意: 求 \[\prod\limits_{i=1}^n\prod\limits_{j=1}^mfib_{\gcd(i,j)} \] \(T\) 组测试数据,\(1 \leq T \leq ...
- 洛咕 P3704 [SDOI2017]数字表格
大力推式子 现根据套路枚举\(\gcd(i,j)\) \(ans=\Pi_{x=1}^nfib[x]^{\sum_{i=1}^{n/x}\sum_{j=1}^{n/x}[\gcd(i,j)=1]}\) ...
- Luogu 3704 [SDOI2017]数字表格
列一下式子: $\prod_{i = 1}^{n}\prod_{j = 1}^{m}fib_{gcd(i, j)}$ 很套路的变成这样: $\prod_{d = 1}^{min(n, m)}fib_{ ...
- 洛谷P3704 [SDOI2017]数字表格(莫比乌斯反演)
传送门 yyb大佬太强啦…… 感觉还是有一点地方没有搞懂orz //minamoto #include<cstdio> #include<iostream> #include& ...
- 并不对劲的bzoj4816:loj2000:p3704[SDOI2017]数字表格
题目大意 有函数\(f(x)\),\(f(0)=0,f(1)=1,f(x)=f(x-1)+f(x-2)\) \(t\)(\(t\leq1000\))组询问,每次给定\(n,m\)(\(n,m\leq1 ...
随机推荐
- 基于Jenkins,docker实现自动化部署(持续交互)
前言 随着业务的增长,需求也开始增多,每个需求的大小,开发周期,发布时间都不一致.基于微服务的系统架构,功能的叠加,对应的服务的数量也在增加,大小功能的快速迭代,更加要求部署的快速化,智能化.因此 ...
- Django cookie相关操作
Django cookie 的相关操作还是比较简单的 首先是存储cookie #定义设置cookie(储存) def save_cookie(request): #定义回应 response = Ht ...
- dispatchEvent(AWTEvent) 分派事件
点一个按钮,显示的分派一个指定的事件给系统. 下面是一个例子,当点击close按钮时,分派一个new WindowEvent(this,WindowEvent.WINDOW_CLOSING)事件给系统 ...
- Spring2
简介:1.Aop编程.2.AspectJ基于xml文件.3.AspectJ基于注解. 4.JdbcTemplate. 5.配置properties文件 1 AOP 1.1 AOP介绍 ...
- Luogu P2292 [HNOI2004]L语言
题目链接 \(Click\) \(Here\) 好久没写\(DP\)了真是水平下降不少,一眼把这个题搞成贪心了,然后一发交上只有\(37\)分\(QwQ\) 这个题好像还可以\(AC\)自动机胡搞?不 ...
- linux基本
一.初识 Linux与windows相比的优点是:长期稳定的运行,避免了因为系统的问题导致的项目运行中断:占用资源少:开源软件多. Centos(Community Enterprise Operat ...
- ansible小结常用模块
根据官方的分类,将模块按功能分类为:云模块.命令模块.数据库模块.文件模块.资产模块.消息模块.监控模块.网络模块.通知模块.包管理模块.源码控制模块.系统模块.单元模块.web设施模块.window ...
- C++ 二维数组作为形参传递使用实例
在线代码编辑器: http://codepad.org/ 1.*指针 void display(int *arr, const int row, const int col) { ; i < r ...
- python django基础三 模版渲染
request对象 当一个页面被请求时,Django就会创建一个包含本次请求原信息的HttpRequest对象.Django会将这个对象自动传递给响应的视图函数,一般视图函数约定俗成地使用 reque ...
- canvas svg webgl threejs d3js 的区别
canvas 绘制2D位图. Echarts是基于Canvas技术的可视化工具,底层封装了原生的JavaScript的绘图 API. canvas里面绘制的图形不能被引擎抓取,canvas中我们绘制图 ...