Cod sursa(job #734870)

Utilizator desoComan Andrei deso Data 15 aprilie 2012 03:42:57
Problema Perb Scor 10
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.35 kb
#include <iostream>
#include <fstream>
#include <cstring>
#include <cstdlib>
#include <string>
#include <vector>
#include <algorithm>
#include <set>
using namespace std;

#define max(a,b) (a>b ? a : b)
#define min(a,b) (a<b ? a : b)
#define INFILE "perb.in" 
#define OUTFILE "perb.out"
#define MAX 602
ifstream fin (INFILE);
ofstream fout (OUTFILE);

int c[MAX][MAX], n, v[MAX/2][4];
char s[MAX];

int main()
{
   int m;
   fin >> n >> m >> s;

   for(int i=0; i<n; i++)
      if( s[i]=='A' ) s[i] = 0;
      else if( s[i]=='C' ) s[i] = 1;
      else if( s[i]=='G' ) s[i] = 2;
      else if( s[i]=='T' ) s[i] = 3;

   memset(c, 0x3f, sizeof(c));

   for(int i=0; i<n-1; i++)
   for(int p=1; i+2*p<n && p<=n/2; p++)
   {
      for(int j=i, k=0; j<i+p; j++, k++)
      {
         v[k][0] = v[k][1] = v[k][2] = v[k][3] = 0;
         v[k][s[j]]++;
      }

      for(int j=i+p, k=0, ps=2, changes=0; j<n; j++, k++)
      {
         v[k][s[j]]++;
         changes += ps - max(max(v[k][0], v[k][1]), max(v[k][2], v[k][3]));

         if( k==p-1 )
         {
            k = -1;    
            c[i][j] = min(c[i][j], changes);
            changes = 0;
            ps++;
         }
      }
   }

   for(int i=0; i<m; i++)
   {
      int x, y;
      fin >> x >> y;
      if( x>y ) swap(x, y);
      if( x==y )  fout << "0\n";
      else        fout << c[x-1][y-1] << endl;
   }
	
   return 0;
}