luogu P1587 [NOI2016]循环之美
首先要知道什么样的数才是"纯循环数".打表可以发现,这样的数当且仅当分母和\(k\)互质,这是因为,首先考虑除法过程,每次先给当前余数\(*k\),然后对分母做带余除法,那么出现循环就要使的某一次除完后的余数在前面出现过.并且有欧拉定理\(a^{\varphi(n)}\equiv 1 (\mod n)(\gcd(a,n)=1)\),这样可以使得在计算小数点后一位时的余数在若干次后再次出现
然后要使得数值不同,所以其实要求的是这个东西$$\sum_{i=1}{n}\sum_{j=1}{m}[\gcd(i,j)=1][\gcd(j,k)=1]$$
先把只有\(j\)和\(k\)的提前$$\sum_{j=1}{m}[\gcd(j,k)=1]\sum_{i=1}{n}[\gcd(i,j)=1]$$
然后把\([\gcd(j,k)=1]\)转化一下 $$\sum_{j=1}^{m}
\sum_{d|j,d|k}\mu(d)\sum_{i=1}^{n}[\gcd(i,j)=1]$$
把\(d\)提前$$\sum_{d=1}{m}\mu(d)\sum_{j=1}{\lfloor\frac{m}{d}\rfloor}
\sum_{i=1}{n}[\gcd(i,jd)=1]$$$$\sum_{d=1}{m}\mu(d)\sum_{j=1}^{\lfloor\frac{m}{d}\rfloor}
\sum_{i=1}^{n}[\gcd(i,j)=1][\gcd(i,d)=1]$$
后面那个是不是有点眼熟?我们如果记\(s(i,j,k)=\sum_{i=1}^{n}\sum_{j=1}^{m}[\gcd(i,j)=1][\gcd(j,k)=1]\),那么就能得到$$s(i,j,k)=\sum_{d|k}\mu(d)s(\lfloor\frac{m}{d}\rfloor,n,d)$$
这个递归处理就好了,边界就是\(m=0\)或者\(n=0\)时为\(0\),\(k=1\)时为\(\sum_{i=1}^{n}\sum_{j=1}^{m}[\gcd(i,j)=1]\),数论分块+杜教筛即可.复杂度大概是\(O(\log n\log m \sqrt{n}+n^{\frac{2}{3}})\),这个还是比较慢的,加了记忆化都要1800多ms 也可能是我写丑了
// luogu-judger-enable-o2
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<vector>
#include<cmath>
#include<ctime>
#include<queue>
#include<map>
#include<set>
#define LL long long
#define db double
using namespace std;
const int N=5e6+10,M=2000+10;
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;
}
int n,m,k,prm[N],tt;
bool pp[N];
LL mu[N];
map<int,LL> f;
int lm;
LL siv(int nn)
{
if(nn<=N-10) return mu[nn];
if(f.count(nn)) return f[nn];
LL &an=f[nn];
an=1;
for(int i=2,j;i<=nn;i=j+1)
{
j=nn/(nn/i);
an-=1ll*siv(nn/i)*(j-i+1);
}
return an;
}
LL ff(int nn,int mm)
{
lm=min(nn,mm);
LL an=0;
for(int i=1,j;i<=lm;i=j+1)
{
j=min(nn/(nn/i),mm/(mm/i));
an+=1ll*(siv(j)-siv(i-1))*(nn/i)*(mm/i);
}
return an;
}
vector<int> dd[M];
struct node
{
int n,m,k;
bool operator < (const node &bb) const {return n!=bb.n?n<bb.n:(m!=bb.m?m<bb.m:k<bb.k);}
};
map<node,LL> g;
LL sov(int n,int m,int k)
{
if(n<=0||m<=0) return 0;
node nw=(node){n,m,k};
if(g.count(nw)) return g[nw];
if(k==1) return g[nw]=ff(n,m);
LL an=0;
vector<int>::iterator it;
for(it=dd[k].begin();it!=dd[k].end();++it)
{
int i=*it;
an+=1ll*sov(m/i,n,i)*(mu[i]-mu[i-1]);
}
return g[nw]=an;
}
int main()
{
mu[1]=1;
for(int i=2;i<=N-10;++i)
{
if(!pp[i]) prm[++tt]=i,mu[i]=-1;
for(int j=1;j<=tt&&i*prm[j]<=N-10;++j)
{
pp[i*prm[j]]=1;
if(i%prm[j]==0) break;
mu[i*prm[j]]=-mu[i];
}
}
for(int i=2;i<=N-10;++i) mu[i]+=mu[i-1];
n=rd(),m=rd(),k=rd();
for(int i=1;i<=k;++i)
if(k%i==0) dd[k].push_back(i);
int nn=dd[k].size();
for(int i=0;i<nn-1;++i)
{
int x=dd[k][i];
vector<int>::iterator it;
for(it=dd[k].begin();it!=dd[k].end();++it)
{
int y=*it;
if(x<y) break;
if(x%y==0) dd[x].push_back(y);
}
}
printf("%lld\n",sov(n,m,k));
return 0;
}
luogu P1587 [NOI2016]循环之美的更多相关文章
- luogu 1587 [NOI2016]循环之美
LINK:NOI2016循环之美 这道题是 给出n m k 求出\(1\leq i\leq n,1\leq j\leq m\) \(\frac{i}{j}\)在k进制下是一个纯循环的. 由于数值相同的 ...
- 洛谷P1587 [NOI2016]循环之美
传送门 不会,先坑着 https://kelin.blog.luogu.org/solution-p1587 //minamoto #include<cstdio> #include< ...
- 并不对劲的bzoj4652:loj2085:uoj221:p1587:[NOI2016]循环之美
题目大意 对于已知的十进制数\(n\)和\(m\),在\(k\)进制下,有多少个数值上互不相等的纯循环小数,可以用\(x/y\)表示,其中 \(1\leq x\leq n,1\leq y\leq m\ ...
- [UOJ#221][BZOJ4652][Noi2016]循环之美
[UOJ#221][BZOJ4652][Noi2016]循环之美 试题描述 牛牛是一个热爱算法设计的高中生.在他设计的算法中,常常会使用带小数的数进行计算.牛牛认为,如果在 k 进制下,一个数的小数部 ...
- 题解 P1587 【[NOI2016]循环之美】
知识点:莫比乌斯反演 积性函数 杜教筛 废话前言: 我是古明地恋,写这篇题解的人已经被我 请各位读者自行无视搞事的恋恋带有删除线的内容,谢谢茄子. 这道题目本身并不难,但是公式推导/代码过程中具有迷惑 ...
- bzoj4652 [Noi2016]循环之美
Description 牛牛是一个热爱算法设计的高中生.在他设计的算法中,常常会使用带小数的数进行计算.牛牛认为,如果在k进制下,一个数的小数部分是纯循环的,那么它就是美的.现在,牛牛想知道:对于已知 ...
- [NOI2016]循环之美
Description 牛牛是一个热爱算法设计的高中生.在他设计的算法中,常常会使用带小数的数进行计算.牛牛认为,如果在 k 进制下,一个数的小数部分是纯循环的,那么它就是美的.现在,牛牛想知道:对 ...
- BZOJ4652: [Noi2016]循环之美(莫比乌斯反演,杜教筛)
Description 牛牛是一个热爱算法设计的高中生.在他设计的算法中,常常会使用带小数的数进行计算.牛牛认为,如果在 k 进制下,一个数的小数部分是纯循环的,那么它就是美的.现在,牛牛想知道:对 ...
- BZOJ4652 [Noi2016]循环之美 【数论 + 莫比乌斯反演 + 杜教筛】
题目链接 BZOJ 题解 orz 此题太优美了 我们令\(\frac{x}{y}\)为最简分数,则\(x \perp y\)即,\(gcd(x,y) = 1\) 先不管\(k\)进制,我们知道\(10 ...
随机推荐
- ...扩展运算符+rest参数+call/apply/bind
之前在set,map里面有提过扩展运算符的概念,但是今天偶然遇到一个问题,类似于扩展运算符的经典用法,突然发现对其了解不是很深,所以再来整理一下扩展运算符的相关知识. 重点:扩展运算符内部调用的是数据 ...
- linux安装提示权限 不足 ,原来是锁了
lsattr 查看 然后就是chattr -i /etc -i 是解锁 +i 就是上锁
- C++入门经典-例6.2-将二维数组进行行列对换
1:一维数组的初始化有两种,一种是单个逐一赋值,一种是使用聚合方式赋值.聚合方式的例子如下: int a[3]={1,2,3}; int a[]={1,2,3};//编译器能够获得数组元素的个数 in ...
- Watir单元库
http://www.cnblogs.com/Javame/p/4045229.html test: #require 'net/http' #require 'uri' #url = URI.par ...
- 后盾网lavarel视频项目---lavarel多表关联一对多操作实例
后盾网lavarel视频项目---lavarel多表关联一对多操作实例 一.总结 一句话总结: 1.一对多中多那个部分的数据前端通过json弄到服务器 2.所有通过一操作多的时候,都要用上模型中定义的 ...
- 【转】JS正则验证邮手机、箱等格式
function test() { var temp = document.getElementById("text1"); //对电子邮件的验证 var myreg = /^([ ...
- Libvirt 版本降级过程记录 4.5.0 to 3.9.0
前言 搞 OpenStack 开发 Libvirt 版本会随着 OpenStack 版本切来切去的,记录一下 Libvirt 从 4.5 降级到 3.9.0 的过程. 步骤 直接 downgrade ...
- Selenium 2自动化测试实战2(数组与字典)
一.数组与字典 1.数组 数组用方括号([])表示,里面的每一项用逗号(,)隔开 Prthon允许在数组里面任意地放置数字或字符串.需要注意的是,数组下标是从0开始的,所以,lists[0]会输出数组 ...
- JavaScript编程精解 初读笔记
1.1 值 JavaScript里有6种基本类型的值: number.string.Boolean.object.function和undefined. 1.3.2 prompt和confirm 浏览 ...
- python学习之模块-模块(五)
5.10 包 5.10.1 包的概念 [官网解释] Packages are a way of structuring Python's module namespace by using " ...