摘要：题意：题意：给出n和m，求满足条件gcd(x, n)>=m的x的gcd(x, n)的和，其中1 2 #include 3 #define ll long long 4 using namespace std; 5 6 /* 7 int euler(int n){ 8 int ans=1; 9 for(int i=2; i*i1){ 20 ans*... 阅读全文

摘要：euler(x)公式能计算小于等于x的并且和x互质的数的个数； 我们再看一下如何求小于等于n的和n互质的数的和， 我们用sum(n)表示； 若gcd(x, a)=1，则有gcd(x, x-a)=1； 证明：假设gcd(x, x-a)=k (k>1)，那么有(x-a)%k=0 1式，x%k=0 2式； 阅读全文

摘要：题目链接：http://acm.nyist.net/JudgeOnline/problem.php?pid=775 关于扩展欧几里得算法还是推一遍好啦； 有方程：a*x+b*y=d=gcd(a, b) 1式（只要a, b不全为0则此方程必有解，不过我不会证明，望大神路过时教一下）； 又有gcd(a, 阅读全文

摘要：题目：noyj774 用代数余子式求逆矩阵方法： 若现有矩阵A，要求其逆矩阵； 若|A|==0，则其不存在逆矩阵； 若|A|！=0，其逆矩阵A^-1==*A/|A|；其中*A为其伴随矩阵； 伴随矩阵的求法： *A[j][i]==|M[i][j]|，其中M[i][j]为A[i][j]的代数余子式； 即 阅读全文

摘要：//*************快速幂计算**************************************** 朴素算法实现： 快速幂算法： 原理： 二分： 假设我们现在要计算pow(x,n），那么有当n为偶数时pow(x, n)==pow(x*x, n/2)，当n为奇数时，pow(x, 阅读全文

摘要：题意：求R-L区间满足x=a1*k+b1=a2*l+b2的x的个数; 思路：求出最小的满足条件的x0，则ans=(L-x)/(a1/gcd(a1, a2)*a2)+1; 注意剪枝，不然会超时； 代码： 据说这题应该用拓展欧里几德解。。可惜我看了好久也没弄懂。。诶。。继续看吧。。。。 阅读全文