Cod sursa(job #988734)

Utilizator misu007Pogonaru Mihai misu007 Data 23 august 2013 19:05:07
Problema Suma si numarul divizorilor Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 1.55 kb
#include <cstdio>
using namespace std;

int nr=1,a[78498];
char p[62501];

int pow(int x)
{
    return (x*x)%9973;
}

int putlg(int y,int x)
{
    if(y==1) return x;
    else if(y&1) return (pow(putlg(y>>1,x))*x)%9973;
    return pow(putlg(y>>1,x));
}

void ciur(int n)
{
  int i, j;
  for (i = 1; ((i * i) << 1) + (i << 1) <= n; i += 1) {
    if ((p[i >> 3] & (1 << (i & 7))) == 0) {
      for (j = ((i * i) << 1) + (i << 1); (j << 1) + 1 <= n; j += (i << 1) + 1) {
        p[j >> 3] |= (1 << (j & 7));
      }
    }
  }
  for (i = 1; 2 * i + 1 <= n; ++i)
       if ((p[i >> 3] & (1 << (i & 7))) == 0)
           a[nr++]=2*i+1;
  a[0]=2;
}

int main()
{
    freopen("ssnd.in","r",stdin);
    freopen("ssnd.out","w",stdout);
    int i,t,s;
    long long n,k,p;
    ciur(1000000);
    scanf("%d",&t);
    while(t)
    {
        --t;
        scanf("%lld",&n);
        if(n==1) printf("1 1\n");
        else
        {
            nr=1;s=1;
            for(i=0;n>1&&a[i]<=n/2;++i)
            {
                k=0;
                while(n%a[i]==0)
                {
                    ++k;
                    n=n/a[i];
                }
                if(k)
                {
                    nr*=k+1;
                    p=putlg(k+1,a[i]%9973);
                    s=(s*(p-1)/((a[i]-1)%9973))%9973;
                }
            }
            if(n>1)
            {
                nr*=2;
                s=(s*((n-1)%9973))%9973;
            }
            printf("%d %d\n",nr,s);
        }
    }
    return 0;
}