1951: [Sdoi2010]古代猪文

Time Limit: 1 Sec  Memory Limit: 64 MB
Submit: 2194  Solved: 919
[Submit][Status][Discuss]

Description

“在那山的那边海的那边有一群小肥猪。他们活泼又聪明,他们调皮又灵敏。他们自由自在生活在那绿色的大草坪,他们善良勇敢相互都关心……” ——选自猪王国民歌 很久很久以前,在山的那边海的那边的某片风水宝地曾经存在过一个猪王国。猪王国地理位置偏僻,实施的是适应当时社会的自给自足的庄园经济,很少与外界联系,商贸活动就更少了。因此也很少有其他动物知道这样一个王国。 猪王国虽然不大,但是土地肥沃,屋舍俨然。如果一定要拿什么与之相比的话,那就只能是东晋陶渊明笔下的大家想象中的桃花源了。猪王勤政爱民,猪民安居乐业,邻里和睦相处,国家秩序井然,经济欣欣向荣,社会和谐稳定。和谐的社会带给猪民们对工作火红的热情和对未来的粉色的憧憬。 小猪iPig是猪王国的一个很普通的公民。小猪今年10岁了,在大肥猪学校上小学三年级。和大多数猪一样,他不是很聪明,因此经常遇到很多或者稀奇古怪或者旁人看来轻而易举的事情令他大伤脑筋。小猪后来参加了全猪信息学奥林匹克竞赛(Pig Olympiad in Informatics, POI),取得了不错的名次,最终保送进入了猪王国大学(Pig Kingdom University, PKU)深造。 现在的小猪已经能用计算机解决简单的问题了,比如能用P++语言编写程序计算出A + B的值。这个“成就”已经成为了他津津乐道的话题。当然,不明真相的同学们也开始对他刮目相看啦~ 小猪的故事就将从此展开,伴随大家两天时间,希望大家能够喜欢小猪。 题目描述 猪王国的文明源远流长,博大精深。 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

有且仅有一行:两个数N、G,用一个空格分开。

Output

有且仅有一行:一个数,表示答案除以999911659的余数。

Sample Input

4 2

Sample Output

2048

HINT

10%的数据中,1 <= N <= 50;
20%的数据中,1 <= N <= 1000;
40%的数据中,1 <= N <= 100000;
100%的数据中,1 <= G <= 1000000000,1 <= N <= 1000000000。


1A太爽了!!!!!!

很简单啊不想写了我还要出题.......

[update 2017-02-16]

简单说一下吧,本题没必要用扩展Lucas,因为999911659分解质因子后是四个质数的乘积,直接单独用Lucas然后CRT合并就行了。唯一不好的地方是命名有点困难....

答案很显然,$G^{\sum\limits_{k|N}{N\choose k}}$

小心$G>P$导致不互质,不互质就是0啊(整除),所以g%P[0]然后快速幂的时候特判a==0

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
typedef long long ll;
const int N=4e4+5;
inline int read(){char c=getchar();int x=0,f=1;while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}return x*f;
}int n,g;
int P[5]={999911659,2,3,4679,35617};
int Pow(ll a,int b,int P){if(a==0) return 0;ll re=1;for(;b;b>>=1,a=a*a%P)if(b&1) re=re*a%P;return re;
}
int Inv(int a,int P){return Pow(a,P-2,P);}
int inv[N][5],fac[N][5],facInv[N][5];
void ini(){for(int j=1;j<=4;j++){int MOD=P[j];inv[1][j]=fac[0][j]=facInv[0][j]=1;for(int i=1;i<=MOD;i++){if(i!=1) inv[i][j]=-MOD/i*inv[MOD%i][j]%MOD;if(inv[i][j]<0) inv[i][j]+=MOD;fac[i][j]=fac[i-1][j]*i%MOD;facInv[i][j]=facInv[i-1][j]*inv[i][j]%MOD;}}
}
inline int C(int n,int m,int j){int p=P[j];return fac[n][j]*facInv[m][j]%p*facInv[n-m][j]%p;
}
int lucas(int n,int m,int j){int MOD=P[0]-1,a=1,p=P[j];for(;m;m/=p,n/=p) a=a*C(n%p,m%p,j)%p;return (ll)a*(MOD/p)%MOD*Inv(MOD/p,p)%MOD;
}
ll Lucas(ll n,ll m){ll re=0,MOD=P[0]-1;for(int i=1;i<=4;i++)re=(re+lucas(n,m,i))%MOD;return re;
}
ll solve(){int m=sqrt(n);ll re=0;for(int i=1;i<=m;i++) if(n%i==0){//printf("hi %d\n",i);re=(re+Lucas(n,i))%(P[0]-1);//printf("Lucas %d\n",Lucas(n,i));if(i*i!=n) re=(re+Lucas(n,n/i))%(P[0]-1);//printf("hi %d\n",n/i);
    }return re;
}
int gcd(int a,int b){return b==0?a:gcd(b,a%b);}
int main(){freopen("in","r",stdin);n=read();g=read();ini();//printf("solve %d\n",solve());printf("%d",Pow(g%P[0],solve(),P[0]));
}

BZOJ 1951: [Sdoi2010]古代猪文 [Lucas定理 中国剩余定理]相关推荐

  1. bzoj 1951 [Sdoi2010]古代猪文 ——数学综合

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=1951 数学综合题. 费马小定理得指数可以%999911658,又发现这个数可以质因数分解.所 ...

  2. 1951: [Sdoi2010]古代猪文

    http://www.lydsy.com/JudgeOnline/problem.php?id=1951 题意:求G^(C(N,N/K))%mod ( K|N) 1951: [Sdoi2010]古代猪 ...

  3. 【BZOJ1951】【SDOI2010】古代猪文 Lucas定理、中国剩余定理、exgcd、费马小定理

    Description "在那山的那边海的那边有一群小肥猪.他们活泼又聪明,他们调皮又灵敏.他们自由自在生活在那绿色的大草坪,他们善良勇敢相互都关心--" --选自猪王国民歌 很久 ...

  4. P2480 [SDOI2010]古代猪文 Lucas+CRT合并

    \(\color{#0066ff}{ 题目描述 }\) 猪王国的文明源远流长,博大精深. iPig在大肥猪学校图书馆中查阅资料,得知远古时期猪文文字总个数为N.当然,一种语言如果字数很多,字典也相应会 ...

  5. bzoj1951 [Sdoi2010]古代猪文 lucas+CRT+exgcd

    本来以为要套exlucas的,结果不用,码量就少了很多 注: CRT要  <=0  +=P 区分清楚加减号 码: #include<iostream> #include<cst ...

  6. P2480 [SDOI2010]古代猪文

    P2480 [SDOI2010]古代猪文 题目背景 "在那山的那边海的那边有一群小肥猪.他们活泼又聪明,他们调皮又灵敏.他们自由自在生活在那绿色的大草坪,他们善良勇敢相互都关心--" ...

  7. P2480 [SDOI2010]古代猪文(数论好题)

    P2480 [SDOI2010]古代猪文 题意: 给你n和g,求g∑d∣nCndmodpg^{\sum_{d|n}C_{n}^{d}}\bmod pg∑d∣n​Cnd​modp p=999911659 ...

  8. [bzoj1951] [Sdoi2010]古代猪文 费马小定理+Lucas定理+CRT

    Description "在那山的那边海的那边有一群小肥猪.他们活泼又聪明,他们调皮又灵敏.他们自由自在生活在那绿色的大草坪,他们善良勇敢相互都关心--" --选自猪王国民歌 很久 ...

  9. BZOJ1951 [Sdoi2010]古代猪文 【费马小定理 + Lucas定理 + 中国剩余定理 + 逆元递推 + 扩展欧几里得】...

    题目 "在那山的那边海的那边有一群小肥猪.他们活泼又聪明,他们调皮又灵敏.他们自由自在生活在那绿色的大草坪,他们善良勇敢相互都关心--" --选自猪王国民歌 很久很久以前,在山的那 ...

最新文章

  1. 手把手教你用Python模拟登录淘宝
  2. php留言板源码免mysql_PHPMYSQL留言板源码(终极完整版).doc
  3. 六十、第一个SpringBoot的 helloworld程序
  4. IOS OC项目的单例模式
  5. 深度学习之卷积神经网络(4)LeNet-5实战
  6. HDU1846 - Brave Game【巴什博弈】
  7. 了解下广告计费模式CPC、CPA和CPM
  8. 抛硬币正面期望_如果抛硬币,正面的数量多于反面的可能性
  9. 工作32:get之前打印
  10. Linux如何检查目录inode占用,linux – 如何确定哪个文件/ inode占用给定扇区
  11. python语言写九九乘法表_怎么使用Python语言写一个九九乘法表?
  12. Android笔记:触摸事件的分析与总结----TouchEvent处理机制
  13. Wireguard笔记
  14. Xshell7免费版详细下载安装教程
  15. ubuntu安装wechat
  16. 当前电子计算机发展的局限性,工程测量技术发展现状与趋势
  17. leetcode-排列组合问题
  18. 网站一直不被收录怎么办,让网站收录更好的五个方法
  19. 键入一个字母,如果小写字母输出大写,大写字母输出小写字母
  20. php 上传微信客服头像,微信多客服上传头像

热门文章

  1. Python学习笔记6(列表生成式)
  2. iOS App如何连接外设
  3. UIView的属性使用
  4. 【Daily Scrum】11-26
  5. Hadoop系列之七:分布式文件系统HDFS(2)
  6. eclipse中的感叹号和x号解决方法
  7. Datafix_for_arinvoice_dist_move
  8. mysql 运行模式_MySQL的运行模式及一些特性,引擎、事务、并发控制、优化总结...
  9. MySQL 高级repeat循环
  10. Nginx的安全控制及SSL加密介绍