![](https://img-blog.csdnimg.cn/20201014180756930.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
caioj数论
zsyz_lb2003
这个作者很懒,什么都没留下…
展开
-
Area[POJ1654]
题面描述传送门思路瞎枚举算面积就好了。(一点技术含量都没有)#include<cstdio>#include<algorithm>#include<cstdlib>#include<cmath>#include<cstring>#define ll long longusing namespace std;cons...原创 2019-05-28 13:12:31 · 94 阅读 · 0 评论 -
GCD4[caioj1283]\[Zoj3435]
欢迎大家访问我的老师的OJ———caioj.cn题面描述传送门思路算是仪仗队的进阶版吧。因为这题的左下角为(1,1,1)(1,1,1)(1,1,1),我们不妨将其变为(0,0,0)(0,0,0)(0,0,0),右上角为(a−1,b−1,c−1)(a-1,b-1,c-1)(a−1,b−1,c−1)方便讨论。首先我们可以观察到,三条棱上的(1,0,0),(0,1,0),(0,0,1)(1,...原创 2019-04-30 10:57:06 · 146 阅读 · 0 评论 -
Fibonacci[POJ3070]
题面描述传送门思路原创 2019-04-25 14:22:07 · 130 阅读 · 0 评论 -
矩阵乘法原理
欢迎大家访问我的老师的OJ———caioj.cn一个n∗mn*mn∗m的矩阵可看作一个n∗mn*mn∗m的二维数组。矩阵的加法和减法,仅满足于两个同样大小(同为n∗mn*mn∗m)的矩阵,就是把矩阵对应位置上的数相加减,即C=A+B↔∀i∈[1,n],∀j∈[1,m],Ci,j=Ai,j+Bi,jC=A+B \leftrightarrow\forall i \in[1,n],\forall...原创 2019-04-25 14:21:20 · 2287 阅读 · 0 评论 -
GCD3[caioj1282]\[BZOJ2301]莫反
题面描述caiojbzoj注意,caioj题面略有不同思路与Zap很像。设D(a,b,k)D(a,b,k)D(a,b,k)表示满足x≤a,y≤bx\le a,y\le bx≤a,y≤b且k∣gcd(x,y)k \mid \gcd(x,y)k∣gcd(x,y)的二元组有多少对。由于D(a,b,k)=⌊a/k⌋∗⌊b/k⌋D(a,b,k)=\left\lfloor a/k\right\...原创 2019-04-30 09:08:49 · 187 阅读 · 0 评论 -
GCD2[caioj1281]\[BZOJ2818]莫反
欢迎大家访问我的老师的OJ———caioj.cn题面描述caiojBZOJ思路这是一道不用莫函的莫反(伪)题啊。我们观察一下,1≤x,y≤n,gcd(x,y)=p,p是质数1\le x,y\le n,\gcd(x,y)=p,p是质数1≤x,y≤n,gcd(x,y)=p,p是质数可以转化为:1≤x,y≤n/p,gcd(x,y)=1,也就是x,y互质1\le x,y\le n/p...原创 2019-04-30 07:55:25 · 129 阅读 · 0 评论 -
Strange Way to Express Integers[POJ2891]
欢迎大家访问我的老师的OJ———caioj.cn题目描述传送门思路详细参考这道题输入有坑,且多组数据。#include<cstdio>#include<cstring>#include<algorithm>#include<cstdlib>#include<cmath>#define ll long longusi...原创 2019-04-25 10:43:14 · 121 阅读 · 0 评论 -
[计算几何]旋转卡壳
题面描述传送门思路凸包可以打SCY的。去看wyjjj的卡壳吧SCY的std有点问题的。对于卡壳这一部分while(mul(h[i+1],h[now],h[i])<mul(h[i+1],h[now+1],h[i]))可以随意变号,就是不要加等于号!因为等于的东西实际上是等效的,所以不用加。加了反而会T。#include<cstdio>#include<c...原创 2019-05-08 09:42:22 · 302 阅读 · 0 评论 -
GCD[caioj1280 莫反]
欢迎大家访问我的老师的OJ———caioj.cn前情提要建议先去看看Zap,数据加强了,但思路较为简单。题面描述传送门思路与Zap不同的就是多了注意:2,3和3,2是一种情况那么我们就要去重了。设D(a,b,k)D(a,b,k)D(a,b,k)表示满足x≤a,y≤bx\le a,y\le bx≤a,y≤b且k∣gcd(x,y)k \mid \gcd(x,y)k∣gcd(x,y)...原创 2019-04-29 20:07:28 · 557 阅读 · 0 评论 -
同余方程[NOIP2012/CH3301]
题面描述传送门思路详细参考#include<cstdio>#include<cstring>#include<algorithm>#include<cmath>#include<cstdlib>#define ll long longusing namespace std;ll exgcd(ll a,ll b,ll ...原创 2019-04-25 10:33:13 · 134 阅读 · 0 评论 -
Sumdiv[POJ1845]
欢迎大家访问我的老师的OJ———caioj.cn题面描述传送门思路一句话题意:求ABA^BAB的所有约数之和mod9901(1≤A,B≤5∗107)\operatorname{mod}9901(1\le A,B\le 5*10^7)mod9901(1≤A,B≤5∗107)提要:等比数列求和公式:Sn=a1∗1−qn1−qS_n=a_1*\frac{1-q^n}{1-q}Sn=a1...原创 2019-04-25 09:23:32 · 102 阅读 · 0 评论 -
The Luckiest Number[POJ3696]
题面描述传送门思路一句话题意:求至少多少个888连在一次组成的正整数是L\operatorname{L}L的倍数。∵x个8\because x个8∵x个8实际上就是(10x−1)/9∗8(10^x-1)/9*8(10x−1)/9∗8,现在我们可以转化一下题意:求最小的xxx,满足L∣(10x−1)/9∗8\operatorname{L} \mid (10^x-1)/9*8L∣(10...原创 2019-04-25 07:38:03 · 118 阅读 · 0 评论 -
[计算几何]Wall[POJ1113]/[UVA1303]
欢迎大家访问我的老师的OJ———caioj.cn题面描述POJUVA思路简要题意大概就是:给定n,l, 然后给定平面中n个点的横纵坐标,求一个最小周长的封闭图形(可以包含曲线),要求该图形包裹住所有点,并且该封闭图形上的任何点距离任意给定点不小于l输出这个最小周长就以样例为例吧:很明显,这个封闭图形的最小周长就为凸包周长加2πl2\pi l2πl。AC code此...原创 2019-05-08 13:35:09 · 220 阅读 · 0 评论 -
NIM博弈证明
古代的Nim\operatorname{Nim}Nim取石子游戏是由两个人面对若干堆石子进行的游戏。设有n≥2n\ge 2n≥2堆石子,各堆分别含有A1、A2⋯AnA_1、A_2\cdots A_nA1、A2⋯An个石子。游戏的目的就是选取最后剩下的石子。游戏规则如下:游戏人交替进行游戏当轮到每个游戏人取子时,选择这些石子中的一堆,并从所选的堆中取走至少一个石子(游戏人也可以取...原创 2019-05-05 19:31:55 · 504 阅读 · 0 评论 -
博弈SG(模版):Nim取石子游戏[caioj1166]
题面描述传送门思路我们先得到sum=A1xorA2xorA3xorA4xor⋯xorAnsum=A_1\operatorname{xor}A_2\operatorname{xor}A_3\operatorname{xor}A_4\operatorname{xor}\cdots\operatorname{xor}A_nsum=A1xorA2xorA3xorA4xor⋯xor...原创 2019-05-05 19:39:51 · 294 阅读 · 2 评论 -
土地划分[fzu1015]
题面描述传送门思路线段是否有交点,之后答案为交点数+原有点数。注意,这里由于相连线段之间有交点,所以与跨立实验的有些不同。AC code#include<cstdio>#include<cstring>#include<algorithm>#include<cmath>#include<cstdlib>using ...原创 2019-05-30 13:42:22 · 168 阅读 · 0 评论 -
The Fortified Forest[POJ1873]
题面描述传送门思路#include<cstdio>#include<cstring>#include<cmath>#include<algorithm>#include<cstdlib>using namespace std;const int N=17;const double inf=1e9;struct nod...原创 2019-05-17 13:20:30 · 266 阅读 · 0 评论 -
Cows[POJ3348]
题目描述传送门思路水题直接算凸包面积/50,向下取整。AC code#include<cstdlib>#include<algorithm>#include<cmath>#include<cstdio>#include<cstring>using namespace std;const int N=1e4+10;...原创 2019-05-13 13:11:12 · 114 阅读 · 0 评论 -
[计算几何]Grandpa's Estate[POJ1228]
欢迎大家访问我的老师的OJ———caioj.cn题目描述网上代码大部分都过不了这个数据461 1 2 2 3 3 4 4 5 5 6 611 180 0 0 1 0 2 0 3 1 1 2 2 3 3 2 380 0 0 2 1 0 1 2 2 2 3 0 3 1 3 2好题可惜数据太...原创 2019-05-10 13:30:05 · 148 阅读 · 0 评论 -
[计算几何]求半平面交的面积
题面描述传送门思路先看一下SCY的视频吧,讲得真的很清(luo)楚(suo)。忽略掉一些奇奇怪怪的语言首先我们先来了解一下SCY开头说的$$原创 2019-05-07 09:54:30 · 488 阅读 · 1 评论 -
[计算几何]Scrambled Polygon[POJ2007]
欢迎大家访问我的老师的OJ———caioj.cn题面描述传送门思路停课集训(真的不想回课内)最后一篇题解了。这道题直接按叉积小优先排序就行了,(0,0)(0,0)(0,0)为基准点。注意:排序时注意若叉积相同,判断距离,距离小的更优先。AC code#include<cstdio>#include<cstring>#include<algorith...原创 2019-05-09 14:02:43 · 141 阅读 · 0 评论 -
[计算几何]凸包
题面描述传送门思路SCY的方法真鬼畜sort看一会才能懂。具体参照一下SCY的视频,此处省略。继续%wjyyy的凸包将一些细节再理一理吧。首先把所有点以横坐标为第一关键字,纵坐标为第二关键字排序。这样可以保证p1p_1p1一定是凸包的左下角的那个点,且一定要被选中。这种做法,需要进行两次操作。第一次,先弄出下凸壳(我的理解:就是从左下点p1p_1p1,到右上点pnp_n...原创 2019-05-06 13:30:50 · 142 阅读 · 0 评论 -
[计算几何]面积
题面描述传送门思路吹爆wjyyy三角剖分求面积这里其实就运用了向量的外积的大小为平行四边形的面积。但是有一点要修改的,我们以a1a_1a1为原点,所以我们mul(a1,a2,a1)mul(a_1,a_2,a_1)mul(a1,a2,a1)并不用计算,从333开始计算。而且博文中描述最终面积仍需要除以2才为答案。AC code#include<cstdio>#...原创 2019-05-06 10:40:16 · 240 阅读 · 0 评论 -
[计算几何]判断线段相交(跨立实验)
欢迎大家访问我的老师的OJ———caioj.cn题面描述传送门思考看完视频之后,发现SCY讲得真(mo)好(hu)。下面我来总结一下规律:我们以p1p_1p1为原点,判断p3p_3p3到p2p_2p2,p2p_2p2到p4p_4p4是否是同一个方向旋转(注意是p3p_3p3到p2p_2p2,p2p_2p2到p4p_4p4),可以利用叉乘来理解。以p3p_3p3...原创 2019-05-06 10:06:00 · 367 阅读 · 0 评论 -
正方形计数[luoguP1665]/caioj[1221]
欢迎大家访问我的老师的OJ———caioj.cn题面描述弱化传送门传送门思路已知AAA、HHH点。易证△ABC≌△EBD\vartriangle ABC≌\vartriangle EBD△ABC≌△EBD。之后我们求对角线交点BBB坐标,用交点坐标求EEE点坐标解即可。上面那个点类似。#include<cstdio>#include<cstring>...原创 2019-05-05 21:04:57 · 400 阅读 · 0 评论 -
[中山市选2011]完全平方数[BZOJ2440]
题面描述传送门思路很明显地,一个数不属于完全平方数或完全平方数的倍数,它的μ\muμ一定不等于0.因此我们只要二分枚举xxx,计算[1,x][1,x][1,x]之中的非完全平方数根据容斥原理:xxx以内非完全平凡数的个数为:xxx-[1,x][1,x][1,x]中一个质数(2,3,5,...2,3,5,...2,3,5,...)的平方的倍数的数量+[1,x][1,x][1,x]中两个...原创 2019-04-30 12:51:57 · 168 阅读 · 0 评论 -
高中数学必备知识
∑i=1n\huge\sum_{i=1}^ni=1∑n即int ans=0;for(int i=1;i<=n;i++)ans+=i;∏i=1n\Huge\prod_{i=1}^ni=1∏n即int ans=0;for(int i=1;i<=n;i++)ans*=i;x∣y x整除y,即...原创 2019-04-19 19:53:02 · 415 阅读 · 4 评论 -
Devu and Flowes[Codeforces451E]
题面描述传送门(由于CF要翻墙,所以推荐用Luogu)思路经过观察,就是一道多重集的组合数的模板题CN+M−1N−1−∑i=1NCN+M−ni−2+∑1≤i<j≤NCN+M−ni−nj−3N−1−⋯+(−1)NCN+M−∑i=1Nni−(N+1)NC_{N+M-1}^{N-1}-\sum_{i=1}^NC_{N+M-n_i-2}+\sum_{1\le i<j\l...原创 2019-04-29 08:52:44 · 96 阅读 · 0 评论 -
开关问题[POJ1830]
题面描述传送门思路设xix_ixi表示第i个开关的操作情况,xi=1x_i=1xi=1表示按了这个开关,xi=0x_i=0xi=0表示没按。再统计ai,ja_{i,j}ai,j表示第i个开关和第j个开关的联系情况,ai,j=1a_{i,j}=1ai,j=1表示按下j会影响i的状态,ai,j=0a_{i,j}=0ai,j=0表示不会影响,特别地,令ai,i=1a_{i,i}=1ai...原创 2019-04-26 11:07:02 · 242 阅读 · 0 评论 -
欧拉定理&费马小定理
若gcd(a,n)gcd(a,n)gcd(a,n)=1,,则aφ(n)=1(mod n)a^{\varphi(n)}=1(mod~n)aφ(n)=1(mod n)铺垫:同余类与剩余系证明:对于nnn的一个简化剩余系{a1‾,a2‾,......aφ(n)‾\overline {a_1},\overline {a_2},......\overline {a_{\varphi(...原创 2019-04-17 14:54:13 · 3427 阅读 · 4 评论 -
球形空间产生器[BZOJ1013]
题面描述传送门思路因为一个球体的球面上的所有点到圆心的距离相等,因此只需求出一个点(x1,x2,x3,......,xnx_1,x_2,x_3,......,x_nx1,x2,x3,......,xn),使得:∑j=1n(ai,j−xj)2=C(i∈[1,n+1])\large\sum_{j=1}^{n}(a_{i,j}-x_j)^2=C (i\in[1,n+1])j=1∑n(a...原创 2019-04-26 08:50:29 · 84 阅读 · 0 评论 -
同余类与剩余系[caioj数论]
剩余系就是指对于某一个特定的正整数nnn,N\mathbb{N}N集中的元素模nnn所得的余数域。完全剩余系对于正整数mmm,有r1,r2,r3,......,rm−2,rm1,rmr_1,r_2,r_3,......,r_{m-2},r_{m_1},r_{m}r1,r2,r3,......,rm−2,rm1,rm个整数元素,且两两模m不同余,这些元素所组成的集合叫做模m的一个...原创 2019-04-17 12:51:40 · 1474 阅读 · 0 评论 -
[欧拉函数]原根
在数论,特别是整除理论中,原根是一个很重要的概念。对于a,ma,ma,m,若aaa与mmm互质,即gcd(a,m)=1gcd(a,m)=1gcd(a,m)=1欧拉定理:若正整数a,n互质,则aφ(n)=1(mod m)a^{\varphi(n)}=1(mod~ m)aφ(n)=1(mod m)使ad=1(mod m)a^d=1(mod~m)ad=1(mod&nb...原创 2019-04-17 11:42:25 · 400 阅读 · 0 评论 -
线性筛
代码#include<cstdio>#include<cstring>#include<algorithm>using namespace std;typedef long long ll;ll prime[1110000],pr=0;bool v[19100000];void getprime(){ memset(v,true,siz...原创 2019-04-17 10:31:56 · 77 阅读 · 0 评论 -
欧拉函数[caioj1158]
111~NNN中与NNN互质的数的个数被成为欧拉函数,记为φ(N)\varphi(N)φ(N)。在算数基本定理中,N=p1c1p2c2...pmcmN=p_1^{c_1}p_2^{c_2}...p_m^{c_m}N=p1c1p2c2...pmcm则:φ(N)=N∗(p1−1p1)∗(p2−1p2)∗......∗(pm−1pm)=N∗∏质数p∣N(1−1p)\varphi(N)=N...原创 2019-04-17 09:12:35 · 223 阅读 · 0 评论 -
同余方程组[caioj1155]
题面描述同余方程是这样的:已知a,b,n,求x的最小正整数解,使得ax=b(mod m)同余方程组是这样:也是求x的最小正整数解,但已知b数组和m数组的情况下,x=b[1](mod m[1]),x=b[2](mod m[2]),x=b[3](mod m[3]),……x=b[n](mod m[n])【输入格式】一行一个整数 n(1<=n<=?)下来n行每行两个整数b[...原创 2019-04-16 20:45:54 · 307 阅读 · 0 评论 -
同余方程[caioj1154]
题面描述已知a,b,m,求x的最小正整数解,使得ax=b(mod m)【输入格式】一行三个整数 a,b,m。 1 ≤ a,b,m ≤ 10^9【输出格式】一行一个整数x,无解输出"no solution!"【样例输入】2 5 7【样例输出】6思路从题目的ax=b(mod m)ax=b(mod~m)ax=b(mod m),可以转化为ax+my=bax+my=...原创 2019-04-16 19:02:47 · 188 阅读 · 0 评论 -
exgcd(扩展欧几里德算法)
题面描述解不定方程Ax+By=K(得到的x和y只是其中一组解)给出A、B、K,求出x和y,满足Ax+By=K。【输入格式】一行三个整数 A,B,K。 1 ≤ A,B,K ≤ 10^9【输出格式】一行两个整数 x,y。如果无解,输出"no solution!"思路令ax+by=gcd(a,b)令ax+by=gcd(a,b)令ax+by=gcd(a,b)∵gcd(b,a m...原创 2019-04-16 17:01:02 · 11824 阅读 · 3 评论 -
欧几里得gcd证明
欧几里得算法:∀a,b∈N,b≠0,gcd(a,b)=gcd(a,mod b)\forall a,b\in N,b\neq0,gcd(a,b)=gcd(a,mod \ b)∀a,b∈N,b̸=0,gcd(a,b)=gcd(a,mod b)若a<ba<ba<b,则gcd(a,b)=gcd(b,a mod b)gcd(a,...原创 2019-04-16 15:42:59 · 266 阅读 · 0 评论 -
[欧拉函数]求原根个数
题面描述【输入格式】第一行n(1<=n<=10000)下来n行,每行一个整数p,p保证为质数,表示求p的原根的数目。(1<=p<=2000 0000)【输出格式】每次询问输出一行一个整数,即x的原根的数目。【样例输入】3233179【样例输出】10824思路原根个数是phi[phi[i]]个,我也不知道为什么?代码#include<...原创 2019-04-17 16:42:58 · 1387 阅读 · 0 评论