Cod sursa(job #87749)

Utilizator cos_minBondane Cosmin cos_min Data 28 septembrie 2007 21:31:47
Problema Cast Scor 40
Compilator cpp Status done
Runda Arhiva de probleme Marime 4.55 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][8193], L[dim];
bool Sel[8192][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);
    
    memset(Sel,0,sizeof(Sel));
    
    for ( int k1 = 1; k1 <= 8192; k1++ )
    {
        for ( int i = 1; i <= 13; i++ )
            if ( (k1>>(i-1)) & 1 ) Sel[k1][i] = 1;
    }
    
    scanf("%d", &T);
    for ( ; T >= 1; --T )
    {
        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) );
    //memset(Sel, 0, sizeof(Sel) );
    
    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]) );
                       }                   
                   }
              }
          }
    } 
    */
    
    int size=0, Q=0;
    
    for ( int k1 = 1; k1 <= S; k1++ )
    {
        for ( int i = 1; i <= N; i++ )
        {
            if ( Sel[k1][i] )
            {
                 size = 0;
                 
                 for ( int j = 1; j <= N; j++ )
                 {
                     if ( j != i && Sel[k1][j] ) L[++size] = j;  
                 }
                 
                 for ( int k2 = 1; k2 < (1<<size); k2++ )
                 {
                     Q = 0;
                     
                     for ( int t = 1; t <= size; t++ )
                     {
                         if ( Sel[k2][t] ) Q += 1<<(L[t]-1);
                     }
                     
                     for ( int t = 1; t <= size; t++ )
                         B[i][k1] = Minim( B[i][k1], A[i][L[t]] + Maxim(B[i][k1-Q],B[L[t]][Q]) ); 
                 }
            }
        }
    }
           // size = 0;
            
   
    
   /* 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]);
}