Cod sursa(job #1016650)

Utilizator Dddarius95Darius-Florentin Neatu Dddarius95 Data 26 octombrie 2013 16:03:46
Problema Suma divizorilor Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 0.95 kb
#include <fstream>
#define X 9901
using namespace std;
ifstream f("sumdiv.in");
ofstream g("sumdiv.out");

int A,B;
long long S;


long long Putere(long long N,long long P)
{

    long long m=1;
    while (P!=1)
    if(P % 2==0)
    {
        N=(N*N) % X;
        P/=2;
    }
    else
    {
        m=(m*N) % X;
        P--;
    }
    return (m*N)% X;
}



int main()
{
    f>>A>>B;
    if(A % 2==0)
    {
        int k=0;
        while(A % 2==0)A/=2,++k;
        int p1=(Putere(2,k*B+1)-1)% X;
        S=(1LL*S*p1 % X);
    }
    for(int i=3;i*i<=A;i+=2)
    if(A % i==0)
    {
        int k=0;
        while(A % i==0)A/=i,++k;
        int p1=(Putere(i,k*B+1)-1)% X;
        int p2=Putere(i-1,X-2);
        S=(((1LL*S*p1) % X)*p2)% X;
    }
    if(A>1)
    {
        int p1=(Putere(A,B+1)-1)% X;
        int p2=Putere(A-1,X-2);
        S=(((1LL*S*p1) % X)*p2)% X;
    }
    f.close();g.close();
    return 0;
}