6

是否有一种算法可以通过给定的最大允许位置数量（最大不匹配，最大汉明距离）生成字符串（DNA序列）的所有可能字符串组合？

字母表是 {A,C,T,G}。

字符串AGCC和最大数量的 Mismatches示例2：

Hamming distance is 0
  {AGCC}
Hamming distance is 1
  {CGCC, TGCC, GGCC, AACC, ACCC, ATCC, AGAC, AGTC, ..., AGCG}
Hamming distance is 2
  {?}

一种可能的方法是生成一个具有给定字符串的所有排列的集合，迭代它们并删除所有具有更大汉明距离的字符串。

这种方法非常消耗资源，给定的 20 个字符的字符串和 5 的最大汉明距离。

是否有另一种更有效的方法/实现？

4

1 回答 1

8

只需使用正常的排列生成算法，除了你绕过距离，当你有不同的字符时减少它。

static void permute(char[] arr, int pos, int distance, char[] candidates)
{
   if (pos == arr.length)
   {
      System.out.println(new String(arr));
      return;
   }
   // distance > 0 means we can change the current character,
   //   so go through the candidates
   if (distance > 0)
   {
      char temp = arr[pos];
      for (int i = 0; i < candidates.length; i++)
      {
         arr[pos] = candidates[i];
         int distanceOffset = 0;
         // different character, thus decrement distance
         if (temp != arr[pos])
            distanceOffset = -1;
         permute(arr, pos+1, distance + distanceOffset, candidates);
      }
      arr[pos] = temp;
   }
   // otherwise just stick to the same character
   else
      permute(arr, pos+1, distance, candidates);
}

致电：

permute("AGCC".toCharArray(), 0, 1, "ACTG".toCharArray());

性能说明：

对于 20 的字符串长度、5 的距离和 5 个字符的字母表，已经有超过 1700 万个候选者（假设我的代码是正确的）。

上面的代码在我的机器上花费不到一秒钟的时间（不打印），但不要指望任何生成器能够在合理的时间内生成更多的东西，因为有太多的可能性.

于 2013-10-09T10:46:15.713 回答