Cod sursa(job #978754)

Utilizator andrettiAndretti Naiden andretti Data 29 iulie 2013 17:03:04
Problema Suma si numarul divizorilor Scor 70
Compilator cpp Status done
Runda Arhiva educationala Marime 1.75 kb
#include<fstream>
#include<math.h>
#include<cstring>
#define maxd 1000005
#define maxp 80000
#define mod 9973
using namespace std;
typedef long long ll;

ifstream fin("ssnd.in");
ofstream fout("ssnd.out");

int t,p,nrd,sum;
int prime[maxp];
ll n;
bool sieve[maxd];

void preproc()
{
    int i,lim=maxd-5;
    prime[++p]=2;
    for(i=3;i*i<=lim;i+=2)
      if(!sieve[i])
      {
          prime[++p]=i;
          for(int j=i;i*j<=lim;j++)
           sieve[i*j]=true;
      }
    for(;i<=lim;i+=2)
     if(!sieve[i])
      prime[++p]=i;
}

void gcd(int a,int b,int &x,int &y)
{
    if(b==0) {x=1; y=0; return;}
    int x0,y0;
    gcd(b,a%b,x0,y0);
    x=y0%mod;
    y=(x0-a/b*y0)%mod;
}

int modular(int a,int b)
{
    int x,y;
    gcd(a,b,x,y);
    while(x<0) x+=mod;
    return x;
}

void update(int nr,int prod,int x)
{
    nrd*=(nr+1);
    sum=(sum*((prod-1)%mod))%mod;
    sum=(sum*x)%mod;
}

void solve()
{
    int nr,prod,inv;
    ll sq;

    fin>>t;
    for(int i=1;i<=t;i++)
    {
        fin>>n;
        sq=sqrt(n); nrd=1; sum=1;

        for(int j=1;prime[j]<=sq;j++)
         if(n%prime[j]==0)
         {
             inv=modular(prime[j]-1,mod);

             nr=0; prod=prime[j]%mod;
             while(n%prime[j]==0)
             {
                 nr++;
                 prod=(prod*(prime[j]%mod))%mod;
                 n/=prime[j];
             }
             update(nr,prod,inv);
         }

        if(n>1)
        {
            inv=modular(n-1,mod);
            prod=( (n%mod)*(n%mod) )%mod;
            update(1,prod,inv);
        }

        fout<<nrd<<" "<<sum<<'\n';
    }
}

int main()
{
    preproc();
    solve();

    fclose(stdin);
    fclose(stdout);
    return 0;
}