Cod sursa(job #630047)

Utilizator ChallengeMurtaza Alexandru Challenge Data 4 noiembrie 2011 16:25:06
Problema Indep Scor 100
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.75 kb
#include <fstream>

using namespace std;

const char InFile[]="indep.in";
const char OutFile[]="indep.out";
const int HUGEINT_SIZE=256;
const int HUGEINT_BASE=10;
const int MaxX=1024;
const int MaxN=512;

ifstream fin(InFile);
ofstream fout(OutFile);

typedef int HUGEINT[HUGEINT_SIZE];

int N,x,D[MaxX],P[MaxX],CNT[MaxX];
HUGEINT SOL,POW[MaxN];

inline void add(HUGEINT A, HUGEINT B)
{
	int i=1,t=0;
	for(;i<=A[0] || i<=B[0] || t;++i,t/=HUGEINT_BASE)
	{
		A[i]=(t+=(A[i]+B[i]))%HUGEINT_BASE;
	}
	A[0]=i-1;
}

inline void cpy(HUGEINT A, HUGEINT B)
{
	for(register int i=0;i<=B[0];++i)
	{
		A[i]=B[i];
	}
}

inline void write(HUGEINT A)
{
	if(A[0]==0)
	{
		fout<<"0";
	}
	else
	{
		for(register int i=A[0];i>0;--i)
		{
			fout<<A[i];
		}
	}
}

inline void sub(HUGEINT A, HUGEINT B)
{
	int i=1,t=0;
	for(;i<=A[0];++i)
	{
		A[i]-=B[i]+t;
		A[i]+=(t=A[i]<0)*10;
	}
	while(!A[A[0]] && A[0])
	{
		--A[0];
	}
}

int main()
{
	fin>>N;
	for(register int i=1;i<=N;++i)
	{
		fin>>x;
		for(register int j=2;j<=x;++j)
		{
			if(x%j==0)
			{
				++D[j];
			}
		}
	}
	fin.close();
	
	for(register int i=2;i<MaxX;++i)
	{
		P[i]=1;
	}
	for(register int i=2;i<MaxX;++i)
	{
		if(P[i]==1)
		{
			for(register int j=i;j<MaxX;j+=i)
			{
				P[j]*=i;
				++CNT[j];
			}
		}
	}
	
	POW[0][0]=POW[0][1]=1;
	for(register int i=1;i<=N;++i)
	{
		cpy(POW[i],POW[i-1]);
		add(POW[i],POW[i-1]);
	}
	
	for(register int i=1;i<=N;++i)
	{
		sub(POW[i],POW[0]);
	}
	sub(POW[0],POW[0]);
	
	cpy(SOL,POW[N]);
	for(register int i=2;i<MaxX;++i)
	{
		if(P[i]==i)
		{
			if(!(CNT[i]&1))
			{
				if(D[i])
				{
					add(SOL,POW[D[i]]);
				}
			}
		}
	}
	for(register int i=2;i<MaxX;++i)
	{
		if(P[i]==i)
		{
			if(CNT[i]&1)
			{
				if(D[i])
				{
					sub(SOL,POW[D[i]]);
				}
			}
		}
	}
	
	write(SOL);
	fout.close();
	return 0;
}