3.用高级毒把质因子相乘

理论https://blog.nowcoder.net/n/09c423c7181546ba97d257de92f17ce3 #include<bits/stdc++.h> using namespace std;   const double eps=1e-6; const int N=105;   int c,r,n,t,flag; double a[N][N];   int guass(){     for(c=0, r=0;c<n;c++){//c列 r行       &nbsp...

模板代码https://blog.nowcoder.net/n/e2dfbe8bd9844910a43f8e8f43c0769f

给出条件： 一：m1,m2...mk两两互质 二： .... 求x 令 令(与mi互质)（用exgcd求mi的逆） 所以可求 Mi模mi的逆（用exgcd求Mi*infMi=1（mod p）） 则 证明上式：题意有 x=a1（mod m1） 按上式：x=a1*M1*infM1=a1 (mod m1) （其余的含m1被约掉）

a*a-1=1（mod m）只有当a与m互质时a才有逆元 求逆元 如果p是质数，快速幂费马小定理a^p-2 如果p不是质数，求a*a-1=1（mod p），用欧几里得exgcd（a，p，a-1，y）    a和p必须是正数，（用exgcd求逆元时候a-1可能是负数，这个时候强转就行了a-1=（a-1%p+p）%p） 分数mod p 对于（n/m）%mod，应该写为：n*（m的逆）%mod 因为，n/m可能是非整数，不可以直接按顺序计算 逆元的积性 即 预处理逆元 ...

😀在c++中 a mod n= 😀模运算性质 a0=a mod n,    b0=b mod n 😀当k|a时有：

求a*b%n a*b会爆ull typedef unsigned long long ull; ull mul(ull a,ull b,ull n){     ull q=(long double)a*b/n;//long double a*b会记录下a*b的前几位，后面几位有误差，但除n只关心前几位，最后转为ull相当于向下取整   return (a*b-q*n+n)%n;//a*b相当于a*b%ull ...

//a*b%m ll mul(ll a,ll b){ ll s=0; for(;b;b>>=1){   if(b&1) s=(s+a)%m;     if(b-1) a=(a+a)%m; } return s; }

明显最小正多边形在三角形外接圆上，我们只要求圆心角的一小块，再x份数 且最小就是要边数越多（圆心角越小），可以看出多边形就越接近圆（面积越大） 数据会给三点坐标。三条边就出来了 三条边出来了我们就可以算出每条边对应的圆心角 然后就gcd（圆心角）g 圆心角小份面积/ 份数 #include <bits/stdc++.h> using namespace std; const double pi=3.1415926535, ACC=2...

就是选择c 首先由c种选择 时间=wa+最后ac+前中ac wa：5*c*a 最后ac：5*c（a+b+1） 前中ac： 由于前中每个位置挑选没有特殊性，所以每个位置分得的次数是一样的 平均每个位置分得的次数：b/(a+b) 平均每个位置的罚时：（a+b+1）（a+b）/2 前中ac所占的罚时： 所有情况的时间相加/情况数： 所以只要判断是否可以被2除 #include <bits/stdc++.h> using namespace std; const in...

线每次穿过一个格子都会经过一条边 就可以得到要走r+c-1，但是x/y=r/c处经两个边但只穿了一个格子所以 穿过格子数N=R+C-gcd（R,C） 我们现在知道N，求满足该式的解的个数           R和C是肯定可以被gcd（R,C）整除，所以等式右边可以被gcd整除，所以左边——N可以被gcd整除 此时r和c互质，n是N的因数 即gcd（r，c）=1 gcd（r，c）=gcd（r+c，c）=gcd（n+1，c）=1 求 &n...

把这个看成是一个坐标，左下是原点，任何一个（x，y）和原点连线就是斜率k 能看见的人x和y都是互质的（除x，y轴） 所以只要遍历1~n-1求和欧拉函数phi，然后加上3个没加到的点 注意加上三个没加到的是当n>2时候，所以要特判 #include <bits/stdc++.h> using namespace std; const int N=4e4+4; long long n,phi[N],cnt,prime[N],ans; bool st[N]; int&nb...

gcd（a，b）=g a=k1 g        ,        b=k2 g 且gcd（k1，k2）=1 gcd（a，b）max就要是n/（k1+k2）最大 得k1=1，k2=2 即k1*2=k2 #include <iostream> using namespace std; int t,n; int main(int a...