Cod sursa(job #894788)

Utilizator veleanduAlex Velea veleandu Data 26 februarie 2013 23:33:54
Problema Fractii Scor 30
Compilator cpp Status done
Runda Arhiva de probleme Marime 0.81 kb
#include <cstdio>
#include <iostream>
using namespace std;
#define maxn 1001000
 
    int V[maxn],Rez[maxn];
    int n,i,j,l;
 
 
int main()
{
    freopen ("fractii.in","r",stdin);
    freopen ("fractii.out","w",stdout);
    scanf ("%d", &n);
 
    for ( i=1; i<=n; i++ )
        Rez[i]=1;
	long long rez=0;
    for ( i=2; i<=n; i++ )
    {
        if ( !V[i] )
        {
			rez--;
            for ( j=i*i; j<=n; j+=i )
                V[j]=1;
			Rez[i]*=i;
            for ( j=2*i; j<=n; j+=i )
                Rez[j]*=(i-1);
            for ( l=i*i; l<=n; l*=i )
                for ( j=l; j<=n; j+=l )
                    Rez[j]*=i;
        }
    }
    if ( !n )
    {
        printf("0");
        return 0;
    }
    for ( i=2; i<=n; i++ )
        rez+=(long long)Rez[i];
	rez<<=1;
    rez++;
    cout<<rez<<"\n";
    return 0;
}