Cod sursa(job #894793)

Utilizator veleanduAlex Velea veleandu Data 26 februarie 2013 23:35:28
Problema Fractii Scor 30
Compilator cpp Status done
Runda Arhiva de probleme Marime 0.82 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] )
        {
            for ( j=1LL*i*i; 1LL*j<=1LL*n; j+=i )
                V[j]=1;
            for ( j=i; 1LL*j<=1LL*n; j+=i )
                Rez[j]*=(i-1);
            for ( l=i*i; 1LL*l<=1LL*n; l*=i )
                for ( j=l; 1LL*j<=1LL*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;
}