gym/100633/J Ceizenpok’s formula

题意：给出n,m,p输出C(n,m)%p p不一定是素数

#include <iostream>//ex_Lucas模板
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <map>
//#define IO  ios::sync_with_stdio(false),cin.tie(0), cout.tie(0);
//#pragma comment(linker, "/STACK:1024000000,1024000000")
typedef long long ll;
void ex_gcd(ll a, ll b, ll &d, ll &x, ll &y)
{ if (!b) { x = 1; y = 0; d = a; } else { ex_gcd(b, a%b, d, y, x); y -= x * (a / b); }; }
int gcd(int a, int b) { return b ? gcd(b, a%b) : a; }
int lcm(int a,int b){return a/gcd(a,b)*b;}//ÏÈ³ýºó³Ë·ÀÒç³ö
//int inv_exgcd(int a, int m) { int d, x, y;ex_gcd(a, m, d, x, y);return d == 1 ? (x + m) % m : -1; }
ll reverse(ll a, ll m) { ll d, x, y;ex_gcd(a, m, d, x, y);return d == 1 ? (x + m) % m : -1; }
typedef long long ll;
const int maxn=1e5;
using namespace std;
ll cnt=0;
int a[maxn];
int primer[maxn];
void isprime()
{cnt=0;for(int i=2;i<maxn;++i){if(!a[i]){primer[cnt++]=i;for(int j=i*2;j<maxn;j+=i)a[j]=1;}}
}
ll Pow(ll a,ll n,ll mod)
{ll ans=1;while(n){if(n&1)ans=ans*a%mod;a=a*a%mod;n>>=1;}return ans;
}
ll C(ll n,ll p,ll pk)
{if(n==0)return 1;ll ans=1;for(ll i=2;i<=pk;++i)if(i%p)ans=ans*i%pk;ans=Pow(ans,n/pk,pk);for(ll k=n%pk,i=2;i<=k;++i)if(i%p)ans=ans*i%pk;return ans*C(n/p,p,pk)%pk;
}
ll ex_lucas(ll n,ll m,ll p,ll pi,ll pk) //C(n,m)%p (p非素数) pi:p分解的一个素因子,pk：pi素因子对应的幂次
{ll i,j,k=0,a,b,c,ans;a=C(n,pi,pk),b=C(m,pi,pk),c=C(n-m,pi,pk);for(i=n;i;i/=pi)k+=i/pi;for(i=m;i;i/=pi)k-=i/pi;for(i=n-m;i;i/=pi)k-=i/pi;ans=a*reverse(b,pk)%pk*reverse(c,pk)%pk*Pow(pi,k,pk)%pk;return ans*(p/pk)%p*reverse(p/pk,pk)%p;
}
ll solve(ll n,ll m,ll x)
{ll ans=0,p=x;isprime();for(ll i=0;primer[i]*primer[i]<=x;++i){if(x%primer[i]==0){ll pk=1;while(x%primer[i]==0)x/=primer[i],pk*=primer[i];ans=(ans+ex_lucas(n,m,p,primer[i],pk))%p;}}if(x>1)ans=(ans+ex_lucas(n,m,p,x,x))%p;return  ans;
}
int main()
{ll n,m,p;cin>>n>>m>>p;ll ans=solve(n,m,p);printf("%lld\n",ans);return 0;
}

gym/100633/J Ceizenpok’s formula相关推荐

2015 ICL, Finals, Div. 1 Ceizenpok’s formula（组合数取模，扩展lucas定理）
J. Ceizenpok's formula time limit per test 2 seconds memory limit per test 256 megabytes input stand ...
[CodeForces gym 101630 J] 过路费（最短路）
problem 给定一张图 nnn 个点 mmm 条边,并给定阈值 kkk,以及起终点 s,ts,ts,t. 然后每条边经过都需要支付 www 的花费,形如 (u,v,w)(u,v,w)(u,v,w) ...
2019 CCPC秦皇岛 J 题 MUV LUV EXTRA【KMP 求最小循环节】
题意: 原题意是给出一个有理数的前一部分,求这个有理数的无限循环部分是什么.有一个值来评估猜的准确度.转换一下就成了下面的题意: 给出一个字符串s,有某一个子串,设p为该子串在s的某个后缀中的匹配长度 ...
python在化学方面的应用-用Python解析化学公式
问题是:给定类似C(Mg2(OH)4)2的公式,返回包含元素及其计数的哈希表.元素名称始终以大写字母开头,并且可以以小写字母开头. 我认为我首先要解决最简单的情况:没有括号. def bracket_ ...
周志华《机器学习》课后习题（第三章）：线性模型
作者 | 我是韩小琦链接 | https://zhuanlan.zhihu.com/p/43270830 3.1 试分析在什么情况下,在以下式子中不比考虑偏置项b. 答: 在样本中有某一个属性 ...
北邮OJ 1022. 16校赛-Saber's Board
时间限制 5000 ms 内存限制 65536 KB 题目描述 In a parallel universe, Saber has won the champion of every kind of ...
LeetCode 726. 原子的数量
给定一个化学式formula(作为字符串),返回每种原子的数量. 原子总是以一个大写字母开始,接着跟随0个或任意个小写字母,表示原子的名字. 如果数量大于 1,原子后会跟着数字表示原子的数量.如果数量 ...
DP! | 不要怂!
跟一个博客刷: http://blog.csdn.net/cc_again/article/details/25866971 一.简单基础dp 1.递推 HDU 2084 1 #include < ...
Qt4_Shreadsheet电子表应用
简介从Qt4--快速设计框开始到Qt4--子类化QTableWidgetItem,讲解了该应用的中每个函数的构建以及应用.详情前往这些章节里面阅读. main.cpp #include <QA ...

gym/100633/J Ceizenpok’s formula

gym/100633/J Ceizenpok’s formula相关推荐

最新文章

热门文章