![](https://img-blog.csdnimg.cn/20201014180756780.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数论
TheWolfWhistlingSong
蜜の夜明け
展开
-
Bzoj4403序列统计:Lucas初探,组合数学
Lucas初探,组合数学原创 2016-03-18 07:27:38 · 984 阅读 · 0 评论 -
Bzoj1004:[HNOI2008]Cards:置换群+dp
题目链接:1004:[HNOI2008]Cards 首先看到题目就知道是个有关群论的题 Burnside:真正染色方案数=(Σ每种置换下不变的染色方案数)/置换总数; 发现这里卡片的数量有限制所以Pólya并不能用了 于是要用dp求出不变的方案数 发现两个在某置换下可以互相到达的序列如果本质相同那么必须涂上相同颜色 可以处理处循环节长度,然后裸dp了 注意不洗牌也是一种方案 最后求原创 2016-04-21 10:38:37 · 460 阅读 · 0 评论 -
Bzoj2142礼物:组合数取模
题目链接:2142:礼物 组合数学推公式是很简单的,关键是取模 注意p并不是质数,而且p可能很大,所以lucas在这道题上并没有什么卵用 存在这样一个事实:设P分解质因数后有一项为pi^ai,那么C(x,y)%p%(pi^ai)=x%(pi^ai),设为xi 而题目中已经给出pi^ai 现在问题转化为怎么求xi 考虑把n!拆开得到1*2*3*4*5*6*...*n,显然和pi^ai不互原创 2016-04-20 17:31:44 · 737 阅读 · 0 评论 -
Bzoj2820:YY的GCD:莫比乌斯反演
题目链接:YY的GCD 公式和上一题一样地推 最后得到答案为 令T=pd',则 设 枚举素数贡献答案即可 #include #include #include #include #include #include #include #include #include #include #define LL long long using namespace std; const i原创 2016-04-19 10:34:30 · 488 阅读 · 0 评论 -
Bzoj2301: [HAOI2011]Problem b:莫比乌斯反演
题目链接:[HAOI2011]Problem b 妈妈我终于会反演啦【雾 其实我就是来码个公式复习一下QwQ 前置技能: 那么按照题意,我们令 首先明确我们为什么要反演,因为我们要求一个不好直接算的东西,而他和一个容易计算的东西有关系,所以我们要同过那个容易计算的东西来反向计算出我们真正要求的东西 于是我们令,其中?是整数,整个式子表示k的倍数 根据前置技能得 这时候发现d/k求起原创 2016-04-19 09:10:28 · 456 阅读 · 0 评论 -
Bzoj2242:[SDOI2011]计算器:BSGS
题目链接:2242:[SDOI2011]计算器 第一问快速幂,第二问扩展欧几里得,第三问BSGS 注意BSGS时要判断y为p的倍数时如果z不为0是无解的,wa了好久 为什么我跑得这样慢QAQ…… #include #include #include #include #include #include #define LL long long using namespace std; co原创 2016-04-12 14:09:45 · 381 阅读 · 0 评论 -
Bzoj3834:[Poi2014]Solar Panels:数论,分块
题目链接:3834: [Poi2014]Solar Panels 我们可以枚举G=gcd(i,j),并判断是否可行,但是这样一定会T 我们可以对G进行分块,使在每个块内对于a,b,c,d,中的任何一个它除以G的值相等 这样对于每个块我们只考虑最后一个G就可以了 #include #include #include #include #include using namespace std;原创 2016-03-22 20:57:45 · 653 阅读 · 0 评论 -
Bzoj2186:[Sdoi2008]沙拉公主的困惑:欧拉函数+乘法逆元
题目链接:2186:[Sdoi2008]沙拉公主的困惑 首先有 因为m可知有(n!-m!)/m!个 所以答案为 因为,其中pi是m!的质因子,所以答案是,pi是m!的质因子 其中m!的质因子就是m的质因子 还要预处理出pi的逆元,可以递推,公式ine[i]=(mod-mod/i*ine[mod%i]%mod); 离线处理即可,神犇说卡常数然而我并没有卡QwQ #include #i原创 2016-03-30 16:48:26 · 507 阅读 · 0 评论 -
Bzoj2190:[SDOI2008]仪仗队:欧拉函数
题目链接:2190:[SDOI2008]仪仗队 在网格图中,两个不同点之间如果没有其他点当且仅当这两个点的横纵坐标差的绝对值互质 那么这道题我们每加一层点就加一个欧拉函数的值即可,水得不行QwQ #include #include #define LL long long using namespace std; const int maxn=40010; int n,p[maxn];原创 2016-03-30 15:08:06 · 614 阅读 · 0 评论 -
Bzoj3782:上学路线:组合数学+Lucas+CRT+DP
题目链接:上学路线 设dp[i]为到了第i个坏点且是第一次到达坏点的方案数 那么dp[i]=C(x[i]+y[i],x[i])-∑C(x[i]+y[i]-x[j]-y[j],x[i]-x[j])dp[j](x[j] 对于mod=100003Lucas直接上 对于mod=1019663265分解质因数得到3*5*6973*10007然后每个因数分别Lucas后CRT合并 #include原创 2016-07-13 15:37:09 · 927 阅读 · 0 评论