king9666的博客

根据威尔逊定理,可得如果p是素数,那么就能找到这样的m值可使得 (p-1)! = (p-1) + mp或者可以说成 (p-1)! = (p-1) mod p举个例子6! = 1 * 2 * 3 * 4 * 5 * 6原式 = 1 * (2 * 4) * (3 * 5) * 6原式 = (1 * 1 * 1* 6) mod 7原式 = 6 mod 7PS: 也可以用模乘验证费马小定理....

拉格朗日的作用就是给你 n+1 个点，让你构造出一个 n 次函数，使得这个函数的图像经过坐标轴上的 n+1 个点。重心拉格朗日插值法优化到O(n)O(n*n)#include<iostream>#include<cstdio>#define ll long longusing namespace std;const int M=2111;const ll...

拉格朗日插值公式指的是在节点上给出节点基函数，然后做基函数的线性组合，组合系数为节点函数值的一种插值多项式。线性插值也叫两点插值，已知函数y = f (x)在给定互异点x0, x1上的值为y0= f (x0)，y1=f (x1)线性插值就是构造一个一次多项式：P1(x) = ax + b，使它满足条件：P1 (x0) = y0， P1 (x1) = y1其几何解释就是一条直线，通过已知点A ...

https://www.cnblogs.com/zwfymqz/p/10063039.html

Number Gamedescription：Alice and Bob are playing a game on a line of N squares. The line is initially populated with one of each of the numbers from 1 to N. Alice and Bob take turns removing a singl...

我还不会证…后补

莫比乌斯反演用于解决排列组合问题的数论知识.应用例子: HDU1695后补

POJ 3904入门莫比乌斯反演(数论)

传送门后补题

莫比乌斯反演莫比乌斯反演是数论数学中很重要的内容，可以用于解决很多组合数学的问题。

传送门题意:给定n，求满足1/a+1/b=1/c，求a,b,c互质的三元组(a,b,c)数目。（a,b,c≤n）n≤10的12次方.后补

后补 …

????????????????????费马小定理欧拉定理阶乘逆元ll j[maxn];ll in[maxn];void pre_solve(int n,ll p){ jc[0]=1; for(int i=1;i<=n;i++) j[i]=(j[i-1]*i)%p; in[n]=qpow(j[n],p-2,p); // 快速幂 for(int i=n-1;i...

ll qmul(ll a,ll b,ll p){ ll ans=0,base=a%p; while(b) { if(b&1) ans=(ans+base)%p; base=(base+base)%p; b>>=1; } return ans;}

Lucas定理是用来求 c(n,m) mod p，p为素数的值。它是用来求组合数 C(n, m) % p 值的定理，这里的p是素数。所以，它是一个解决大组合数求模的算法。卢卡斯定理适用于模数 ???? 为质数的情况首先我们知道组合数的计算公式????????????=????!/????!(????−????)!所以我们可以先递推出阶乘再用快速幂或 ???????????????????? 计算逆元也可以直接递推阶乘逆元然后就可以计算 ????了ll qpow(l...

三角形数1 3 6 10 15 21n * ( n + 1 ) / 2长方形数2 6 12 20 30 42(n+1)* nGnomon数1 3 5 7 92n - 1正方形数前N个Gnomon数的求和公式就是正方形数n的平方

GCM最大公因数差不多,最大公度量line_segment gcm(line_segment a, line_segment b){if(a == b) {return a;}if(b < a){return gcm(a - b, b);}/*if(a < b)return gcm(a, b-a);*/}...

// 机器学习与深度学习 C语言实现// 基于归纳学习的股价预测程序#include <stdio.h>#include <stdlib.h>#define OK 1#define NG 0#define SETSIZE 100#define CNO 10#define GENMAX 100000#define SEED 32767int calc...

// 蚁群优化算法#include "stdio.h"#include "stdlib.h"#include "math.h"#define NOA 10#define ILIMIT 50#define Q 3#define RHO 0.8#define STEP 10#define EPSILON 0.15#define SEED 32768void up...

稍后补

１． FCFS 算法实现C++ 代码如下实现的方法： 模拟简单的数据结构代码功能：实现现有数据的ＣＰＵ调度算法　　 #include <iostream>#include <string>#include <cstdio>#include <algorithm>#include <set>#include <vec...

Ping pongTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 6476 Accepted Submission(s): 2398Problem DescriptionN(3<=N<=20000) ping po...

Otvio is a calm guy with a guilty pleasure: kicking PCs. Now he is participating in a programming competition where the computers are arranged as a tree, that is, the N PCs are connected by N - 1 cabl...

A - Nicoleta and the circle of kidsNicoleta is a happy kid, he goes every day to school and is always complaining about waking up early and having to complete his homework assignments. His favorite p...

Sum Of GcdTime Limit: 10000/5000 MS (Java/Others) Memory Limit: 65535/32768 K (Java/Others)Total Submission(s): 986 Accepted Submission(s): 465Problem DescriptionGiven you a sequence of numb...

D. Powerful arraytime limit per test5 secondsmemory limit per test256 megabytesinputstandard inputoutputstandard outputAn array of positive integers a1, a2, …, an is given. Let us consider its ar...

Given a sequence of n numbers a1, a2, …, an and a number of d-queries. A d-query is a pair (i, j) (1 ≤ i ≤ j ≤ n). For each d-query (i, j), you have to return the number of distinct elements in the su...

题目描述作为一个生活散漫的人，小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿。终于有一天，小Z再也无法忍受这恼人的找袜子过程，于是他决定听天由命……具体来说，小Z把这N只袜子从1到N编号，然后从编号L到R(L 尽管小Z并不在意两只袜子是不是完整的一双，甚至不在意两只袜子是否一左一右，他却很在意袜子的颜色，毕竟穿两只不同色的袜子会很尴尬。你的任务便是告诉小Z，他有多大的概率抽到两只...

Problem Description给定一个含有n个数字的数列，每个数字都有一个值a[i]（下标从1开始）。定义第i个数字和第j个数字间的距离dis(i,j)=abs(i-j)。接下来给出q个询问，每次询问一个区间[l,r]，要求求出一对数字(i,j)(l<=i<=j<=r),使得a[i]=a[j]并且dis(i,j)最大，由于这样的数对可能有多个，因此答案只要输出dis。...

GroupThere are n men ,every man has an ID(1…n).their ID is unique. Whose ID is i and i-1 are friends, Whose ID is i and i+1 are friends. These n men stand in line. Now we select an interval of men to...

稍后补上

奇异局势（必败态）博弈是不公平的游戏 因为只要双方足够聪明 从游戏开始就已经确定了结果在我们的博弈游戏中 想获取胜利 就要寻找必输状态要寻找必败状态 首先要知道什么情况下算输在游戏规则下 轮到你了却无法进行操作 就认定为输了（前提）然后 我们寻找必败态想一下 如果你面对必败态 那么你做任何操作 都将会把局面转化为非必败态（如果你能转化为必败态 你对手岂不是输了~~~~...

错排公式核心递推公式：D(n) = (n-1) [D(n-2) + D(n-1)]特殊地，D(1) = 0, D(2) = 1.问题: 十本不同的书放在书架上。现重新摆放，使每本书都不在原来放的位置。有几种摆法?这个问题推广一下，就是错排问题，是组合数学中的问题之一。考虑一个有n个元素的排列，若一个排列中所有的元素都不在自己原来的位置上，那么这样的排列就称为原排列的一个错排。 n个元素的...

费马小定理(Fermat’s little theorem)费马小定理是数论中的一个重要定理，在1636年提出。如果p是一个质数，而整数a不是p的倍数，则有a^（p-1）≡1（mod p）。费马小定理是初等数论四大定理（威尔逊定理，欧拉定理（数论中的欧拉定理），中国剩余定理(又称孙子定理）之一，在初等数论中有着非常广泛和重要的应用。实际上，它是欧拉定理的一个特殊情况（即 ，见于词条“欧拉函数...

不要62Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 66538 Accepted Submission(s): 26439Problem Description杭州人称那些傻乎乎粘嗒嗒的人为62（音：laoer）。杭州交...

埃氏筛（埃拉特斯特尼筛法）int visit[maxn];void Prime(){mem(visit,0); //初始化都是质数visit[0] = visit[1] = 1; //0 和 1不是质数for (int i = 2; i <= maxn; i++) {if (!visit[i]) { //如果i是素数，让i的所有倍数都不是质数...

https://linxi99.gitee.io/20190211/ACM几何基础篇/

https://www.bilibili.com/video/av27156563

BombTime Limit: 2000/1000 MS (Java/Others) Memory Limit: 131072/65536 K (Java/Others)Total Submission(s): 26039 Accepted Submission(s): 9876Problem DescriptionThe counter-terrorists found a ...

快速幂对大数时间复杂度的优化，具体操作是利用二进制操作11的二进制是1011，11 = 2³×1 + 2²×0 + 2¹×1 + 2º×1，因此，我们将a¹¹转化为算 a(20)*a(21)*a(23) ，看出来快的多了吧原来算11次，现在算三次& 运算还可以判断奇偶x&1== 0为偶，x&1== 1为奇int ksm ( int a, int b){ i...