Cod sursa(job #883600)

Utilizator RobertSSamoilescu Robert RobertS Data 20 februarie 2013 10:35:04
Problema Suma si numarul divizorilor Scor 40
Compilator cpp Status done
Runda Arhiva educationala Marime 1.31 kb
#include<iostream>
#include<fstream>
#include<vector>
using namespace std;
const int MAX_N = 1000005;
const int MOD = 9973;

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

struct pereche{
	int div, nr;
	pereche(){
	 nr = 0;
	}
};

int t;
bool prim[MAX_N];
vector<long>vec;

//---==== PENTRU DETERMINAREA NUMERELOR PRIME ===----//
void ciur(){
	prim[1] = true;
	for(int i=2; i<=MAX_N; i++){
		if(!prim[i]){
			vec.push_back(i);
			for(int j=i; j<=MAX_N; j+=i)
				prim[j] = true; 
		}
	}
}


inline int pow(int x, int p) {
    int rez = 1;
    x %= MOD;
 
    for(; p; p >>= 1) {
        if(p & 1) {
            rez *= x;
            rez %= MOD;
        }
        x *= x;
        x %= MOD;
    }
 
    return rez;
}


void solve(int n){
	long nr_div = 1;
	long suma_div = 1;
	for(size_t i=0; i<vec.size() && n>1; i++){
		if(n%vec[i] == 0){
			pereche p; p.div = vec[i];
			while(n%vec[i] == 0){
				p.nr ++;
				n/= vec[i];
			}
			
			int p1 = (pow(p.div, p.nr+1)-1)%MOD;
			int p2 = pow(p.div-1, MOD-2)%MOD;
			
			nr_div *=(p.nr+1);
			suma_div *= (p1*p2)%MOD;
			
		}
	}
	
	
	fout << nr_div << " " <<suma_div % MOD<< '\n';
}


int main(){

	
	ciur();
	fin >> t;
	for(int i=1; i<=t; i++){
		long n; fin >> n;
		solve(n);
	}
	

return 0;
}