\(\color{#0066ff}{ 题目描述 }\)

猪王国的文明源远流长,博大精深。

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的余数就可以了。

\(\color{#0066ff}{输入格式}\)

输入文件ancient.in有且仅有一行:两个数N、G,用一个空格分开。

\(\color{#0066ff}{输出格式}\)

输出文件ancient.out有且仅有一行:一个数,表示答案除以999911659的余数。

\(\color{#0066ff}{输入样例}\)

4 2

\(\color{#0066ff}{输出样例}\)

2048

\(\color{#0066ff}{数据范围与提示}\)

10%的数据中,1 <= N <= 50;

20%的数据中,1 <= N <= 1000;

40%的数据中,1 <= N <= 100000;

100%的数据中,1 <= G <= 1000000000,1 <= N <= 1000000000。

\(\color{#0066ff}{ 题解 }\)

一句话题意,求\(G^{\sum_{k|n}C_n^k} \mod 999911659\)

根据欧拉定理

$G{\sum_{k|n}C_nk} \equiv G{\sum_{k|n}C_nk \mod 999911658}\mod 999911659 $

现在我们要求\(\sum_{k|n}C_n^k \mod 999911658\)

把\(999911658\)质因数分解

\(999911658=2*3*4679*35617\)

分别算出上面组合数对四个模数的ans,然后用CRT合并

最后快速幂出答案

值的注意的是

1、阶乘和逆元要预处理

2、要考虑不互质的情况,用扩展欧拉定理

#include<bits/stdc++.h>
#define LL long long
LL in() {
char ch; LL x = 0, f = 1;
while(!isdigit(ch = getchar()))(ch == '-') && (f = -f);
for(x = ch ^ 48; isdigit(ch = getchar()); x = (x << 1) + (x << 3) + (ch ^ 48));
return x * f;
}
bool flag;
const int mix = 999911659;
const int mox = 999911658;
LL b[] = {0LL, 2LL, 3LL, 4679LL, 35617LL};
LL a[10];
LL fac[50505][5], inv[50505][5];
LL ksm(LL x, LL y, LL mod) {
if(x == 0) return 0;
LL re = 1LL;
while(y) {
if(y & 1) re = re * x % mod;
x = x * x % mod;
y >>= 1;
}
return re % mod;
}
LL C(LL n, LL m, LL mod, int pos) {
if(m > n) return 0;
if(n >= mod || m >= mod) return C(n / mod, m / mod, mod, pos) * C(n % mod, m % mod, mod, pos) % mod;
return ((fac[n][pos] * inv[m][pos] % mod) * inv[n - m][pos]) % mod;
}
void exgcd(LL a, LL b, LL &x, LL &y) {
if(!b) return(void)(x = 1, y = 0);
exgcd(b, a % b, x, y);
LL t = x - a / b * y;
x = y, y = t;
}
LL work(LL n, LL k) {
LL ans = 0;
for(int i = 1; i <= 4; i++) {
a[i] = C(n, k, b[i], i);
LL x, y;
exgcd(b[i], mox / b[i], x, y);
x = ((mox / b[i] * y) % mox + mox) % mox;
(ans += (a[i] * x % mox)) %= mox;
}
return ((ans % mox) + mox) % mox;
}
void predoit() {
for(int i = 1; i <= 4; i++) fac[0][i] = 1;
for(int i = 1; i <= 4; i++)
for(int j = 1; j < b[i]; j++)
fac[j][i] = 1LL * fac[j - 1][i] * j % b[i];
for(int i = 1; i <= 4; i++) inv[b[i] - 1][i] = ksm(fac[b[i] - 1][i], b[i] - 2, b[i]);
for(int i = 1; i <= 4; i++)
for(int j = b[i] - 2; j >= 0; j--)
inv[j][i] = 1LL * (j + 1) * inv[j + 1][i] % b[i];
} int main() {
LL n = in(), g = in();
LL ans = 0;
predoit();
for(LL i = 1; i * i <= n; i++) {
if(n % i == 0) {
ans = ans + work(n, i);
if(ans >= mox) ans %= mox, flag = true;
if(i * i != n) {
ans = ans + work(n, n / i);
if(ans >= mox) ans %= mox, flag = true;
}
}
}
printf("%lld\n", ksm(g, (flag? ans + mox : ans), mix));
return 0;
}

P2480 [SDOI2010]古代猪文 Lucas+CRT合并的更多相关文章

  1. BZOJ1951:[SDOI2010]古代猪文(Lucas,CRT)

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

  2. 洛咕 P2480 [SDOI2010]古代猪文

    洛咕 P2480 [SDOI2010]古代猪文 题目是要求\(G^{\sum_{d|n}C^d_n}\). 用费马小定理\(G^{\sum_{d|n}C^d_n\text{mod 999911658} ...

  3. 洛谷 P2480 [SDOI2010]古代猪文 解题报告

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

  4. P2480 [SDOI2010]古代猪文

    P2480 [SDOI2010]古代猪文 比较综合的一题 前置:Lucas 定理,crt 求的是: \[g^x\bmod 999911659,\text{其中}x=\sum_{d\mid n}\tbi ...

  5. 【题解】P2480 [SDOI2010]古代猪文 - 卢卡斯定理 - 中国剩余定理

    P2480 [SDOI2010]古代猪文 声明:本博客所有题解都参照了网络资料或其他博客,仅为博主想加深理解而写,如有疑问欢迎与博主讨论✧。٩(ˊᗜˋ)و✧*。 题目描述 猪王国的文明源远流长,博大精 ...

  6. 【BZOJ1951】[Sdoi2010]古代猪文 Lucas定理+CRT

    [BZOJ1951][Sdoi2010]古代猪文 Description 求$X=\sum\limits_{d|n}C_n^d$,$Ans=G^X (\mod 999911659)$. Input 有 ...

  7. BZOJ 1951: [Sdoi2010]古代猪文 [Lucas定理 中国剩余定理]

    1951: [Sdoi2010]古代猪文 Time Limit: 1 Sec  Memory Limit: 64 MBSubmit: 2194  Solved: 919[Submit][Status] ...

  8. 洛谷 P2480 [SDOI2010]古代猪文 题解【欧拉定理】【CRT】【Lucas定理】

    数论综合题. 题目背景 题目背景与题目无关因此省略.题目链接 题目描述 猪王国的文明源远流长,博大精深. iPig 在大肥猪学校图书馆中查阅资料,得知远古时期猪文文字总个数为 \(N\).当然,一种语 ...

  9. Luogu P2480 [SDOI2010]古代猪文 卢卡斯+组合+CRT

    好吧刚开始以为扩展卢卡斯然后就往上套..结果奇奇怪怪又WA又T...后来才意识到它的因子都是质数...qwq怕不是这就是学知识学傻了.. 题意:$ G^{\Sigma_{d|n} \space C_n ...

随机推荐

  1. Delphi XE2 新控件 布局Panel TGridPanel TFlowPanel

    Delphi XE2 新控件 Firemonkey 布局Panel Windows平台VCl TGridPanel TFlowPanel FMX 跨平台 TLayout TGridLayout TFl ...

  2. vsftp部署和优化错误

    ftp登录失败 vim /etc/vsftpd/vsftpd.conf 添加虚拟机配置的时候有空行,删除空行解决

  3. SQL 实现行列互换

    Oracle:不过大多数是采用 oracle 数据库当中的一些便捷函数进行处理,比如 ”pivot”: MySql:目前没有找到更好的方法 题目:数据库中有一张如下所示的表,表名为sales. 年 季 ...

  4. 基于cookie实现用户验证

    #!/usr/bin/env python import tornado.ioloop import tornado.web class IndexHander(tornado.web.Request ...

  5. 使用网络用户命令行工具的/passwordreq:yes

    提示:"新建域时,本地administrator帐户将成为域administrator账户.无法新建域,因为本地administrator账户密码不符合要求.目前,本地administrat ...

  6. 1-1+zookeeper简介

     zookeeper是中间件,可以为分布式系统提供协调服务.如果没有分布式系统,zookeeper也发挥不了它的优势.

  7. webfrom 做项目的注意事项

    1.展示细节 如男女显示问题 不能显示true   false  时间转换成中文  民族显示汉字  不能直接显示代码2.用户名重复验证 从数据库中查询验证4.日期判断 判断年份有点问题   var y ...

  8. schedule与scheduleAtFixedRate比较

    schedule与scheduleAtFixedRate: 不延时: schedule(TimerTask, Date runDate, long period)方法任务不延时----Date类型 i ...

  9. 【摘自lvs官网】lvs介绍

    Linux Virtual Server项目的目标 :使用集群技术和Linux操作系统实现一个高性能.高可用的服务器,它具有很好的可伸缩性(Scalability).可靠性(Reliability)和 ...

  10. 关于const指针的误解

    转自 http://blog.csdn.net/yingxunren/article/details/3968800 const char*, char const*, char*const的区别问题 ...