摘要：/* 设a≡b[i](mod w[i) i=1->n; 我们记k=w[1]*w[2]*w[3]*……*w[n] 则方程组在mod k 同余的意义下有唯一解 我们令x=(k/w[i])*yi 那么方程等价于(k/w[i])y≡1(mod w[i]) 那么方程组的解x0=b1x1+b2x2+...+bn 阅读全文

摘要：若a*x≡1(mod b) a，b互质，则称x为a的逆元，记为a-1。在计算(t/a)mod b 时等价于t*a-1mod b。 显然可转化为ax+by=1,用扩欧来求解，但相对来说是比较慢的，这里还有线性算法: 首先1-1≡1(mod p) 设p=k*i+r，那么k*i+r≡0(mod p)两边同 阅读全文

摘要：φ//如果不熟悉欧拉函数可转超链接<-- 正题： 1.仪仗队： 求从(0,0)点可以看到的点，我们考虑正比例函数的斜率，同一斜率上只能看到一个点，我们要知道对于斜率0~1在一个n*n的点阵中有多少可能的斜率使得有若干点在函数上。 观察规律： 1*1：显然答案为0（自己看到自己当然不算了） 2*2：有 阅读全文

摘要：我们规定φ（p）表示1~p-1中与p互质的数的个数，规定φ（1）=1。 有以下性质： 1.当p为素数是φ（p）=p-1 2.设m>1,(a,m)=1,则： aφ(m)≡1(mod m). (欧拉定理) 3.设p为素数，(a,p)=1，则： ap-1≡1(mod p).(费马小定理) 4.若i mod 阅读全文

摘要：#include #include #include using namespace std; bool b[1000010]; int prime[100010],point; int n; int main() { scanf("%d",&n); b[0]=b[1]=1; for(int i=2;in)break; b[i*prime[j]]=1; if(i%prime[j... 阅读全文