10.22

数论，

筛素数：

bool vis[maxn];
int prime[maxn];
int prim(int n)
{
memset(vis,0,sizeof(vis));
int cnt=0;
for(i=2;i<n;i++)
{
if(!vis[i])
prime[cnt++]=i;
for(j=0;j<cnt&&i*prime[j]<n;j+=i)
{
vis[i*prime[j]=1;
if(i%prime[j]==0)
break;
}
}
return cnt;//小于n的素数的个数
}

欧拉函数:小于等于n的正整数中与n互质的数的个数

int euler(int n)
{
int ret=n;
for(i=2;i<=sqrt(n);i++)
{
if(n%i==0)
ret=ret/i*(i-1);
while(n%i==0)
n/=i;
}
if(n>1)
ret=ret/n*(n-1);
return ret;
}

int euler[maxn];
void eule()
{
    euler[1]=1;
    for(i=2;i<n;i++)
        euler[i]=i;
    for(i=2;i<n;i++)
    {
        if(euler[i]=i)
        {
            for(j=i;j<n;j+=i)
            {
                euler[j]=euler[j]/i*(i-1);
            }
        }
    }
}

欧拉定理，费马定理 Miller-Rabin算法，伪随机数还没弄懂，密码学

根据日期数算星期几有公式

w=(d+2*m+3*(m+1)/5+y+y/4-y/100+y/400)%7;

注意一月份，二月份，是上年的十三月份，十四月份

如何计算a^bmodc

也就是常说的快速幂取模

int quickmd(int a,int b,int c)
{
    a=a%c;
    int cnt=1;
    while(b)
    {
        if(b&1)
            cnt=cnt*a%c;
        a=a*a%c;
        b>>=1;
    }
    return cnt;
}