Pagini recente » Cod sursa (job #1491851) | Cod sursa (job #3326457)
#include <iostream>
#include <fstream>
using namespace std;
ifstream f("mins.in");
ofstream g("mins.out");
const int NMAX = 5001;
int ciurn[NMAX], ///Vector al numarului de divizori primi
ciurp[NMAX], ///Ciurul numerelor divizibile cu patrate de numere prime
c, d, n; ///(i.e. la care nu toti divizorii primi sunt distincti)
long long card;
void calcul(int i)
{
if(ciurn[i] == 0)
{
for(int j = i; j <= n; j += i)
ciurn[j]++;
int k = i * i;
for(int jj = k; jj <= n; jj += k)
ciurp[jj] = 1;
}
if(ciurp[i] == 0)
{
int t = (c / i) * (d / i);
if(ciurn[i] % 2 == 0)
card -= t;
else
card += t;
}
}
int main()
{
f >> c >> d;
--c, --d;
n = min(c, d);
for(int i = 2; i <= n; i++)
calcul(i);
g << 1LL * c * d - card;
f.close();
g.close();
return 0;
}