摘要：
#include #include #include #include using namespace std;const int MaxM=11;int a[MaxM],b[MaxM];void exgcd(int a,int b,int &d,int &x,int &y){ if(b==0){ ... 阅读全文
摘要：
特判r1=0时的情况，因为0是不能模的。#include #include #include #include using namespace std;const int MaxM=11;int a[MaxM],b[MaxM];void exgcd(int a,int b,int &d,int &x... 阅读全文
摘要：
和POJ 1061一样要求最小解，就尽可能的把ax的值附到by上去，所以可以有ax=b*k+a*v（因为附到by上后必须仍上a*x的形式）。两边同除a就可得到结果。但是，我们知道，（a,b）=1。所以k|a，也就是说，ans=(x%b+b)%b。后来加上b是为了防止负数。#include #incl... 阅读全文
摘要：
线性同余方程组，模板了。但要注意读完数据才跳出循环啊#include #include #include #include using namespace std;void exgcd(__int64 a,__int64 b,__int64 &d,__int64 &xx,__int64 &yy){ ... 阅读全文
摘要：
#include #include #include using namespace std;int main(){ int a,b; while(scanf("%d%d",&a,&b)!=EOF){ if(!a&&!b) break; int ans=1; for(int i=1;i<=b;... 阅读全文
摘要：
注释掉那段会超时，呃。。#include #include #include #include #include using namespace std;const int Max=1000000;bool val[Max+1];int p[350];int main(){ int t,n; sca... 阅读全文
摘要：
由算术基本定理，直接使用公式就好#include #include #include #include using namespace std;const int Maxp=1000;bool isprime[Maxp];int prime[Maxp],nprime;void Doprime(){ ... 阅读全文
摘要：
很容易由算术基本定理知道，完全平方数就是所有质因子指数为偶数的数。而求得N以下的质因子，可由前两篇的公式知，由N!与p的关系求得。对于指数为p的，用N!除去就可以，因为p必定属于N以内，且无重复。至于除法，在下实在不会，学得别人的，记录一下。MOD数除法，可以由费马小定理a^(p-1)=1 (mod... 阅读全文