Cod sursa(job #2610720)

Utilizator dancu_mihai13Dancu Mihai dancu_mihai13 Data 5 mai 2020 15:51:42
Problema Suma si numarul divizorilor Scor 0
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.3 kb
#include <bits/stdc++.h>

using namespace std;

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

const int NMAX = 1000005;
const int MOD = 9973;
int primes[NMAX], c, d, sum_div, nr_div, power;
bool ciur[NMAX / 2];
long long val;

void init()
{
	for (int i = 2; i * i < NMAX; i++)
		if (!ciur[i])
		{
			primes[++c] = i;
			for (int j = i * i; j < NMAX; j += i)
				ciur[j] = 1;
		}
}

int exp_log(int n, int p)
{
    int r = 1;
    n %= MOD;

    while(p)
    {
        if(p % 2)
            r = (r * n) % MOD;
        n = n * n;
        n %= MOD;
        p /= 2;
    }
    return r;
}

void divizori(long long numar)
{
	nr_div = sum_div = 1;
	int j = 1;
	d = primes[j];
	while (numar > 1)
	{
		power = 0;
		while (!(numar % d))
		{
			numar /= d;
			power++;
		}
		if (power)
		{
			nr_div *= (power + 1);
			sum_div = ((sum_div * (exp_log(d, power + 1) - 1) % MOD) / (d - 1)) % MOD;
		}
		j++;
		d = primes[j];
		if(d * d > numar && numar > 1)
        {
            nr_div *= 2;
            sum_div = (sum_div * (numar + 1) * 1LL) % MOD;
            numar = 1;
        }
	}
	fout << nr_div << ' ' << sum_div << '\n';
}

int main()
{
	int nr; fin >> nr;
	init();
	for (int i = 0; i < nr; i++)
	{
		fin >> val;
		divizori(val);
	}
	return 0;
}