https://www.zybuluo.com/ysner/note/1317548---

题面

给出\(n\),用所有长为\(a\)、宽为\(b\)\((1\leq a,b\leq n)\)的长方形拼成正方形,最少需多少块?

多组数据。

  • \(30pts\) \(n\leq100,T\leq100\)
  • \(60pts\) \(n\leq3*10^4,T\leq300\)
  • \(100pts\) \(a,b\leq10^5,T\leq1000\)

解析

显然答案是$$\frac{lcm(a,b)}{a}*\frac{lcm(a,b)}{b}$$

暴力复杂度\(O(Tn^2logn)\),可以通过\(30pts\)。

考虑推推柿子。

\[\prod_{a=1}^n\prod_{b=1}^n\frac{lcm(a,b)}{a}*\frac{lcm(a,b)}{b}
\]

\[=\prod_{a=1}^n\prod_{b=1}^n\frac{lcm^2(a,b)}{ab}
\]

\[=\prod_{a=1}^n\prod_{b=1}^n\frac{ab}{gcd^2(a,b)}
\]

\[=\frac{(n!)^{2n}}{\prod_{a=1}^n\prod_{b=1}^ngcd^2(a,b)}
\]

现在问题是\(\prod_{a=1}^n\prod_{b=1}^ngcd(a,b)\)。

这个复杂度\(O(n^2)\),很不划算。

考虑枚举最大公约数的值。

则化为

\[\prod_{d=1}^nd^{\sum_{a=1}^n\sum_{b=1}^n[gcd(a,b)==d]}
\]

\[\prod_{d=1}^nd^{\sum_{a=1}^{\lfloor\frac{n}{d}\rfloor}\sum_{b=1}^{\lfloor\frac{n}{d}\rfloor}[gcd(a,b)==1]}
\]

这个指数怎么化呢?

我们可以强制\(a>b\),那么指数为\(\sum_{a=1}^{\lfloor\frac{n}{d}\rfloor}\varphi(a)\)

如果不强制,考虑到\(\varphi(1)=1\)的特殊情况,柿子可化为:

\[\prod_{d=1}^nd^{[\sum_{a=1}^{\lfloor\frac{n}{d}\rfloor}2*\varphi(a)]-1}
\]

然后线性预处理一下欧拉函数前缀和,这样复杂度\(O(Tn)\),\(60pts\)稳了。

然后吗,注意到\(\lfloor\frac{n}{d}\rfloor\)在一段区间内是相同的,可以数论分块。

于是就做完了。

预处理逆元后,复杂度\(O(T\sqrt n\))。

#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define ll long long
#define re register
#define il inline
#define fp(i,a,b) for(re int i=a;i<=b;i++)
#define fq(i,a,b) for(re int i=a;i>=b;i--)
using namespace std;
const int mod=19260817,N=1e6+100;
int a,b,jc[N],pri[N],ol[N],n,tot,inv[mod+100];
ll ans,gu;
bool vis[N];
il ll gi()
{
re ll x=0,t=1;
re char ch=getchar();
while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
if(ch=='-') t=-1,ch=getchar();
while(ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar();
return x*t;
}
il ll ksm(re ll S,re ll n)
{
re ll T=S;S=1;
if(n<0) return 0;
while(n)
{
if(n&1) S=S*T%mod;
T=T*T%mod;
n>>=1;
}
return S;
}
il void Pre(re int n)
{
ol[1]=1;
fp(i,2,n)
{
if(!vis[i]) pri[++tot]=i,ol[i]=i-1;
for(re int j=1;j<=tot&&i*pri[j]<=n;++j)
{
vis[i*pri[j]]=1;
if(i%pri[j]) ol[i*pri[j]]=ol[i]*(pri[j]-1);
else {ol[i*pri[j]]=ol[i]*pri[j];break;}
}
}
fp(i,1,n) (ol[i]+=ol[i-1])%=(mod-1);
jc[0]=1;fp(i,1,n) jc[i]=1ll*jc[i-1]*i%mod;
inv[0]=inv[1]=1;fp(i,2,mod) inv[i]=1ll*(mod-mod/i)*inv[mod%i]%mod;
}
int main()
{
re int T=gi();
Pre(1e6);
while(T--)
{
n=gi();
ans=ksm(jc[n],2*n);gu=1;
re int L=0;
for(re int i=1;i<=n;i=L+1)
{
L=n/(n/i);
(gu*=ksm(1ll*jc[L]*inv[jc[i-1]]%mod,2*ol[n/i]-1))%=mod;
}
(ans*=inv[gu*gu%mod])%=mod;
printf("%lld\n",ans);
}
return 0;
}

luogu4917天守阁的地板的更多相关文章

  1. 【洛谷】4917:天守阁的地板【欧拉函数的应用】【lcm与gcd】【同除根号优化】

    P4917 天守阁的地板 题目背景 在下克上异变中,博丽灵梦为了找到异变的源头,一路打到了天守阁 异变主谋鬼人正邪为了迎击,将天守阁反复颠倒过来,而年久失修的天守阁也因此掉下了很多块地板 异变结束后, ...

  2. Luogu 4917 天守阁的地板(莫比乌斯反演+线性筛)

    既然已经学傻了,这个题当然是上反演辣. 对于求积的式子,考虑把[gcd=1]放到指数上.一通套路后可以得到∏D∏d∏i∏j (ijd2)μ(d) (D=1~n,d|D,i,j=1~n/D). 冷静分析 ...

  3. [ Luogu 4917 ] 天守阁的地板

    \(\\\) \(Description\) 定义二元函数\(F(x,y)\)表示,用 \(x\times y\) 的矩形不可旋转的铺成一个任意边长的正方形,所需要的最少的矩形个数. 现在\(T\)组 ...

  4. 省选前的th题

    沙茶博主终于整完了知识点并学完了早该在NOIP之前学的知识们 于是终于开始见题了,之前那个奇怪的题单的结果就是这个了 题目按沙茶博主的做题顺序排序 个人感觉(暂时)意义不大的已被自动忽略 洛谷 491 ...

  5. X000011

    P1890 gcd区间 \(\gcd\) 是满足结合律的,所以考虑用 ST 表解决 时间复杂度 \(O((n\log n+m)\log a_i)\) 考虑到 \(n\) 很小,你也可以直接算出所有的区 ...

  6. Hihocoder 太阁最新面经算法竞赛18

    Hihocoder 太阁最新面经算法竞赛18 source: https://hihocoder.com/contest/hihointerview27/problems 题目1 : Big Plus ...

  7. hihoCoder太阁最新面经算法竞赛15

    hihoCoder太阁最新面经算法竞赛15 Link: http://hihocoder.com/contest/hihointerview24 题目1 : Boarding Passes 时间限制: ...

  8. vijos1144(小胖守皇宫)

    也是ural1039 描述 huyichen世子事件后,xuzhenyi成了皇上特聘的御前一品侍卫. 皇宫以午门为起点,直到后宫嫔妃们的寝宫,呈一棵树的形状:某些宫殿间可以互相望见.大内保卫森严,三步 ...

  9. zz 圣诞丨太阁所有的免费算法视频资料整理

    首发于 太阁实验室 关注专栏   写文章     圣诞丨太阁所有的免费算法视频资料整理 Ray Cao· 12 小时前 感谢大家一年以来对太阁实验室的支持,我们特地整理了在过去一年中我们所有的原创算法 ...

随机推荐

  1. 简单的Fleury算法模板

    假设数据输入时采用如下的格式进行输入:首先输入顶点个数n和边数m,然后输入每条边,每条边的数据占一行,格式为:u,v,表示从顶点u到顶点v的一条有向边 这里把欧拉回路的路径输出了出来: 手写栈: #i ...

  2. Vim command handbook

    /* 本篇文章已经默认你通过了vimtuor训练并能熟练使用大部分命令.此篇文章主要是对于tutor命令的总结和梳理.适合边学习边记忆 tutor那个完全是在学习中记忆 符合认知规律但是练习有限.所以 ...

  3. Java 实体-实体的映射框架

    一.Object mapping 的技术分类: 运行期 反射调用set/get 或者是直接对成员变量赋值 . 该方式通过invoke执行赋值 *,实现时一般会采用beanutil, Javassist ...

  4. T1503 愚蠢的宠物 codevs

    http://codevs.cn/problem/1503/  时间限制: 1 s  空间限制: 128000 KB  题目等级 : 黄金 Gold 题目描述 Description 大家都知道,sh ...

  5. mysql中如何查询最近24小时、top n查询

    MySQL中如何查询最近24小时. where visittime >= NOW() - interval 1 hour; 昨天. where visittime between CURDATE ...

  6. 【stl学习笔记】vector

    vector是定义于namespace std内的template: namespace std { template<class T, class Allocator = allocator& ...

  7. php monolog 的写日志到unix domain socket 测试终于成功

    在另外一个客户端执行 php s.php后, 通过nc -lU /tmp/tg.sck 建立的unix domain socket 有接收到消息. <?php require 'vendor/a ...

  8. EXTJS 4 树形表格组件使用演示样例

    EXTJS 4 树形表格组件使用演示样例 一.总体效果图 version=1&modificationDate=1412058826000&api=v2" alt=" ...

  9. Redis入门教程(三)— Java中操作Redis

    在Redis的官网上,我们可以看到Redis的Java客户端众多 其中,Jedis是Redis官方推荐,也是使用用户最多的Java客户端. 开始前的准备 使用jedis使用到的jedis-2.1.0. ...

  10. VC2010 利用 def 文件生成 dll 文件的方法

    近期有个需求,要生成一个dll 文件.文件里的函数都是採用 stdcall 函数调用约定,可是不希望函数名被修饰(add 被修饰成 add@8). 这时就要用def 文件了. 比方我有以下两个函数: ...