BZOJ 1951 浅谈猪王国古代文字及中国剩余定理合并半拓展LuCas

世界真的很大

本来只是想找一道水题热热手，结果调了好久，听说大佬都是1A太强了

前几天做了一道拓展Lucas以为这道题也是但是其实没有那么复杂

如果是拓展Lucas那么mod数一般会给出，像这种直接告诉你的肯定有猫腻，仔细一分析就能发现

看题先：

：

“在那山的那边海的那边有一群小肥猪。他们活泼又聪明，他们调皮又灵敏。他们自由自在生活在那绿色的大草坪，他们善良勇敢相互都关心……” ——选自猪王国民歌 很久很久以前，在山的那边海的那边的某片风水宝地曾经存在过一个猪王国。猪王国地理位置偏僻，实施的是适应当时社会的自给自足的庄园经济，很少与外界联系，商贸活动就更少了。因此也很少有其他动物知道这样一个王国。 猪王国虽然不大，但是土地肥沃，屋舍俨然。如果一定要拿什么与之相比的话，那就只能是东晋陶渊明笔下的大家想象中的桃花源了。猪王勤政爱民，猪民安居乐业，邻里和睦相处，国家秩序井然，经济欣欣向荣，社会和谐稳定。和谐的社会带给猪民们对工作火红的热情和对未来的粉色的憧憬。 小猪iPig是猪王国的一个很普通的公民。小猪今年10岁了，在大肥猪学校上小学三年级。和大多数猪一样，他不是很聪明，因此经常遇到很多或者稀奇古怪或者旁人看来轻而易举的事情令他大伤脑筋。小猪后来参加了全猪信息学奥林匹克竞赛(Pig in , POI)，取得了不错的名次，最终保送进入了猪王国大学(Pig , 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觉得这个数字可能是天文数字，所以你只需要告诉他答案除以的余数就可以了。

input：

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

：

有且仅有一行：一个数，表示答案除以的余数。

题目巨长系列

说白了就是枚举n的所有因数，算一个C（n，k）加起来，在算G的这么多次方

考虑前面算的组合数多的位置都是在G的指数位置，所以说可以用欧拉定理来优化一下，指数mod一个phi（m）

那么求组合数时的mod数就是这个phi了

问题在于是一个质数，那他的phi就是，恰恰是一个合数，对合数求组合数，这不是拓展Lucas吗？

考虑这个，作为题目给出的mod数却不是1e9+7之类的司空见惯的数，有点怪

-1得到

这个东西质因数一分解，就是 2,3,4679,35617，全是质数，对于每一个质数可以使用Lucas，算出来5个值再用中国剩余定理合起来就好啦

然而我一开始先是打表没有打出35617，调了好久才发现

之后又是Lucas的时候mod写成了而不是mi，orz

两个小时就这么过去了

完整代码：

#include
#include
using namespace std;
typedef long long dnt ;const dnt mod=999911659;dnt n,G;
dnt saber[5][100050],inv[5][100050],m[5]={0,2,3,4679,35617};void init(dnt m,int typ)
{saber[typ][0]=inv[typ][0]=inv[typ][1]=1;for(int i=1;i<=m;i++) saber[typ][i]=saber[typ][i-1]*i%m;for(int i=2;i<=m;i++) inv[typ][i]=(m-m/i)*inv[typ][m%i]%m;for(int i=1;i<=m;i++) inv[typ][i]=inv[typ][i-1]*inv[typ][i]%m;
}void exgcd(dnt a,dnt b,dnt &x,dnt &y)
{if(b==0){x=1,y=0;return ;}dnt x0,y0;exgcd(b,a%b,x0,y0);x=y0,y=x0-a/b*y0;
}dnt inverse(dnt a,dnt m)
{dnt x,y;exgcd(a,m,x,y);return (x%m+m)%m;
}dnt Misaka(dnt a,dnt b,dnt m,int typ)
{if(areturn 0;if(a<m) return saber[typ][a]*inv[typ][b]%m*inv[typ][a-b]%m;return Misaka(a%m,b%m,m,typ)*Misaka(a/m,b/m,m,typ)%m;
}dnt CRT(dnt a,dnt b)
{dnt rt=0,Mod=mod-1;for(int i=1;i<=4;i++){dnt tmp=Misaka(a,b,m[i],i);rt=(rt+(Mod/m[i])*tmp%Mod*inverse(Mod/m[i],m[i])%Mod)%Mod;}   return rt%Mod;
}dnt cal()
{dnt rt=0;for(dnt i=1;i*i<=n;i++){if(n%i) continue ;rt=(rt+CRT(n,i))%(mod-1);if(i!=n/i) rt=(rt+CRT(n,n/i))%(mod-1);}return rt%(mod-1);
}dnt quickmub(dnt a,dnt b)
{dnt rt=1;while(b){if(b&1) rt=rt*a%mod;a=a*a%mod,b>>=1;}return rt;
}int main()
{for(int i=1;i<=4;i++)init(m[i],i);scanf("%lld%lld",&n,&G);dnt idx=cal();printf("%lld\n",quickmub(G,idx+mod-1));return 0;
}
/*
EL PSY CONGROO
*/

嗯，就是这样