题链:

http://www.lydsy.com/JudgeOnline/problem.php?id=3129

题解:

容斥,扩展Lucas,中国剩余定理

先看看不管限制,只需要每个位置都是正整数时的方案数的求法。
假设有 N 个未知数,加起来的和为 M。
转化一下问题变为:"小球分配"
有 M 个相同的小球,放在 N 个盒子里,且每个盒子至少有一个的方案数。
那么方案数为 ${C}_{M-1}^{N-1}$
怎么理解这个式子呢?"插隔板法"。
使 M个小球放在一排,考虑可以在相邻小球的空隙间(共 M-1个空隙)插入一个隔板,总共插入 N-1个隔板。
把相邻的两个隔板(把首尾也看作另外2个隔板)中间的区域看做一个个的盒子,区域内的小球即为该盒子里的小球。
则任意一种插隔板的方法都对应一种把小球放入盒子的方法。
所以方案数为
${C}_{M-1}^{N-1}$

对于题目给的两种限制,不要被吓到了。
其实大于等于(>=W)这一个限制很好处理,直接把 M-=(W-1),即事先给这些位置分配 W-1个小球。
因为接下来按照上面的 "小球分配" 方式,这些限制一定可以满足。
然后对于小于等于(<=W) 就需要用到容斥了。
定义 f[S] 表示至少 S 集合里面的这些"小于等于"限制都不能满足,那么方案数的求法就是:
先强制给这些盒子分配 Wi 个小球,
使得接下来按照上面的 "小球分配" 方式求出方案数,这些盒子一定会超过限制。
所以按照容斥的做法:
答案 ANS = 至少 0 个盒子超出限制的方案数
                - 至少 1 个盒子超出限制的方案数
                +至少 2 个盒子超出限制的方案数
                -...+...

但是求那个方案数 C(N,M) 就比较麻烦了,因为 N,M很大,且模数不保证为质数,
所以用到 扩展Lucas(+国剩) 来求。

代码:

#include<cstdio>
#include<cstring>
#include<iostream>
#define _ % P
#define filein(x) freopen(#x".in","r",stdin);
#define fileout(x) freopen(#x".out","w",stdout);
using namespace std;
int h[300],k[300],d[300],pi[10],ppi[10],fc[10][100500],pnt;
int T,N,N1,N2,M,ANS;
void init(){
ANS=0;
memset(h,0,sizeof(h));
memset(k,0,sizeof(k));
memset(d,0,sizeof(d));
}
void pre_prime(int P){
static bool np[100500];
static int prime[100500],cnt;
for(int i=2;i<=100000;i++){
if(!np[i]){
prime[++cnt]=i;
if(P%i==0){
pi[++pnt]=i; ppi[pnt]=1;
while(P%i==0) ppi[pnt]*=i,P/=i;
}
}
for(int j=1;j<=cnt&&i<=100000/prime[j];j++){
np[i*prime[j]]=1;
if(i%prime[j]==0) break;
}
}
for(int i=1;i<=pnt;i++){
fc[i][0]=1;
for(int j=1;j<=100000;j++)
fc[i][j]=(1ll*fc[i][j-1]*(j%pi[i]?j:1))%ppi[i];
}
}
int pow(int a,int b,int P){
int ret=1; a=(a)_;
while(b){
if(b&1) ret=(1ll*ret*a)_;
a=(1ll*a*a)_; b>>=1;
}
return ret;
}
void exEuclidean(int a,int b,int &g,long long &x,long long &y){
if(!b) g=a,x=1,y=0;
else exEuclidean(b,a%b,g,y,x),y-=(a/b)*x;
}
int inv(int a,int P){
int g; long long x,y;
exEuclidean(a,P,g,x,y);
x=((1ll*x)_+P)_;
return (int)x;
}
int fac(int n,int I,int Pi,int P){
if(!n) return 1; int ret=1,r=n%P;
//for(int i=2;i<=P;i++) if(i%Pi) ret=(1ll*ret*i)_;
ret=fc[I][P];
ret=pow(ret,n/P,P);
//for(int i=2;i<=r;i++) if(i%Pi) ret=(1ll*ret*i)_;
ret=(1ll*ret*fc[I][r])_;
return (1ll*ret*fac(n/Pi,I,Pi,P))_;
}
int C(int n,int m,int I,int Pi,int P){
int x,y,z,c=0;
x=fac(n,I,Pi,P); y=fac(m,I,Pi,P); z=fac(n-m,I,Pi,P);
for(int i=n;i;i/=Pi) c+=i/Pi;
for(int i=m;i;i/=Pi) c-=i/Pi;
for(int i=n-m;i;i/=Pi) c-=i/Pi;
return ((((1ll*x*inv(y,P))_)*inv(z,P))_*pow(Pi,c,P))_;
}
int exLucas(int n,int m,int P){
int ret=0,tmp; if(n<m) return 0;
for(int i=1;i<=pnt;i++){
tmp=C(n,m,i,pi[i],ppi[i]);
tmp=((1ll*tmp*(P/ppi[i]))_*inv(P/ppi[i],ppi[i]))_;
ret=(1ll*ret+tmp)_;
}
return ret;
}
int main()
{
int P;
scanf("%d%d",&T,&P);
pre_prime(P);
while(T--){
init();
scanf("%d%d%d%d",&N,&N1,&N2,&M);
for(int i=1;i<=N1;i++) scanf("%d",&d[1<<(i-1)]);
for(int i=1,x;i<=N2;i++) scanf("%d",&x),M-=(x-1);
for(int s=1,p;p=s&-s,s<(1<<N1);s++)
h[s]=h[s^p]+1,k[s]=k[s^p]+d[p];
for(int s=0,m,tmp;s<(1<<N1);s++){
m=M-k[s];
tmp=exLucas(m-1,N-1,P);
if(h[s]&1) tmp=(-1ll*tmp+P)_;
ANS=((1ll*ANS+tmp)_+P)_;
}
printf("%d\n",ANS);
}
return 0;
}

●BZOJ 3129 [Sdoi2013]方程的更多相关文章

  1. [BZOJ 3129] [Sdoi2013] 方程 【容斥+组合数取模+中国剩余定理】

    题目链接:BZOJ - 3129 题目分析 使用隔板法的思想,如果没有任何限制条件,那么方案数就是 C(m - 1, n - 1). 如果有一个限制条件是 xi >= Ai ,那么我们就可以将 ...

  2. BZOJ 3129 [SDOI2013]方程 (拓展Lucas)

    题目大意:给定一个方程$X_{1}+X_{2}+X_{3}+X_{4}+...+X_{n}=M$,$\forall X_{i}<=A_{i} (i<=n1)$ $\forall X_{i} ...

  3. BZOJ 3129 SDOI2013 方程

    如果没有限制,答案直接用隔板法C(m-1,n-1) 对于>=x的限制,我们直接在对应位置先放上x-1即可,即m=m-(x-1) 对于<=x的限制,由于限制很小我们可以利用容斥原理将它转化为 ...

  4. bzoj千题计划267:bzoj3129: [Sdoi2013]方程

    http://www.lydsy.com/JudgeOnline/problem.php?id=3129 如果没有Ai的限制,就是隔板法,C(m-1,n-1) >=Ai 的限制:m减去Ai &l ...

  5. bzoj3129[Sdoi2013]方程 exlucas+容斥原理

    3129: [Sdoi2013]方程 Time Limit: 30 Sec  Memory Limit: 256 MBSubmit: 582  Solved: 338[Submit][Status][ ...

  6. 【BZOJ3129】[SDOI2013]方程(容斥,拓展卢卡斯定理)

    [BZOJ3129][SDOI2013]方程(容斥,拓展卢卡斯定理) 题面 BZOJ 洛谷 题解 因为答案是正整数,所先给每个位置都放一个就行了,然后\(A\)都要减一. 大于的限制和没有的区别不大, ...

  7. BZOJ_3129_[Sdoi2013]方程_组合数学+容斥原理

    BZOJ_3129_[Sdoi2013]方程_组合数学+容斥原理 Description 给定方程     X1+X2+. +Xn=M 我们对第l..N1个变量进行一些限制: Xl < = A ...

  8. bzoj 3129

    非常好的一道数学题,考察了大量数论和组合数学的知识 在做本题之前强烈建议先完成下列两个背景知识: ①: bzoj 2142礼物 因为本题的一部分数据需要利用到拓展卢卡斯定理,而礼物是拓展卢卡斯定理的裸 ...

  9. [BZOJ 3198] [Sdoi2013] spring 【容斥 + Hash】

    题目链接:BZOJ - 3198 题目分析 题目要求求出有多少对泉有恰好 k 个值相等. 我们用容斥来做. 枚举 2^6 种状态,某一位是 1 表示这一位相同,那么假设 1 的个数为 x . 答案就是 ...

随机推荐

  1. XCode Build Settings中几种Search Paths

    Header search path:去查找头文件的路径,同在在你需要使用第三方库的时候,在这里设置你的头文件路径目录,如图 <code><span class="str& ...

  2. Tornado介绍及自定义组件

    Tornado 的性能是相当优异的,因为它试图解决一个被称之为"C10k"问题,就是处理大于或等于一万的并发.一万呀,这可是不小的量 条件:处理器为 AMD Opteron, 主频 ...

  3. Java如何调取创蓝253短信验证码

    基于创蓝253短信服务平台的Java调用短信接口API package com.bcloud.msg.http; import java.io.ByteArrayOutputStream; impor ...

  4. 测试驱动开发实践3————从testList开始

    [内容指引] 运行单元测试: 装配一条数据: 模拟更多数据测试列表: 测试无搜索列表: 测试标准查询: 测试高级查询. 一.运行单元测试 我们以文档分类(Category)这个领域类为例,示范如何通过 ...

  5. 学习phalcon框架按照官网手册搭建第一个项目注册功能

    中文手册官网:http://phalcon.ipanta.com/1.3/tutorial.html#bootstrap 官网提供http://www.tutorial.com项目源码github地址 ...

  6. 由路由器AP隔离引起的WEB服务不能访问的问题

    最近在自己的电脑上运行了一个Web服务,发现它只能被本机访问,无法在同一个局域网的其它设备访问. 于是在网上搜索了各种资料,通通都是在操作系统层面来分析解决问题的. 而我怎么尝试都不成功的情况下,差点 ...

  7. Python中字符串颜色

    格式:\033[显示方式;前景色;背景色m 说明: 前景色 背景色 颜色 --------------------------------------- 30 40 黑色 31 41 红色 32 42 ...

  8. Codeforces Round #436 (Div. 2) B. Polycarp and Letters

    http://codeforces.com/contest/864/problem/B 题意: 给出一个字符串,要求找到一个集合S,使得从S中选出的所有数,在这些数的位置上的字母全部为小写且是不同的字 ...

  9. codeforces round 425 div2

    A. Sasha and Sticks 水题,判断一下次数的奇和偶就可以的. B. Petya and Exam 赛上的时候没有写出来,orz,记录一下吧. 题意:给出一个模式串,可能会有?和*两种符 ...

  10. [转]Python爬虫框架--pyspider初体验

    标签: python爬虫pyspider 2015-09-05 10:57 9752人阅读 评论(0) 收藏 举报  分类: Python(8)  版权声明:本文为博主原创文章,未经博主允许不得转载. ...