Cod sursa(job #562454)

Utilizator SadmannCornigeanu Calin Sadmann Data 23 martie 2011 08:44:10
Problema Suma si numarul divizorilor Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.53 kb
#include<stdio.h>
#include<math.h>
#include<bitset>
#define Nmx 1<<17
using namespace std;

const int MAX_N = 1000005;
const int MOD = 9973;

FILE *in,*out;

long long N;
int T, K=1, P[MAX_N];
unsigned char V[Nmx];
/*bitset<MAX_N> viz;

void ciur() {
	for(int i = 2; i < MAX_N; ++i) {
		if(viz[i] == 0) {
			P[++K] = i;

			for(int j = i+i; j < MAX_N; j += i) {
				viz[j] = 1;
			}
		}
	}
}*/

void ciur()
{
    int i2,j;
    for(int i=3;i<MAX_N;i+=2)
    {
        if(V[i>>4] & (1<<( (i>>1) & 7) ) )
            continue;
        P[++K]=i;
        for(j=i+(i2=i+i);j<MAX_N;j+=i2)
            V[j>>4]|= 1<<( (j>>1) & 7);
    }
}

inline int putere(int x, int p) {
	int rez = 1;
	x %= MOD;
    int q;
    for (q = 0; (1<<q) <= p; ++ q)
    {
        if( ((1<<q) & p) )
            rez=(rez*x)%MOD;
        x=(x*x)%MOD;
    }
	return rez;
}

void solve() {
	fscanf(in,"%lld",&N);

	int	nd = 1, sd = 1;

	for(int i = 1; i <= K && 1LL * P[i] * P[i] <= N; ++i) {
		if(N % P[i]) continue;
		int p = 0;

		while(N % P[i] == 0)
		{
			N /= P[i];
			++p;
		}

		nd *= (p+1);

		int p1 = (putere(P[i], p+1) - 1) % MOD;
		int p2 = putere(P[i]-1, MOD-2) % MOD;

		sd = (((sd * p1) % MOD) * p2) % MOD;
	}

	if(N > 1) {
		nd *= 2;
		sd = (1LL*sd*(N+1) % MOD);
	}

	fprintf(out,"%d %d\n",nd,sd);
}

int main() {
    P[1]=2;
	ciur();
    in=fopen("ssnd.in","rt");
    out=fopen("ssnd.out","wt");
    fscanf(in,"%d",&T);
	for(; T; --T)
		solve();
    fclose(in);
    fclose(out);
}