Cod sursa(job #2496529)

Utilizator divianegoescuDivia Negoescu divianegoescu Data 20 noiembrie 2019 23:16:56
Problema Suma si numarul divizorilor Scor 70
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.29 kb
#include <fstream>
#include <bitset>
#define R 9973
using namespace std;
ifstream fin("ssnd.in");
ofstream fout("ssnd.out");
long long n,t,i,j,exp,prime[100000];
long long x,nrdiv,divizor,putere,sumdiv;
long long explog(long long a,long long p){
    if(p==0)return 1;
    long long r=explog(a,p/2);
    if(p%2)return r*r*a;
    return r*r;
}
bitset <1000100> b;
int main(){
     for(i=2;i<=1000000;i++)
        if(b[i]==0){
            prime[++prime[0]]=i;
            for(j=i+i;j<=1000000;j+=i)
                b[j]=1;
        }
     for(fin>>t;t--;){
        fin>>x;
        if(x==1){fout<<"1 1\n";continue;}
        nrdiv=sumdiv=1;
        for(i=1;i<=prime[0] && prime[i]<=x/prime[i];i++){
            if(x%prime[i])continue;
            exp=0;
            divizor=prime[i];
            while(x%prime[i]==0){
                x/=divizor;
                exp++;
            }
            if(exp==1)putere=divizor*divizor;
            else if(exp==1)putere=divizor*divizor*divizor;
            else putere=explog(divizor,exp+1);
            nrdiv*=(exp+1);
            sumdiv=(sumdiv*(putere-1)/(divizor-1))%R;
        }
        if(x>1){
            nrdiv*=2;
            sumdiv=(sumdiv*(x+1))%R;
        }
        fout<<nrdiv<<" "<<sumdiv<<"\n";
    }
    return 0;
}