Cod sursa(job #87725)

Utilizator cos_minBondane Cosmin cos_min Data 28 septembrie 2007 20:03:52
Problema Cast Scor 40
Compilator cpp Status done
Runda Arhiva de probleme Marime 3.35 kb
#include <stdio.h>
#include <fstream>
using namespace std;

#define in "cast.in"
#define out "cast.out"
#define dim 13

int N, T, Tmin=-1;
int B[dim][8192], L[dim];
int A[dim][dim];

int Maxim(int,int);
int Ok(int,int);
void Solve();

int main()
{
    freopen(in,"r",stdin);
    freopen(out,"w",stdout);
    
    scanf("%d", &T);
    for ( int k = 1; k <= T; k++ )
    {
        scanf("%d", &N);
        for ( int i = 1; i <= N; i++ )
            for ( int j = 1; j <= N; j++ )
                scanf("%d", &A[i][j]);
                
        Solve();
    }                
}    

int Maxim(int a, int b)
{
    if ( a > b ) return a;
    return b;
}    

int Minim(int a, int b)
{
    if ( a < b ) return a;
    return b;
    //return a < b ? a : b;
}

int Ok(int k1, int k2)
{
    for ( int i = 1; i <= N; i++ )
        if ( (k2>>(i-1))&1 ) 
             if ( !((k1>>(i-1))&1) ) return 0;
    
    return 1;  
}    

void Solve()
{
    int S = 1<<N;
    S -= 1;
    
    memset(B, 0x3f, sizeof(B) );
    for ( int i = 1; i <= N; i++ )
        B[i][ (1<<(i-1)) ] = 0;
   
   /* for ( int k1 = 1; k1 <= S; ++k1 )
    {
        for ( int i = 1; i <= N && 1<<(i-1) <= k1; ++i )
        {
            if ( (k1>>(i-1))&1 == 0 ) continue;
            for ( int k2 = 1; k2 < k1; ++k2 )
            {
                if ( !Ok(k1,k2) ) continue;
                for ( int j = 1; j <= N && 1<<(j-1) <= k2; ++j )
                {
                    if ( (k2>>(j-1))&1 == 0 || i == j ) continue;
                    
                    B[i][k1] = Minim( B[i][k1], A[i][j] + Maxim(B[i][k1-k2],B[j][k2]) ); 
                }             
            }
        }
    }*/
    
    int size, Q, k, i, k1, k2, j, t;
    
    int pow=0;
    
    for ( k1 = 1; k1 <= S; ++k1 )
    {
          for ( i = 1; i <= N; ++i )
          {
              if ( (k1>>(i-1))&1 )
              {
                   if ( k1 == pow ) continue;
                   
                   size = 0;
                   pow = 1;
                   
                   for ( j = 1; pow <= k1 && j <= N; ++j )
                   {
                       pow = (1<<(j-1));
                       if ( i != j && (k1>>(j-1))&1 ) L[++size] = j;
                   }
                   
                   pow = (1<<size);
                   
                   for ( k2 = 1; k2 < pow; ++k2 )
                   {
                       Q = 0;
                       for ( t = 1; t <= size; ++t )
                           if ( (k2>>(t-1))&1 ) Q += (1<<(L[t]-1));
                       
                       if ( k1 < Q ) continue;
                       
                       int R = k1-Q;
                       
                       for ( t = 1; t <= size; ++t )
                       {
                           if ( B[i][k1] > A[i][L[t]] + Maxim(B[i][R],B[L[t]][Q])  ) B[i][k1] = A[i][L[t]] + Maxim(B[i][R],B[L[t]][Q]);
                          // B[i][k1] = Minim( B[i][k1], A[i][L[t]] + Maxim(B[i][R],B[L[t]][Q]) );
                       }                   
                   }
              }
          }
    } 
    
   /* printf("\n");
     for ( int i = 1; i <= N; i++, printf("\n") )
        for ( int j = 1; j <= (1 << N) - 1; j++ )
            printf("%d ", B[i][j]);*/
    
    printf("%d\n", B[1][S]);
}