Pagini recente » Cod sursa (job #2295410) | Cod sursa (job #2036173) | Cod sursa (job #2487596) | Cod sursa (job #1529137) | Cod sursa (job #1677635)
#include <fstream>
#include <bitset>
#include <math.h>
#define NM 1000005
using namespace std;
bitset<NM> pr;
void ciur()
{
pr[1] = 1;
for(int i = 2; i < 1000; ++i)
if(!pr[i])
for(int j = i*i; j < NM-5; j += i)
pr[j] = 1;
}
int main()
{
ifstream f("fractii.in");
ofstream g("fractii.out");
long long nr = 0;
int n;
ciur();
f >> n;
for(int i = 2; i <= n; ++i)
{
long long r = 1;
for(int c = 2; c <= i; ++c)
if(!pr[c])
{
int pw = 0, rp = i;
while(!(rp%c))
{
rp /= c;
++pw;
}
if(pw)
r *= (c-1)*pow(c, pw-1);
}
nr += r;
}
g << nr*2+1;
}