【BZOJ1951】【SDOI2010】古代猪文
Background
“在那山的那边海的那边有一群小肥猪。他们活泼又聪明,他们调皮又灵敏。他们自由自在生活在那绿色的大草坪,他们善良勇敢相互都关心……”
——选自猪王国民歌
很久很久以前,在山的那边海的那边的某片风水宝地曾经存在过一个猪王国。猪王国地理位置偏僻,实施的是适应当时社会的自给自足的庄园经济,很少与外界联系,商贸活动就更少了。因此也很少有其他动物知道这样一个王国。
猪王国虽然不大,但是土地肥沃,屋舍俨然。如果一定要拿什么与之相比的话,那就只能是东晋陶渊明笔下的大家想象中的桃花源了。猪王勤政爱民,猪民安居乐业,邻里和睦相处,国家秩序井然,经济欣欣向荣,社会和谐稳定。和谐的社会带给猪民们对工作火红的热情和对未来的粉色的憧憬。
小猪iPig是猪王国的一个很普通的公民。小猪今年10岁了,在大肥猪学校上小学三年级。和大多数猪一样,他不是很聪明,因此经常遇到很多或者稀奇古怪或者旁人看来轻而易举的事情令他大伤脑筋。小猪后来参加了全猪信息学奥林匹克竞赛(Pig Olympiad in Informatics, POI),取得了不错的名次,最终保送进入了猪王国大学(Pig Kingdom University, PKU)深造。
现在的小猪已经能用计算机解决简单的问题了,比如能用P++语言编写程序计算出A + B的值。这个“成就”已经成为了他津津乐道的话题。当然,不明真相的同学们也开始对他刮目相看啦~
小猪的故事就将从此展开,伴随大家两天时间,希望大家能够喜欢小猪。
Description
猪王国的文明源远流长,博大精深。
iPig在大肥猪学校图书馆中查阅资料,得知远古时期猪文文字总个数为N。当然,一种语言如果字数很多,字典也相应会很大。当时的猪王国国王考虑到如果修一本字典,规模有可能远远超过康熙字典,花费的猪力、物力将难以估量。故考虑再三没有进行这一项劳猪伤财之举。当然,猪王国的文字后来随着历史变迁逐渐进行了简化,去掉了一些不常用的字。
iPig打算研究古时某个朝代的猪文文字。根据相关文献记载,那个朝代流传的猪文文字恰好为远古时期的k分之一,其中k是N的一个正约数(可以是1和N)。不过具体是哪k分之一,以及k是多少,由于历史过于久远,已经无从考证了。
iPig觉得只要符合文献,每一种能整除N的k都是有可能的。他打算考虑到所有可能的k。显然当k等于某个定值时,该朝的猪文文字个数为N / k。然而从N个文字中保留下N / k个的情况也是相当多的。iPig预计,如果所有可能的k的所有情况数加起来为P的话,那么他研究古代文字的代价将会是G的P次方。
现在他想知道猪王国研究古代文字的代价是多少。由于iPig觉得这个数字可能是天文数字,所以你只需要告诉他答案除以999911659的余数就可以了。
Input
输入文件ancient.in有且仅有一行:两个数N、G,用一个空格分开。
Output
输出文件ancient.out有且仅有一行:一个数,表示答案除以999911659的余数。
Sample Input
4 2
Sample Output
2048
Hint
10%的数据中,\(1 \leq N \leq 50\);
20%的数据中,\(1 \leq N \leq 1000\);
40%的数据中,\(1 \leq N \leq 100000\);
100%的数据中,\(1 \leq N \leq 10^9\)。
Solution
题目要求求出$$g^{\Sigma_{k|N} C_{N}^{N/k}}$$
首先,根据欧拉定理,有\(a^{\varphi(p)} \equiv 1\ (mod\ p)\)
又999911659为质数,故只需求上述质数对\({\varphi(p)}\)取模的结果
但是999911658不是质数,要怎么办呢,CRT!!!
\(999911658=2*3*4679*35617\),对这4个质数分别求出\(\Sigma_{k|N} C_{N}^{N/k}\)
取模后的值,然后利用CRT合并即可。现在的问题只剩下,如何计算\(C_{N}^{N/k}\) 了,显然N太大无法递推,我们考虑利用Lucas定理计算即可。
如此下来时间复杂度为\(O(\sqrt n log n)\).
由于之前blog都没提到CRT和Lucas 定理,这里简单介绍一下好了。
CRT
CRT说的是这么一件事:对于给定的同余方程组\(x\equiv a_i\ (mod\ m_i)\)若\(m_i\)两两互素,则该方程存在整数解,且在模\(M=\prod_{i=1}^{n} m_i\) 的意义下存在唯一解。
解为\(x \equiv \Sigma_{i=1}^{n} a_i \frac{M}{m_i} (\frac{M}{m_i})^{-1} \ (\ mod\ M)\)
根据定理去算即可,CRT的时间复杂度为\(O(n\log max(m_i))\).
Lucas 定理
Lucas 定理是这样的:对于一个质数p,\(C(n,m) \equiv C(n/p,m/p)*C(n\ mod\ p,m\ mod\ p) \ (mod\ p)\)
用于大质数取模,对于\(n,m<p\) 的情况,利用阶乘以及其乘法逆元计算,可以用O(p)的时间先预处理出,时间复杂度为\(O(\log_p n+p)\)
Code
#include <stdio.h>
#define MN 40000
#define R register
#define ll long long
#define P 999911658
#define file(x) freopen(#x".in","r",stdin),freopen(#x".out","w",stdout)
#define endfile fclose(stdin),fclose(stdout)
inline int read(){
R int x; R bool f; R char c;
for (f=0; (c=getchar())<'0'||c>'9'; f=c=='-');
for (x=c-'0'; (c=getchar())>='0'&&c<='9'; x=(x<<3)+(x<<1)+c-'0');
return f?-x:x;
}
const int pr[]={2,3,4679,35617};
int inv[MN],fc[MN],g,n;
inline int C(int n,int m,int p){
if (n<m) return 0;
if (n<p&&m<p) return (ll)fc[n]*inv[m]%p*inv[n-m]%p;
return (ll)C(n/p,m/p,p)*C(n%p,m%p,p)%p;
}
inline int pw(int x,int k,int p){
if (!x) return 0;R int res=1;
for (; k; k>>=1,x=(ll)x*x%p) if (k&1) res=(ll)res*x%p;
return res;
}
inline int solve(int p){
fc[0]=1;for (R int i=1; i<p; ++i) fc[i]=(ll)fc[i-1]*i%p;
inv[0]=inv[1]=1;for (R int i=2; i<p; ++i) inv[i]=(ll)(p-p/i)*inv[p%i]%p;
for (R int i=1; i<p; ++i) inv[i]=(ll)inv[i]*inv[i-1]%p;
R int res=0;
for (R int i=1; i*i<=n; ++i)
if (n%i==0){
if (i*i==n) res=(res+C(n,i,p))%p;
else res=(res+C(n,i,p)+C(n,n/i,p))%p;
}return res;
}
inline int CRT(){
R int res=0;
for (R int i=0; i<4; ++i){
R int tmp=solve(pr[i]);
tmp=(ll)pw(P/pr[i]%pr[i],pr[i]-2,pr[i])*(P/pr[i])%P*tmp%P;
res+=tmp;if (res>=P) res-=P;
}return res;
}
int main(){n=read(),g=read();printf("%d\n",pw(g%(P+1),CRT(),P+1));}
【BZOJ1951】【SDOI2010】古代猪文的更多相关文章
- [bzoj1951] [Sdoi2010]古代猪文 费马小定理+Lucas定理+CRT
Description "在那山的那边海的那边有一群小肥猪.他们活泼又聪明,他们调皮又灵敏.他们自由自在生活在那绿色的大草坪,他们善良勇敢相互都关心--" --选自猪王国民歌 很久 ...
- BZOJ1951[SDOI2010]古代猪文
Description "在那山的那边海的那边有一群小肥猪.他们活泼又聪明,他们调皮又灵敏.他们自由自在生活在那绿色的大草坪,他们善良勇敢相互都关心--" --选自猪王国民歌 很久 ...
- BZOJ1951:[SDOI2010]古代猪文(Lucas,CRT)
Description “在那山的那边海的那边有一群小肥猪.他们活泼又聪明,他们调皮又灵敏.他们自由自在生活在那绿色的大草坪,他们善良勇敢相互都关心……” ——选自猪王国民歌 很久很久以前,在山的那边 ...
- BZOJ1951 [Sdoi2010]古代猪文 【费马小定理 + Lucas定理 + 中国剩余定理 + 逆元递推 + 扩展欧几里得】
题目 "在那山的那边海的那边有一群小肥猪.他们活泼又聪明,他们调皮又灵敏.他们自由自在生活在那绿色的大草坪,他们善良勇敢相互都关心--" --选自猪王国民歌 很久很久以前,在山的那 ...
- BZOJ1951 [Sdoi2010]古代猪文 中国剩余定理 快速幂 数论
原文链接http://www.cnblogs.com/zhouzhendong/p/8109156.html 题目传送门 - BZOJ1951 题意概括 求 GM mod 999911659 M=∑i ...
- bzoj千题计划323:bzoj1951: [Sdoi2010]古代猪文(Lucas+CRT+欧拉定理)
https://www.lydsy.com/JudgeOnline/problem.php?id=1951 先欧拉降幂 然后模数质因数分解 分别计算组合数的结果,中国剩余定理合并 #include&l ...
- bzoj1951 [Sdoi2010]古代猪文 ——数论综合
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=1951 题意就是要求 G^( ∑(k|n) C(n,k) ) % p,用费马小定理处理指数,卢 ...
- 【BZOJ1951】[SDOI2010]古代猪文
[BZOJ1951][SDOI2010]古代猪文 题面 bzoj 洛谷 题解 题目实际上是要求 $ G^{\sum d|n\;C_n^d}\;mod \; 999911659 $ 而这个奇怪的模数实际 ...
- 【BZOJ1951】[Sdoi2010]古代猪文 Lucas定理+CRT
[BZOJ1951][Sdoi2010]古代猪文 Description 求$X=\sum\limits_{d|n}C_n^d$,$Ans=G^X (\mod 999911659)$. Input 有 ...
- 【bzoj1951】: [Sdoi2010]古代猪文 数论-中国剩余定理-Lucas定理
[bzoj1951]: [Sdoi2010]古代猪文 因为999911659是个素数 欧拉定理得 然后指数上中国剩余定理 然后分别lucas定理就好了 注意G==P的时候的特判 /* http://w ...
随机推荐
- Beta敏捷冲刺每日报告——Day1
1.情况简述 Beta阶段Scrum Meeting 敏捷开发起止时间 2017.11.2 00:00 -- 2017.11.3 00:00 讨论时间地点 2017.11.2 晚9:30,电话会议会议 ...
- NFC驱动调试
1.NFC基本概念: NFC 又称为近场通信,是一种新兴技术,可以在彼此靠近的情况下进行数据交换,是由非接触式射频识别(RFID) 及互连互通技术整合演变而来,通过单一芯片集成感应式读卡器: NFC有 ...
- JAVA_SE基础——47.接口
如果一个抽象类中的所有方法都是抽象的,则可以将这个类用另一种方法来定义,即接口~ 在定义接口时,需要用interface关键字来声明,具体实例如code1 接口的定义格式:interface 接口名{ ...
- machine learning 之 导论 一元线性回归
整理自Andrew Ng 的 machine learnig 课程 week1. 目录: 什么是机器学习 监督学习 非监督学习 一元线性回归 模型表示 损失函数 梯度下降算法 1.什么是机器学习 Ar ...
- volt问题
1./表示当前目录:/college/detail/{{ item.sid }}表示这个路径超链接,url实在不好写就不写,作为开发人员想怎么弄就怎么弄最后发布是项目主管的事 2.不需要服务器给值,直 ...
- 求解数独难题, Sudoku问题(回溯)
Introduction : 标准的数独游戏是在一个 9 X 9 的棋盘上填写 1 – 9 这 9 个数字,规则是这样的: 棋盘分成上图所示的 9 个区域(不同颜色做背景标出,每个区域是 3 X 3 ...
- Appium+python测试app实例
Appium和selenium差不到,只是一个用于测web,一个用于测APP.下面记录一下我搭的测试框架,同样是基于PO模式,用的unittest. 最后测试报告如下: 1.1 代码结构 这 ...
- Oracle 使用pl/sql将表中的数据读出到文件中
(1)在服务器上创建文件路径及文件 [oracle@redhat errormsg]$ touch test01.txt (2)在数据库中创建路径及授权 (3)创建存储过程 CREATE OR REP ...
- kubernetes入门(07)kubernetes的核心概念(4)
一.pod 二.Volume volume可以为容器提供持久化存储,比如 三.私有镜像 在使用私有镜像时,需要创建一个docker registry secret,并在容器中引用.创建docker r ...
- 新概念英语(1-25)Mrs. Smith's Kitchen
What colour is the electric cooker? Mrs. Smith's Kitchen is small. There is a refrigerator in the ki ...