Cod sursa(job #2161868)

Utilizator Daria09Florea Daria Daria09 Data 11 martie 2018 21:32:23
Problema Indep Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 2.53 kb
#include <iostream>
#include <fstream>
#include <math.h>
#define CMAX 100005
#define NMAX 505
#define CIURMAX 1005
using namespace std;
ifstream f("indep.in");
ofstream g("indep.out");
int p2[NMAX][159],ans[CMAX],prim[CIURMAX],rez[CMAX];
int n,a[NMAX],MAX,nr,semn,v[NMAX];
bool ciur[CIURMAX];
void read()
{
    f>>n;
    for(int i=1; i<=n; i++)
        f>>a[i],MAX=max(MAX,a[i]);
    f.close();
}
int verif(int x)
{
    int sol=0,okay=0;
    for(int i=1; i<=nr&&1<x; i++)
    {
        int ok=0;
        while(x%prim[i]==0)
        {
            ok++;
            x/=prim[i];
        }
        if(ok>1)
            return 0;
        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<=a[0]||i<=b[0]||t; i++,t/=10)
        a[i]=(t+=a[i]+b[i])%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 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;
        }
    }
    v[0]=v[1]=1;
    for(int i=0;i<=MAX;i++)
        scadere(p2[i],v);
    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 solve()
{
    int x=n;
    for(int i=2; i<=MAX; i++)
    {
        int ok=verif(i);
        if(ok==0)continue;
        int s=divv(i);
        if(s!=0)
        {
            if(ok%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;
}