Cod sursa(job #2114531)

Utilizator andreinichitaTirziu Nichita andreinichita Data 25 ianuarie 2018 17:11:32
Problema Fractii Scor 10
Compilator cpp Status done
Runda ichb_10 Marime 1.19 kb
#include <cstdio>

using namespace std;
const int NMAX=500000;
int c[NMAX+5],prim[100000],phi[10010];
int putere(int a,int b)
{
    int prod=1,i;
    for(i=1; i<=b-1; i++)
        prod*=a;
    return prod*(a-1);
}
int main()
{
  freopen("fractii.in","r",stdin);
  freopen("fractii.out","w",stdout);
    int n,i,fi,j,aux,d,k=0,put;
    long long sum=0;
    c[1]=1;
    for(i=2; i<=NMAX; i++)
        if(c[i]==0){
            prim[++k]=i;
            for(j=2*i; j<=NMAX; j+=i)
                c[j]=1;
        }
    scanf("%d",&n);
    for(i=2; i<=n; i++)
    {
        aux=i;
        fi=1;
        for(d=1; d<=k; d++)
        {
            if(prim[d]>i)
                break;
            put=0;
            while(aux%prim[d]==0)
            {
                put++;
                aux/=prim[d];
            }
            if(put==1)
                fi*=prim[d]-1;
            if(put>=2)
                fi*=putere(prim[d],put);
            if(aux<=1005&&i>=1006)
            {
              fi*=phi[aux];
              break;
            }
        }
        sum+=2*fi;
        if(i<=1005)
          phi[i]=fi;
    }
    printf("%lld",sum+1);
    return 0;
}