描述
小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<
<