洛谷$P$3301 $[SDOI2013]$方程 $exLucas$+容斥
正解:$exLucas$+容斥
解题报告:
在做了一定的容斥的题之后再看到这种题自然而然就应该想到容斥,,,?
没错这题确实就是容斥,和这题有点儿像
注意下的是这里的大于和小于条件处理方式不同昂$QwQ$
对于大于等于,直接在一开始就先给它那么多就好,就先提前把$m-=\sum_{i=n_{1}+1}^{n_{1}+n_{2}} A_i$,这样就只剩小于等于的条件了
小于等于,一看最多就8个,显然就成了经典容斥套路题了鸭,
于是就枚举哪个爆了,然后可重排列搞下,容斥下,就做完了
放下推出来的式子趴,,,$\sum_{s} [\ |s|\ \&\ 1\ ]\cdot \binom{n-\sum_{i\in S}(x_{i}+1)+m-1}{m}$,大概是这样儿的,如果有问题我打完代码之后会$upd$的$QwQ$
但是还有一个要注意的就,这个$mod$不一定是质数昂,,,所以考虑用个$exLucas$就欧克了鸭
$over$
昂还有一个卡时间的小$tip$,,,就是在$exLucas$里最开始特判下,当$n<0$,$m<0$,$n<m$的时候就可以直接$return\ 0$了,这样就能成功从2571$ms$变成150$ms$,,,
#include<bits/stdc++.h>
using namespace std;
#define il inline
#define gc getchar()
#define int long long
#define ri register int
#define rb register bool
#define rc register char
#define rp(i,x,y) for(ri i=x;i<=y;++i) const int N=1e6+;
int fac_cnt,A[N],B[N],mod,wei[N],tmpp; il int read()
{
rc ch=gc;ri x=;rb y=;
while(ch!='-' && (ch>'' || ch<''))ch=gc;
if(ch=='-')ch=gc,y=;
while(ch>='' && ch<='')x=(x<<)+(x<<)+(ch^''),ch=gc;
return y?x:-x;
}
il void ex_gcd(ri a,ri b,ri &x,ri &y){if(!b){x=;y=;return;}ex_gcd(b,a%b,y,x);y-=(a/b)*x;}
il int power(ri x,ri y){/*if(!x)return 0;*/ri ret=;while(y){if(y&)ret=1ll*ret*x%mod;x=1ll*x*x%mod;y>>=;}return ret;}
struct nod
{
int p,pk,jc[N];
il void init(){jc[]=;rp(i,,pk)jc[i]=1ll*jc[i-]*(i%p?i:)%pk;}
il int power(ri x,ri y){ri ret=;while(y){if(y&)ret=1ll*ret*x%pk;x=1ll*x*x%pk;y>>=;}return ret;}
il int inv(ri d){ri x,y;ex_gcd(d,pk,x,y);return (x%pk+pk)%pk;}
il int multi(ri x){ri ret=;while(x){ret=ret*power(jc[pk],x/pk)%pk*jc[x%pk]%pk;x/=p;}return ret;}
il void getp(ri x,ri op,ri &k){while(x)k+=op*(x/p),x/=p;}
il int C(ri n,ri m)
{
ri a=multi(n),b=multi(m),c=multi(n-m),k=;getp(n,,k);getp(m,-,k);getp(n-m,-,k);
return 1ll*a*inv(b)%pk*inv(c)%pk*power(p,k)%pk;
}
nod(){pk=;}
}fac[];
il int crt()
{
ri p=A[],r=B[];
rp(i,,fac_cnt)
{ri x,y,np=A[i]*p;ex_gcd(p,A[i],x,y);x=(1ll*x*(B[i]-r)%np+np)%np;r=(1ll*x*p%np+r+np)%np;p=np;}
return r;
}
il int ex_lucas(ri n,ri m)
{
if(n< || m< || n<m)return ;
rp(i,,fac_cnt)A[i]=fac[i].pk,B[i]=fac[i].C(n,m);
return crt();
} main()
{
int T=read();tmpp=mod=read();
for(ri i=;i*i<=mod;++i){if(!(mod%i))fac[++fac_cnt].p=i;while(!(mod%i))fac[fac_cnt].pk*=i,mod/=i;}
if(mod>)fac[++fac_cnt].p=mod,fac[fac_cnt].pk=mod;mod=tmpp;
rp(i,,fac_cnt)fac[i].init();
while(T--)
{
ri n=read(),n1=read(),n2=read(),m=read()-n,S=<<n1,as=;
rp(i,,n1+n2)wei[i]=read()-;rp(i,n1+,n1+n2)m-=wei[i];
rp(i,,S-)
{
ri opt=,tmp_m=m;
rp(j,,n1)if(i&(<<(j-)))opt=mod-opt,tmp_m-=wei[j]+;
as=(as+1ll*opt*ex_lucas(tmp_m+n-,n-)%mod)%mod;
}
printf("%lld\n",as);
}
return ;
}
洛谷$P$3301 $[SDOI2013]$方程 $exLucas$+容斥的更多相关文章
- 洛谷P3172 [CQOI2015]选数(容斥)
传送门 首先,进行如下处理 如果$L$是$K$的倍数,那么让它变成$\frac{L}{K}$,否则变成$\frac{L}{K}+1$ 把$H$变成$\frac{H}{K}$ 那么,现在的问题就变成了在 ...
- 洛谷P1447 [NOI2010]能量采集(容斥)
传送门 很明显题目要求的东西可以写成$\sum_{i=1}^{n}\sum_{j=1}^m gcd(i,j)*2-1$(一点都不明显) 如果直接枚举肯定爆炸 那么我们设$f[i]$表示存在公因数$i$ ...
- 洛谷P1450 [HAOI2008]硬币购物 背包+容斥
无限背包+容斥? 观察数据范围,可重背包无法通过,假设没有数量限制,利用用无限背包 进行预处理,因为实际硬币数有限,考虑减掉多加的部分 如何减?利用容斥原理,减掉不符合第一枚硬币数的,第二枚,依次类推 ...
- bzoj 2005 & 洛谷 P1447 [ Noi 2010 ] 能量采集 —— 容斥 / 莫比乌斯反演
题目:bzoj 2005 https://www.lydsy.com/JudgeOnline/problem.php?id=2005 洛谷 P1447 https://www.luogu.org/ ...
- 洛谷P2522 [HAOI2011]Problem b (莫比乌斯反演+容斥)
题意:求$\sum_{i=a}^{b}\sum_{j=c}^{d}[gcd(i,j)==k]$(1<=a,b,c,d,k<=50000). 是洛谷P3455 [POI2007]ZAP-Qu ...
- 模板—点分治A(容斥)(洛谷P2634 [国家集训队]聪聪可可)
洛谷P2634 [国家集训队]聪聪可可 静态点分治 一开始还以为要把分治树建出来……• 树的结构不发生改变,点权边权都不变,那么我们利用刚刚的思路,有两种具体的分治方法.• A:朴素做法,直接找重心, ...
- 洛谷 P6295 - 有标号 DAG 计数(生成函数+容斥+NTT)
洛谷题面传送门 看到图计数的题就条件反射地认为是不可做题并点开了题解--实际上这题以我现在的水平还是有可能能独立解决的( 首先连通这个条件有点棘手,我们尝试把它去掉.考虑这题的套路,我们设 \(f_n ...
- 洛谷 P7360 -「JZOI-1」红包(Min-Max 容斥+推式子)
洛谷题面传送门 hot tea. 首先注意到这个 \(\text{lcm}\) 特别棘手,并且这里的 \(k\) 大得离谱,我们也没办法直接枚举每个质因子的贡献来计算答案.不过考虑到如果我们把这里的 ...
- 洛谷P5206 [WC2019]数树 [容斥,DP,生成函数,NTT]
传送门 Orz神仙题,让我长了许多见识. 长式子警告 思路 y=1 由于y=1时会导致后面一些式子未定义,先抓出来. printf("%lld",opt==0?1:(opt==1? ...
随机推荐
- Hbase数据模型物理视图
- oracle函数 RTRIM(c1,[,c2])
[功能]删除右边出现的字符串 [参数]C1 字符串 c2 追加字符串,默认为空格 [返回]字符型 [示例] SQL> select RTRIM('gao qian jingXXXX','X') ...
- 开发板ping通虚拟机与主机
刚因为虚拟机与主机没法互相ping通的事情,奋战到将近凌晨一点.现在把这个过程总结一下,以方便后加入该行业的广大IT精英. VMWare提供了三种工作模式:bridged(桥接模式).NAT(网络地址 ...
- SDUT-2054_数据结构实验之链表九:双向链表
数据结构实验之链表九:双向链表 Time Limit: 1000 ms Memory Limit: 65536 KiB Problem Description 学会了单向链表,我们又多了一种解决问题的 ...
- Libev源码分析02:Libev中的IO监视器
一:代码流程 在Libev中,启动一个IO监视器,等待该监视器上的事件触发,然后调用该监视器的回调函数.整个的流程是这样的: 首先调用ev_default_loop初始化struct ev_loop ...
- 数据库设计mysql字段不默认为NULL原因搜集
索引不会包含有NULL值的列 只要列中包含有NULL值都将不会被包含在索引中,复合索引中只要有一列含有NULL值,那么这一列对于此复合索引就是无效的.所以我们在数据库设计时不要让字段的默认值为NULL ...
- H3C路由器SSH服务配置命令
- avalon2 第一个demo
<!DOCTYPE html> <html> <head> <title>TODO supply a title</title> <m ...
- 使用jackson转json解决双向关联循环调用
ITOO V1.0的开发算是告一段落了,现在是整理总结交接环节,在这个项目中常见的问题也该好好整理一下和大家分享了,这次主要介绍转json循环调用的问题. 一.问题背景 相信只要使用ORM映射实体关联 ...
- C# 标准性能测试
经常我写一个类,作为一个工具类,小伙伴会问我这个类的性能,这时我就需要一个标准的工具进行测试. 本文告诉大家如何使用 benchmarkdotnet 做测试. 现在在 github 提交代码,如果有小 ...