Cod sursa(job #970713)

Utilizator miu_mik93FMI - Paduraru Miruna miu_mik93 Data 7 iulie 2013 16:52:17
Problema Suma si numarul divizorilor Scor 70
Compilator cpp Status done
Runda Arhiva educationala Marime 1.6 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+1];
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;
		}
	}
}
int isPrim(long long int n)
{
	for (int i=2; i<=sqrt(n); i++)
		if( n % i == 0)
			return 0;
	return 1;
}
void SumaSiNrDiv(long long int n, long int &s, long 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 (copie_n >= NMax && isPrim(copie_n) && copie_n > sqrt(n))
			{
				nrDiv *= 2;
				s = (s *( (pow(copie_n, 2)-1)/ (copie_n-1)));
				s = s % 9973;
			}
			else if (copie_n < NMax && 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 int s = 1, nrDiv = 1, s1 = 1, nrDiv1 = 1;
		fscanf(f, "%lld", &n);
		SumaSiNrDiv(n, s, nrDiv);
		fprintf(g, "%ld %ld \n", nrDiv, s%9973);
	}
	fclose(f);
	fclose(g);
	return 0;
}