Cod sursa(job #1820692)

Utilizator Kln1000Ciobanu Bogdan Kln1000 Data 2 decembrie 2016 01:40:06
Problema Suma si numarul divizorilor Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 2.04 kb
#include <fstream>
#include <vector>
#include <bitset>
#define mod 0x26F5
#define MXP 0xF424A
std::vector <int> prim;
std::ofstream t ("ssnd.out");
class parser{
    public:
        parser() {}
        parser(const char *file_name){
            input_file.open(file_name);
            cursor=0;
            input_file.read(buffer,SIZE);}
        inline parser &operator >>(int64_t &n){
            while(buffer[cursor]<'0' or buffer[cursor]>'9') inc();
            n=0;
            while('0'<=buffer[cursor] and buffer[cursor]<='9')
                n=n*10+buffer[cursor]-'0',inc();
            return *this;}
    private:
        std::ifstream input_file;
        static const int SIZE=0x8800;
        char buffer[SIZE];
        int cursor=0;
        inline void inc(){
            if(++cursor==SIZE)
                cursor=0,input_file.read(buffer,SIZE);
        }
};


int64_t pow(int64_t x,int64_t y)
{
    x%=mod;
    int64_t m=1;
    while(y){
        if(y&1)
            m=(m*x)%mod;
        x=(x*x)%mod;
        y>>=1;
    }
    return m;
}

void ciuruieste()
{
    std::bitset <MXP> prime;
    for (int i=1;((i*i)<<1)+(i<<1)<MXP;++i)
        if(!prime[i])
            for(int j=((i*i)<<1)+(i<<1);(j<<1)+1<MXP;j+=(i<<1)+1)
                prime[j]=true;
    prim.push_back(2);
    for(int i=1;(i<<1)+1<MXP;++i)
        if(!prime[i])
            prim.push_back((i<<1)+1);
}

void solve(int64_t target){
    int how_many=1,sum=1,aux;
    for (std::vector<int>::iterator i=prim.begin();(*i)*(*i)<=target and i!=prim.end();advance(i,1)){
        if (target%(*i)==0){aux=0;
            while (target%(*i)==0){target/=(*i);++aux;}
        how_many*=(1+aux);
        sum=(((sum*((pow((*i),1+aux)-1)%mod))%mod)*(pow((*i)-1,mod-2)%mod))%mod;}
    }
    if (target>1){
        how_many<<=1;
        sum=(sum*(1+target))%mod;}
    t<<how_many<<" "<<sum<<'\n';
}

int main()
{
    parser f ("ssnd.in");
    int64_t aux,q;
    ciuruieste();
    f>>q;
    for (;q;--q)
        f>>aux,solve(aux);
    return 0;
}