数学
Adolphrocs
这个作者很懒,什么都没留下…
展开
-
矩阵翻硬币 数学 蓝桥杯
题目描述: 小明先把硬币摆成了一个 n 行 m列的矩阵。随后,小明对每一个硬币分别进行一次 Q 操作。对第x行第y列的硬币进行 Q 操作的定义:将所有第 i*x 行,第 j*y列的硬币进行翻转。其中i和j为任意使操作可行的正整数,行号和列号都是从1开始。 当小明对所有硬币都进行了一次 Q操作后,他发现了一个奇迹——所有硬币均为正面朝上。 想知道最开始有多少枚硬币是反面朝上的。只需要...原创 2018-03-27 22:17:41 · 289 阅读 · 0 评论 -
欧拉判别法 欧拉准则
若是奇质数且不能整除,则:是模的二次剩余当且仅当:是模的二次非剩余当且仅当:以勒让德符号表示,即为:原创 2019-08-11 16:23:22 · 4272 阅读 · 0 评论 -
莫比乌斯反演的两种形式及证明
莫比乌斯反演的两种形式及证明第一种形式若则证明右边等于===========已知==============所以当且仅当n/k =1,即n=k时,∑μ(d)=1,其余时候等于0。得证==============第二种形式若则证明右边等于===========已知==============得===...原创 2019-08-29 16:41:40 · 284 阅读 · 0 评论 -
BZOJ 2301 莫比乌斯反演模板题 Problem b
#include <bits/stdc++.h>#define FOR(i,s,t) for(int i=(s);i<=(t);i++)#define ROF(i,s,t) for(int i=(s);i>=(t);i--)#define pb push_back#define mp make_pair#define eb emplace_back#...原创 2019-08-29 17:11:28 · 210 阅读 · 0 评论 -
BZOJ 2820:YY的GCD 莫比乌斯反演
BZOJ2301的强化版上面这题是求gcd(i,j)==1的个数这题是求gcd(i,j) ==p的个数 p是素数我就在上面那个式子中额外枚举个素数p#include <bits/stdc++.h>#define FOR(i,s,t) for(int i=(s);i<=(t);i++)#define ROF(i,s,t) for(int i=(s);i>...原创 2019-08-30 14:35:06 · 197 阅读 · 0 评论 -
幂塔问题 拓展欧拉定理 连续b个a的幂次模c 2019南京网络赛B
拓展欧拉定理问题描述(连续b个a的幂次)#include <bits/stdc++.h>#define FOR(i,s,t) for(int i=(s);i<=(t);i++)#define ROF(i,s,t) for(int i=(s);i>=(t);i--)#define pb push_back#define mp make_pair#de...原创 2019-09-02 12:36:48 · 309 阅读 · 0 评论 -
2019上海网络赛 C [FFT 快速傅里叶变换]
题意:给出三组数组,从每个数组里挑出一个数组成一个三元组(a,b,c),使得这三个数可以组成三角形题解:由于题目n>1000的组数不超过20, FFT的大小与数量无关与最大值有关题源:这题是 HDU 4609 的改编https://blog.csdn.net/Adolphrocs/article/details/101555137/* *FFT */#include &l...原创 2019-09-27 19:27:00 · 251 阅读 · 0 评论 -
HDU 4609 3-idiots [FFT 快速傅里叶变换]
3-idiots Time Limit: 10000/5000 MS (Java/Others)Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 10046Accepted Submission(s): 3406 Problem Description ...原创 2019-09-27 19:33:13 · 276 阅读 · 0 评论 -
洛谷P4157 [SCOI2006]整数划分 高精度乘法[FFT 快速傅里叶变换]
题目描述从文件中读入一个正整数n(10≤n≤31000)。要求将n写成若干个正整数之和,并且使这些正整数的乘积最大。例如,n=13,则当n表示为4+3+3+3(或2+2+3+3+3)时,乘积=108为最大。输入格式只有一个正整数:n (10≤n≤31000)输出格式第1行输出一个整数,为最大乘积的位数。第2行输出最大乘积的前100位,如果不足100位,则按实际位数输出最...原创 2019-09-28 10:30:19 · 480 阅读 · 0 评论 -
指数型母函数 生成函数 [2019上海网络赛 E题 Counting Sequences II]
泰勒展开公式常见公式。指数型母函数例题Counting Sequences II题意:构造一盒序列a[1,n],满足1<=a[i] <=m,其中偶数出现的次数是偶数次#include <bits/stdc++.h>using namespace std;typedef long long ll;const int maxn...原创 2019-10-01 16:11:48 · 409 阅读 · 0 评论 -
母函数 生成函数 模板 HDU 1028 Ignatius and the Princess III
题意:给出一个N,问N的划分数题解:显然n的划分数等于(1+x+x^2+...)*(1+x^2+x^4+...)*(1+x^3+x^6+...)的结果中x^n次前的系数#include <bits/stdc++.h>using namespace std;typedef long long ll;const int maxn = 200;int n;int c1[...原创 2019-10-01 17:46:29 · 185 阅读 · 0 评论 -
母函数 生成函数 模板 HDU 1398 Square Coins
题意:给出n,求n只有平方数组成的划分题解:显然n的划分数等于(1+x+x^2+...)*(1+x^4+x^8+...)*(1+x^9+x^16+...)*(1+x^的结果中x^n次前的系数#include <bits/stdc++.h>using namespace std;typedef long long ll;const int maxn = 200;int...原创 2019-10-01 18:01:26 · 191 阅读 · 0 评论 -
二次剩余 数论 勒让德
在数论中,特别在同余理论里,一个整数对另一个整数的二次剩余(英语:Quadratic residue)指的平方除以得到的余数。当存在某个,式子成立时,称“是模的二次剩余”当对任意,不成立时,称“d是模p的二次非剩余”质数的二次剩余对于质数2,每个整数都是它的二次剩余。以下讨论是奇质数的情况:对于,而言,能满足“是模 的二次剩余”的共有个(剩余类),分别为:(0计算在内)...原创 2019-08-10 20:29:23 · 1620 阅读 · 0 评论 -
逆元 阶乘逆元 数论 HDU 5651 xiaoxin juju needs help
逆元用于用于对除法运算取模的运算表示的逆元求解逆元方法一:拓展欧几里得求逆元(求单个数的逆元)已知同时取模得可见是的逆元所以中的int ex_gcd(int a, int b, int &x, int &y){ if (b == 0){ x = 1; y = 0; return ...原创 2019-08-10 17:09:19 · 271 阅读 · 0 评论 -
Covering HDU6185
Bob's school has a big playground, boys and girls always play games here after school. To protect boys and girls from getting hurt when playing happily on the playground, rich boy Bob decided to cov...原创 2018-07-14 21:51:21 · 212 阅读 · 0 评论 -
CS Course HDU 6186
CS CourseTime Limit: 4000/2000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 1964 Accepted Submission(s): 812Problem DescriptionLittle A has come to college and ...原创 2018-07-14 22:30:01 · 181 阅读 · 0 评论 -
polya 定理+ 手敲模板
polya 定理处理的问题 :假设你是个卖项链的,你有m种颜色的珍珠,每串项链有n个珍珠,问一共能做多少种不同的项链(因为客户不是傻子,所以项链通过旋转、反转变相同的 算同一种项链)https://blog.csdn.net/xuzengqiang/article/details/7476671https://blog.csdn.net/liangzhaoyang1/article/d...原创 2018-07-21 18:39:51 · 336 阅读 · 0 评论 -
莫比乌斯反演
感觉网上写的都不太看得懂,自己重新写先看看莫比乌斯函数(不用看得懂)https://blog.sengxian.com/algorithms/mobius-inversion-formulavoid sieve() { fill(isPrime, isPrime + maxn, 1); mu[1] = 1, num = 0; for (int i = 2; i...原创 2018-10-11 15:10:41 · 1725 阅读 · 0 评论 -
欧拉函数 欧拉定理 欧拉降幂
洛谷 P5091φ(n)表示小于n的正整数与n互质数的个数第一种情况:如果n=1, φ(1) = 1第二种情况:如果n是质数,φ(n) = n -1;第三种情况:如果n是质数的某个次方=第四种情况:如果n两个互质数的乘积第五种情况:如果n任意大于1的整数 int phi(int n){ int res=n,a=n; ...原创 2019-08-05 16:40:52 · 329 阅读 · 0 评论 -
stein算法(快速gcd)
设x,y为非0奇数,有以下结论:#include<bits/stdc++.h>using namespace std;#define LL long longLL stein(LL a, LL b) { if(!a) return b; if(!b) return a; if(!(a | 1) && !(b | 1)) r...原创 2019-08-08 15:02:36 · 760 阅读 · 0 评论 -
威尔逊定理 数论
威尔逊定理n为素数时应用YAPTCHA YAPTCHA Time Limit: 10000/5000 MS (Java/Others)Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 1858Accepted Submission(s): 960 ...原创 2019-08-06 11:57:34 · 253 阅读 · 0 评论 -
中国剩余定理 孙子定理(互质与不互质) FZU1402(互质) POJ2891(不互质)
中国剩余定理先看个互质的解得 /* *中国剩余定理 互质 */#include <stdio.h>#include <cstring>#include <iostream>using namespace std;typedef long long LL;///n个mi互质const LL maxn = 20;LL ...原创 2019-08-06 17:40:35 · 425 阅读 · 0 评论 -
洛谷P5091 扩展欧拉函数
扩展欧拉函数/* *扩展欧拉函数+快速幂+(龟速乘) */#include <bits/stdc++.h>using namespace std;typedef unsigned long long ll;ll n, m, phi , t, b;ll Pow(ll x,ll p){ if (p == 0) return 1ll; ll tmp =...原创 2019-08-06 17:44:50 · 225 阅读 · 0 评论 -
第一类斯特林数 First Kind of Stirling Number
题目描述有n个高度分别为1到n的人,从前到后站成一列,站在前面的人会挡住后面比他矮的人。现从前面一共能看见m个人,请问有多少种可能的站法输入格式第一行两个整数n,m含义如图所述输出格式一行一个整数,为答案模1e9+7的值数据范围时间限制1s,空间限制:512M样例1Input3 2Output3样例1解释一共有132、231、213...原创 2019-08-07 10:46:15 · 623 阅读 · 0 评论 -
费马小定理 费马大定理 数论
a与m互质时费马小定理费马小定理降幂求逆元费马大定理不存在 n > 2 使得n = 2时a为奇数时 a 、、a为偶数时 a 、、当 a 为奇数时,则 a = 2 * k + 1 ,解得 k 的值,则 b =2 * k * ( k + 1 ),c =2 * k * ( k + 1 ) + 1;当 a 为偶数时,则 a...原创 2019-08-07 16:59:36 · 518 阅读 · 0 评论 -
母函数模板
(x^(v[0]*n1[0])+x^(v[0]*(n1[0]+1))+x^(v[0]*(n1[0]+2))+...+x^(v[0]*n2[0]))(x^(v[1]*n1[1])+x^(v[1]*(n1[1]+1))+x^(v[1]*(n1[1]+2))+...+x^(v[1]*n2[1]))...(x^(v[K]*n1[K])+x^(v[K]*(n1[K]+1))+x^(v[K]*(n1[K...原创 2019-10-02 23:47:27 · 160 阅读 · 0 评论