做题目遇到的数论用的代码
杨氏筛法求质数
void IsPrime()
{
int sum = 0;
for(int i = 2; i <= 1000000; ++i)
if(Prime[i] == 0)
{
++sum;
Prime[i] = sum;
for(int j = i+i; j <= 1000000; j+=i)
Prime[j] = sum;
}
}
快速幂
int Mod(int a, int b, int c) {
int ans = 1;
while (b) {
if (b & 1) {
ans = ans * a % c;
}
b >>= 1;
a = a * a % c;
}
return ans;
}
加法代替乘法的快速幂,下面两个合在一起可以算大数值的快速幂。
ll mul(ll a,ll b,ll c)//第一步
{
ll ans=0;
while(b)
{
if(b&1)
{
ans=(ans+a)%c;
}
a=(a<<1)%c;
b>>=1;
}
return ans;
}
ll qsm(ll a,ll b,ll c)
{
ll ans=1;
a=a%c;
while(b!=0)
{
if(b&1) ans=mul(ans,a,c);
b>>=1;
a=mul(a,a,c);
}
return ans;
}
最大公约数
int gcd(int a,int b)
{
return (b == 0)?a:gcd(b,a%b);
}
扩展欧基米德
ll solve(ll a,ll b,ll &x,ll &y)
{
if(b ==0)
{
x=1; y=0;
return a;
}
ll d=solve(b,a%b,x,y);
ll xx=x;
x=y;
y=xx-a/by;
return d;
}
欧拉函数
int euler(int x)
{
int ans=x;
for (int i=2;ii<=x;i++)
{
if (x% ==0)
{
ans=ans/i*(i - 1);
while(x%i ==0)
x/=i;
}
}
if (x>1)
ans=ans/x*(x - 1);
return ans;
}
欧拉函数打表的方法
void Euler(){
phi[1] = 1;
for(int i = 2; i < maxn; i++){
phi[i] = i;
}
for(int i = 2; i < maxn; i++){
if(phi[i] == i){
for(int j = i; j < maxn; j += i){
phi[j] = phi[j] / i * (i - 1);
}
}
}
}
升级版
int vis[maxn];
void init(){
int i,j;
memset(vis, 0, sizeof(vis));
for(i=1;i<maxn;i++){
vis[i]=i;
}
for(i=2;i<maxn;i++){
if(vis[i]==i){
for(j=i;j<maxn;j=j+i){
vis[j]=(vis[j]/i)*(i-1);
}
}
}
}
这周大部分题都是围绕上面几个数论要么变换一下,要么就是从题目推导出这些定理的解题。还有的中国剩余定理,线性同余还是懵,继续努力”专研”一下。总体感觉老师上课没怎么听懂,数论非常烧脑题目围绕定理推来推去推去推出一个简化的公式,看题解都一阵头大,除此之外要么是要优化算法防止时间超时,要么就是注意数据范围大小注意代码。总之,刚开始接触就觉得很难,很头大。