Pagini recente » Cod sursa (job #1216811) | Cod sursa (job #844672) | Cod sursa (job #3260481) | Cod sursa (job #2340112) | Cod sursa (job #2161806)
#include <iostream>
#include <fstream>
#include <math.h>
#define NMAX 505
using namespace std;
ifstream f("indep.in");
ofstream g("indep.out");
int p2[NMAX][NMAX],ans[NMAX],prim[NMAX],rez[NMAX];
int n,a[NMAX],MAX,nr,semn;
bool ciur[NMAX];
void prelucrare()
{
p2[0][0]=p2[0][1]=1;
for(int i=1; i<=MAX; i++)
{
for(int j=0; j<=p2[i-1][0]; j++)
p2[i][j]=p2[i-1][j];
int t=0,p;
for(int j=1; j<=p2[i][0]; j++)
{
p=p2[i][j]*2+t;
t=p/10;
p2[i][j]=p%10;
}
if(t!=0)
{
p2[i][0]++;
p2[i][p2[i][0]]=t;
}
}
prim[1]=2;
nr=1;
for(int i=4; i<=MAX; i+=2)
ciur[i]=1;
for(int i=3; i<=MAX; i+=2)
if(ciur[i]==0)
{
prim[++nr]=i;
for(int j=2*i; j<=MAX; j+=i)
ciur[j]=1;
}
}
void read()
{
f>>n;
for(int i=1; i<=n; i++)
f>>a[i],MAX=max(MAX,a[i]);
f.close();
}
int okay;
int verif(int x)
{
int sol=0; okay=0;
for(int i=1; i<=nr&&x>1; i++)
{
int ok=0;
while(x%prim[i]==0&&x>1)
{
ok++;
x/=prim[i];
}
if(ok>1)
okay=1;
else
if(ok!=0)
sol++;
}
return sol;
}
int divv(int x)
{
int sol=0;
for(int i=1; i<=n; i++)
if(a[i]%x==0)
sol++;
return sol;
}
void adunare(int a[],int b[])
{
int t=0,i;
for(i=1; i<=max(a[0],b[0])||t; i++)
{
a[i]=(t+a[i]+b[i])%10;
t/=10;
}
a[0]=i-1;
}
void scadere(int a[],int b[])
{
int t=0;
for(int i=1; i<=a[0]; i++)
{
a[i]=a[i]-b[i]-t;
if(a[i]>=0)
t=0;
else
{
a[i]+=10;
t=1;
}
}
while(a[0]>1&&a[a[0]]==0)
a[0]--;
}
void write();
void solve()
{
int x=n;
for(int i=2; i<=MAX; i++)
{
int sol=verif(i),s=0;
if(okay==0)s=divv(i);
if(sol%2==1)
adunare(ans,p2[s]);
else
scadere(ans,p2[s]);
}
adunare(rez,p2[n]);
scadere(rez,ans);
}
void write()
{
for(int i=ans[0]; i>=1; i--)
g<<rez[i];
g.close();
}
int main()
{
read();
prelucrare();
solve();
write();
return 0;
}