Cod sursa(job #533356)

Utilizator SadmannCornigeanu Calin Sadmann Data 13 februarie 2011 19:53:17
Problema Fractii Scor 10
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.32 kb
#include<stdio.h>
#include<vector>
#include<math.h>
using namespace std;
int cmmdc(int a,int b)
{
    while(b)
    {
        int r=a%b;
        a=b;
        b=r;
    }
    return a;
}


FILE *in,*out;
int N,cntprim,i,j,nr,contscad;
double euler;
bool primviz[1001];  //cele cu false sunt prime
double totprim[1001];
double scadere[1001];
int main()
{
    in=fopen("fractii.in","rt");
    out=fopen("fractii.out","wt");
    fscanf(in,"%d",&N);
	euler=1;
	//prim[1]=2;
	//cntprim=1;
	for(i=2;i<=N;i++)
	{
		if(!primviz[i])
		{
			euler+=2*(i-1);
			for(j=i<<1;j<=N;j+=i)
			{
				primviz[j]=true;
				if(scadere[j]>1)
					nr=j/scadere[j];
				else
				{
					nr=j;
					scadere[j]=1;
				}
				contscad=0;
				if(!totprim[j])
					totprim[j]=j;
				while( !(nr%i) && nr>1 )
				{
					nr/=i;
					contscad++;
				}
				if(nr<=1)
				{
					euler+=2*(totprim[j]-totprim[j]/i);
					totprim[j]=0;
				}
				else
				{
					totprim[j]=totprim[j]-totprim[j]/i;
					scadere[j]*=pow(i,contscad);
				}
			}
		}
	}
   /* for(i=2;i<=N;i++)
    {
		nr=i;
		double prod=i;
		for(j=1;j<=cntprim && prim[j]<=nr;j++)
		{
			if( !(nr%prim[j]) )
			{
				nr/=prim[j];
				prod*=(1-(double)1/prim[j]);
			}
		}
		euler+=2*prod;
    }*/
	fprintf(out,"%.0lf",euler);
	
    return 0;
}