数论模板

const int maxn=1025;
const int maxp=33;//sqrt
const int maxl=11;//log
bool isprime[maxp];
int prime[20],c[maxn];

void makeprime()
{
int i,j;
memset(isprime,1,sizeof(isprime));
for(i=2;i<=maxn;i++)
{
if(isprime[i])
prime[tot++]=i;
for(j=0;j<tot&&i*prime[j]<=maxn;j++)
{
isprime[i*prime[j]]=false;
if(i%prime[j]==0)break;
}
}
}

int power(int x,int n)
{
int ret=1;
while(n--)ret*=x;
return ret;
}

void abundance()
{
int n,m,i,j,k,sum;
memset(c,-1,sizeof(c));
for(n=1;n<maxn;n++)
{
k=n;
sum=1;
for(i=0;i<tot&&prime[i]*prime[i]<=k;i++)
if(k%prime[i]==0)
{
m=1;
while(k%prime[i]==0)
{
k/=prime[i];
m++;
}
sum*=(power(prime[i],m)-1)/(prime[i]-1);
}
if(k>1)sum*=k+1;
if(sum>n*2)c
=sum-n*2;
}
}

void process()
{
int i,j,l;
for(i=0;i<maxn;i++)
h[i][0]=c[i];
for(j=1,l=1;2*l<=maxn;j++,l*=2)
for(i=0;i<=maxn-l*2;i++)
h[i][j]=max(h[i][j-1],h[i+1][j-1]);
}

int query(int be,int ed)
{
int j=0;
int l=1;
while(l*2<=ed-be+1)
{
j++;
l*=2;
}
return (h[be][j]>h[ed+1-l][j]?h[be][j]:h[ed+1-l][j]);
}