Cod sursa(job #3146695)

Utilizator AlexSerban21Serban Alexandru AlexSerban21 Data 22 august 2023 11:12:59
Problema Fractii Scor 100
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 0.8 kb
#include <fstream>
#include <bitset>
using namespace std;
ifstream fin ("fractii.in");
ofstream fout ("fractii.out");
int k,j,n,i,d,x,v[500001];
long long sum,phi;
bitset <1000001> ciur;
int main()
{
    fin>>n;
    for (i=2; i<=n; i++)
    {
        if (ciur[i]==0)
        {
            v[++k]=i;
            for (j=i+i; j<=n; j+=i)
                ciur[j]=1;
        }
    }
    for (i=2; i<=n; i++)
    {
        x=i;
        phi=x;
        for (d=1; d<=k&&v[d]*v[d]<=x; d++)
        {
            if (x%v[d]==0)
            {
                phi=phi/v[d]*(v[d]-1);
                while (x%v[d]==0)
                    x/=v[d];
            }
        }
        if (x!=1)
            phi=phi/x*(x-1);
        sum+=phi;
    }
    sum=sum*2+1;
    fout<<sum;
    return 0;
}