Cod sursa(job #970596)

Utilizator miu_mik93FMI - Paduraru Miruna miu_mik93 Data 7 iulie 2013 14:17:51
Problema Suma si numarul divizorilor Scor 30
Compilator cpp Status done
Runda Arhiva educationala Marime 1.3 kb
#include <stdio.h>
#include <math.h>
#include <iostream>
#include <string.h>
#include <stdlib.h>
using namespace std;
#define NMax 1000000
long long int v[NMax];
void ciur()
{
	for (long long int i=2; i<=NMax; i++)
	{
		if (v[i] == 0)
		{
			for (long long int j=i+i; j<=NMax; j+=i)
				v[j] = 1;
		}
	}
}
void SumaSiNrDiv(long long int n, long long int &s, int &nrDiv)
{
	int prim = 1;
	long long int copie_n = n;
	for (int i=2; i<=sqrt(n); i++)
	{
		int p = 1;
		if (!v[i] && !(n % i))
		{
			int actsum = 1;
			while (copie_n % i == 0)
			{
				actsum =  (actsum * i);
				copie_n = copie_n / i;
				p++;
			}
			if (v[copie_n] == 0 && copie_n > sqrt (n) )
			{
				nrDiv *= 2;
				s = (s *( (pow(copie_n, 2)-1)/ (copie_n-1)));
				s = s % 9973;
			}
			actsum = (actsum * i);
			s = (s *((actsum - 1) / (i-1))); s = s % 9973;
			prim = 0;
		}
		nrDiv *= p;
	}
	if (prim)
	{
		nrDiv = 2;
		s = (1 + n) % 9973;
	}
}
int main()
{
	FILE *f = fopen("ssnd.in", "r");
	FILE *g = fopen("ssnd.out", "w");
	int t; long long int n; 
	fscanf(f, "%d", &t);
	ciur();
	for (int i=0; i<t; i++)
	{
		long long int s = 1;
		int nrDiv = 1;
		fscanf(f, "%lld", &n);
		SumaSiNrDiv(n, s, nrDiv);
		fprintf(g, "%d %lld \n", nrDiv, s%9973);
	}
	fclose(f);
	fclose(g);
	return 0;
}