hihocoder#1052之基因工程



描述

小Hi和小Ho正在进行一项基因工程实验。他们要修改一段长度为N的DNA序列，使得这段DNA上最前面的K个碱基组成的序列

与最后面的K个碱基组成的序列完全一致。  

例如对于序列"ATCGATAC"和K=2,可以通过将第二个碱基修改为"C"使得最前面2个碱基与最后面两个碱基都为"AC"。

当然还存在其他修改方法，例如将最后一个碱基改为"T"，或者直接将最前面两个和最后面两个碱基都修改为"GG"。

小Hi和小Ho希望知道在所有方法中，修改碱基最少的方法需要修改多少个碱基。

输入

第一行包含一个整数T(1 <= T <= 10)，代表测试数据的数量。

每组测试数据包含2行，第一行是一个由"ATCG"4个大写字母组成的长度为N(1 <= N <= 1000)的字符串。

第二行是一个整数K(1 <= K <= N)。

输出

对于每组数据输出最少需要修改的碱基数量。

样例输入 2 
ATCGATAC 
2 
ATACGTCT
6 样例输出 1 
3

#include 
    
    
     
     
#include
     
     
      
      
#include
      
      
       
       

using namespace std;

int main()
{
    int i, j, l, t, n, k, m, mn, cnt, num[4];
    char s[1111];
    cin>>t;
    while(t--)
    {
        m = 0;
        memset(s, 0, sizeof(s));
        cin>>s>>k;
        n = strlen(s);
        if(k <= n / 2)
        {
            for(i = 0; i < k; i++)
            {
                if(s[i] != s[n - k + i])
                    m++;
            }
        }
        else
        {
            for(i = 0; i < n - k; i++)
            {
                mn = INT_MAX;
                memset(num, 0, sizeof(num));
                for(j = i, cnt = 0; j < n; j += n - k, cnt++)
                {
                    if(s[j] == 'A')
                        num[0]++;
                    else if(s[j] == 'C')
                        num[1]++;
                    else if(s[j] == 'G')
                        num[2]++;
                    else
                        num[3]++;
                }
                for(l = 0; l < 4; l++)
                {
                    mn = mn < cnt - num[l] ? mn : cnt - num[l];
                }
                m += mn;
            }
        }
        cout<
       
       
        
        <