2142: 礼物

Time Limit: 10 Sec  Memory Limit: 259 MB
Submit: 2286  Solved: 1009
[Submit][Status][Discuss]

Description

一年一度的圣诞节快要来到了。每年的圣诞节小E都会收到许多礼物,当然他也会送出许多礼物。不同的人物在小E
心目中的重要性不同,在小E心中分量越重的人,收到的礼物会越多。小E从商店中购买了n件礼物,打算送给m个人
,其中送给第i个人礼物数量为wi。请你帮忙计算出送礼物的方案数(两个方案被认为是不同的,当且仅当存在某
个人在这两种方案中收到的礼物不同)。由于方案数可能会很大,你只需要输出模P后的结果。

Input

输入的第一行包含一个正整数P,表示模;
第二行包含两个整整数n和m,分别表示小E从商店购买的礼物数和接受礼物的人数;
以下m行每行仅包含一个正整数wi,表示小E要送给第i个人的礼物数量。

Output

若不存在可行方案,则输出“Impossible”,否则输出一个整数,表示模P后的方案数。

Sample Input

100
4 2
1
2

Sample Output

12
【样例说明】
下面是对样例1的说明。
以“/”分割,“/”前后分别表示送给第一个人和第二个人的礼物编号。12种方案详情如下:
1/23 1/24 1/34
2/13 2/14 2/34
3/12 3/14 3/24
4/12 4/13 4/23
【数据规模和约定】
设P=p1^c1 * p2^c2 * p3^c3 * … *pt ^ ct,pi为质数。
对于100%的数据,1≤n≤109,1≤m≤5,1≤pi^ci≤10^5。

HINT

 

Source

 

[Submit][Status][Discuss]

数论大集合。

https://blog.csdn.net/clove_unique/article/details/54571216

答案显然$C_n^{n - w[1]}C_{n - w[1]}^{w[2]}C_{n - w[1] - w[2]}^{w[3]}......$,模数不互质成为难点。

扩展Lucas:将模数质因数分解,再用CRT合并。

问题就只剩下求$n!\%p_i^{k_i}$了。

首先由于不互质无法找到逆元,把$n!$中的$p_i$全部取出来:$ans=\lfloor\frac{n}{p}\rfloor+\lfloor\frac{n}{p^2} \rfloor+\lfloor\frac{n}{p^3}\rfloor+...$

然后考虑剩下部分怎么做,用$F(n,p_i,p_i^{k_i})$表示答案,$f(n,p_i,p_i^{k_i})$表示$\prod_{j=1,j\perp p_i}^{n}j(mod\ p_i^{k_i})$,则:

$F(n,p_i,p_i^{k_i})=\lfloor\frac{n}{p_i}\rfloor!\times p_i^{\lfloor\frac{n}{p_i}\rfloor}\times f(p_i^{k_i},p_i,p_i^{k_i})^{\lfloor\frac{n}{p_i^{k_i}}\rfloor}\times f(n\%p_i^{k_i})\%p_i^{k_i}$

递归处理即可,要用到exgcd。

复杂度好像是,$O(\sqrt{P}+m\log_{2}^{2}p_i^{k_i}\log{p_i})$不过肯定跑不满。

 #include<cstdio>
#include<algorithm>
#define rep(i,l,r) for (int i=l; i<=r; i++)
typedef long long ll;
using namespace std; const int N=;
ll mod,n,m,w[],ans,x,y,Mod[N],st[N],r[N],num; ll ksm(ll a,ll b,ll p){
ll res;
for (res=; b; a=a*a%p,b>>=)
if (b&) res=res*a%p;
return res;
} ll exgcd(ll a,ll b,ll &x,ll &y){
ll d=a;
if (b) d=exgcd(b,a%b,y,x),y-=a/b*x; else x=,y=;
return d;
} ll inv(ll t,ll p){ ll x,y; exgcd(t,p,x,y); return (x+p)%p; } ll F(ll n,ll pi,ll pk){
if (!n) return ;
ll ans=;
rep(i,,pk) if (i%pi) ans=ans*i%pk;
ans=ksm(ans,n/pk,pk);
rep(i,,n%pk) if (i%pi) ans=ans*i%pk;
return ans*F(n/pi,pi,pk)%pk;
} ll exlucas(ll n,ll m,ll pi,ll pk){
if (m>n) return ;
ll a=F(n,pi,pk),b=F(m,pi,pk),c=F(n-m,pi,pk);
ll k=;
for (ll i=n; i; i/=pi) k+=i/pi;
for (ll i=m; i; i/=pi) k-=i/pi;
for (ll i=n-m; i; i/=pi) k-=i/pi;
return a*inv(b,pk)%pk*inv(c,pk)%pk*ksm(pi,k,pk)%pk;
} ll CRT(ll n,ll r[],ll m[]){
ll M=,res=,w;
rep(i,,n) M*=m[i];
rep(i,,n) w=M/m[i],res=(res+w*inv(w,m[i])*r[i])%M;
return (res+M)%M;
} ll par(ll n,ll m[],ll st[]){
ll num=;
for (ll i=; i*i<=n; i++) if (n%i==){
ll pk=;
while (n%i==) pk*=i,n/=i;
m[++num]=pk; st[num]=i;
}
if (n>) m[++num]=n,st[num]=n;
return num;
} ll excomb(ll n,ll m){
rep(i,,num) r[i]=exlucas(n,m,st[i],Mod[i]);
return CRT(num,r,Mod);
} int main(){
freopen("bzoj2142.in","r",stdin);
freopen("bzoj2142.out","w",stdout);
scanf("%lld\n",&mod); scanf("%lld%lld",&n,&m);
ll sum=; rep(i,,m) scanf("%lld",&w[i]),sum+=w[i];
if (n<sum){ puts("Impossible"); return ; }
num=par(mod,Mod,st); ans=;
rep(i,,m) n-=w[i-],ans=ans*excomb(n,w[i])%mod;
printf("%lld\n",ans);
return ;
}

[BZOJ2142]礼物(扩展Lucas)的更多相关文章

  1. [bzoj2142]礼物(扩展lucas定理+中国剩余定理)

    题意:n件礼物,送给m个人,每人的礼物数确定,求方案数. 解题关键:由于模数不是质数,所以由唯一分解定理, $\bmod  = p_1^{{k_1}}p_2^{{k_2}}......p_s^{{k_ ...

  2. BZOJ2142 礼物 扩展lucas 快速幂 数论

    原文链接http://www.cnblogs.com/zhouzhendong/p/8110015.html 题目传送门 - BZOJ2142 题意概括 小E购买了n件礼物,送给m个人,送给第i个人礼 ...

  3. BZOJ.2142.礼物(扩展Lucas)

    题目链接 答案就是C(n,m1) * C(n-m1,m2) * C(n-m1-m2,m3)...(mod p) 使用扩展Lucas求解. 一个很简单的优化就是把pi,pi^ki次方存下来,因为每次分解 ...

  4. BZOJ - 2142 礼物 (扩展Lucas定理)

    扩展Lucas定理模板题(貌似这玩意也只能出模板题了吧~~本菜鸡见识鄙薄,有待指正) 原理: https://blog.csdn.net/hqddm1253679098/article/details ...

  5. bzoj 2142 礼物——扩展lucas模板

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2142 没给P的范围,但说 pi ^ ci<=1e5,一看就是扩展lucas. 学习材料 ...

  6. BZOJ2142礼物——扩展卢卡斯

    题目描述 一年一度的圣诞节快要来到了.每年的圣诞节小E都会收到许多礼物,当然他也会送出许多礼物.不同的人物在小E 心目中的重要性不同,在小E心中分量越重的人,收到的礼物会越多.小E从商店中购买了n件礼 ...

  7. bzoj2142 礼物——扩展卢卡斯定理

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2142 前几天学了扩展卢卡斯定理,今天来磕模板! 这道题式子挺好推的(连我都自己推出来了) , ...

  8. BZOJ2142 礼物 【扩展Lucas】

    题目 一年一度的圣诞节快要来到了.每年的圣诞节小E都会收到许多礼物,当然他也会送出许多礼物.不同的人物在小E 心目中的重要性不同,在小E心中分量越重的人,收到的礼物会越多.小E从商店中购买了n件礼物, ...

  9. BZOJ_2142_礼物_扩展lucas+组合数取模+CRT

    BZOJ_2142_礼物_扩展lucas+组合数取模 Description 一年一度的圣诞节快要来到了.每年的圣诞节小E都会收到许多礼物,当然他也会送出许多礼物.不同的人物在小E 心目中的重要性不同 ...

随机推荐

  1. SQL 与关系代数

    Table of Contents 前言 关系与表 关系代数的基本运算 投影 选择 并运算 集合差运算 笛卡尔积 更名运算 关系代数的附加运算 集合交运算 连接运算 自然连接 内连接 外连接 结语 前 ...

  2. Ubuntu15.04 python升级到python-3.6.x

    简略记录步骤,容后补充:   sudo add-apt-repository ppa:jonathonf/python-3.6   sudo apt-get update sudo apt-get i ...

  3. corosync.conf

    ##totem定义集群内各节点间是如何通信的,totem本是一种协议,专用于corosync专用于各节点间的协议,协议是有版本的 totem { ##版本号 version: ##安全认证on|off ...

  4. Vue打包app

    前言 公司之前用的app就是一个套壳挂个链接就能用的app,后来需要添加微信分享方便传播,没办法只好做成混合式的app了, 因为之前做.net用vs可以创建cordova项目也试着玩过,就决定用cor ...

  5. Java 打印* 三角形

    package anli1; public class sanjiaoxing { public static void main(String[] agrs){ System.out.println ...

  6. thinkphp3.2 学习笔记 基础篇

    环境要求:PHP5.3以上版本注意:PHP5.3DEV和php6不支持 目录结构 www WEB部署目录(或者子目录)├─index.php 入口文件├─README.md README文件├─App ...

  7. [AT2699]Flip and Rectangles

    题目大意:有一个$n\times m$的$01$矩阵,可以把任意行或列反转,问最大的全为一的子矩阵的面积 题解:有一个结论:若一个子矩形$S$中的任意一个$2\times 2$的子矩形都含有偶数个$1 ...

  8. [洛谷P3805]【模板】manacher算法

    题目大意:给你一个字符串,求出它的最长回文字段 题解:$manacher$算法 卡点:$p$数组未开两倍空间 C++ Code: #include <cstdio> #include &l ...

  9. python代理池的实现

    https://github.com/wangqifan/ProxyPool http://python.jobbole.com/86994/

  10. oracle的隐式游标

    游标的概念:     游标是SQL的一个内存工作区,由系统或用户以变量的形式定义.游标的作用就是用于临时存储从数据库中提取的数据块.在某些情况下,需要把数据从存放在磁盘的表中调到计算机内存中进行处理, ...