Cod sursa(job #1149388)

Utilizator PatrikStepan Patrik Patrik Data 21 martie 2014 19:55:16
Problema Indep Scor 15
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.88 kb
   #include<cstdio>
    #include<cstring>
    using namespace std;
    #define NMAX 501
    #define VMAX 1001
    #define BASE 1000000000
    int N , v[NMAX] , g[VMAX][VMAX];
    struct huge{
        int nr[200];
    }D[2][VMAX] , unu;

    void add(huge &A , huge B)
    {
        int t = 0 , p = 1;
        for( p = 1 ; p <= A.nr[0] || p <= B.nr[0] || t ; ++p , t/=BASE)
        {
            t += A.nr[p] + B.nr[p];
            A.nr[p] = t%BASE;
        }
        A.nr[0] = p-1;
    }

    void cpy(huge &A , huge B)
    {
        for(int i = 1 ; i <= B.nr[0] ; ++i )
            A.nr[i] = B.nr[i];
        A.nr[0] = B.nr[0];
    }


    void seteaza(huge &A , int val)
    {
        A.nr[1] = val;
        A.nr[0] = 1;
    }

    int gcd(int a , int b)
    {
       while(b)
       {
           int r = a%b;
           a = b;
           b = r;
       };
       return a;
    }

    void tipar(huge A)
    {
        printf("%d" , A.nr[A.nr[0]]);
        for(int i = A.nr[0]-1 ; i ; i-- )
            printf("%09d" , A.nr[i]);
    }

    int main()
    {
        freopen("indep.in" , "r" , stdin );
        freopen("indep.out" , "w" , stdout );
        scanf("%d" , &N );
        for(int i = 1 ; i <= N ; ++i )
            scanf("%d" , &v[i]);
        for(int i = 1 ; i < VMAX ; ++i )
            for(int j = 1 ; j < VMAX ; ++j )
            g[i][j]  = g[j][i] = gcd(i,j);

        seteaza (D[0][v[1]] , 1);
        unu.nr[1] = unu.nr[0] = 1;
        for(int i = 2 , j = 1; i <= N ; ++i , j = 1-j )
        {
            for(int k = 1 ; k < VMAX ; ++k )
                    cpy(D[j][k],D[1-j][k]);
            for(int k = 1 ; k < VMAX ; ++k )
               add( D[j][g[k][v[i]]] , D[1-j][k]);
            add(D[j][v[i]],unu);
        }
        if(N%2)
            tipar(D[0][1]);
        else
        tipar(D[1][1]);
        return 0;
    }