动态规划1-----------poj1080

 #include<cstdio>
 #include<cstdlib>
 #include<iostream>
 #include<algorithm>

 using namespace std;
 /*dp，poj1080*/

 int dp[][];//动态规划数据存放
 int map[][];//用来存放原始数据

 void map_init()
 {
     map['A']['A']=map['C']['C']=map['G']['G']=map['T']['T']=;
     map['A']['C']=map['C']['A']=map['A']['T']=map['T']['A']=map['T'][' ']=map[' ']['T']=-;
     map['A']['G']=map['G']['A']=map['C']['T']=map['T']['C']=map['G']['T']=map['T']['G']=map['G'][' ']=map[' ']['G']=-;
     map['A'][' ']=map[' ']['A']=map['G']['C']=map['C']['G']=-;
     map['C'][' ']=map[' ']['C']=-;
 }

 int max_X3(int a,int b,int c)
 {
     if(a>b)
     {
         if(a>c)
             return a;
         else
             return c;
     }
     else
     {
         if(b>c)
             return b;
         else
             return c;
     }
 }

 int    main()
 {
     int y;//全局次数
     int i,j;//循环变量
     int a,b;//用户输入
     char str1[];
     char str2[];

     //初始化
     map_init();

     cin>>y;
     while (y--)
     {
         scanf("%d %s",&a,str1);
         scanf("%d %s",&b,str2);

         //初始化第一行第一列
         dp[][]=;
         for (i = ; i < a; i++)
             dp[][i+] = dp[][i] + map[str1[i]][' '];

         for (j = ; j < b; j++)
             dp[j+][] = dp[j][] + map[str2[j]][' '];

         for (i = ; i <= a; i++)
         {
             for (j = ; j <= b; j++)
             {
                 dp[j][i] = max_X3(dp[j-][i-]+map[str2[j-]][str1[i-]],
                     dp[j-][i]+map[str2[j-]][' '],
                     dp[j][i-]+map[str1[i-]][' ']);
             }
         }

         cout<<dp[b][a]<<endl;
     }
     return ;
 }

先上代码，然后说明。

首先对于动态规划到现在的理解，只是现阶段的理解。动态规划下面用dp代替。

要点：

1、一个问题可以被分成多个相同的子问题，子问题和原问题差别只有数据规模，总结就是大化小。

2、这个问题的解可以由子问题的解得出，总结就是用已经有的小解得出最后的解。

3、经常用一个二维数组去保存已经求出的解，之后要是用到就可以直接取不用计算，总结就是记录已经解过的方程的解，不做相同的无用功。

4、状态转移方程，其中的状态，初始值，这些都要考虑清楚。

对于模型，等我做过10道题目熟练之后再说吧，现在给出还为时过早。

题目：

首先明确状态，问问自己两个字符串对比有几种状态？下面用AB两个字符表示两个字符串

1、A的这个字母和B相同

2、A的这个字母和B不同，A和-对应

3、A的这个字母和B不同，B和-对应

最后出来的两个字符串一定是满足这个规则的。只有这三种状态。

那么状态转移方程呢？

1、明确目的，我们要求的是匹配的最后数字上面越大越好。

2、两个字母相同，那么是+5，其他都要减一个数，那么问题来了，这里千万不要以为相同就是最好的，如果这里认为相同就是最好的，那么你用的是贪心而不是dp了。

3、我们取之前的解为Y1，Y2，Y3，123对应上面三个状态,max（Y1+5，Y2-5，Y3-5）这里5是个虚数，可能是0也可能是负数。

4、要注意的有两点，第一，这个状态的之前那个状态的解，是不同的！如下面所示，如果我们处理第二个字符时候可能出现之前的情况就有下面三种

ATGC

GGGG

ATGC

-GGGG

-ATGC

GGGG

第二，Y1，Y2，Y3，是不同的，所以要取这三个值得max。

那么初始值呢？

初始值一般是二维数组的0行0列，这里要注意的是，这里的初始值不全是0，看看下面的情况

ATGC - - - -

- - - -GGGG

所以初始值应该是每一个字母都和空对应，而且要在之前那个值加上去，之前A和空对应如果是-3，那么T和空对应如果是-1，在二维数组中T这里的值应该是-4

最后输出二维数组中最右下角的值就是最后的解。

这里灵活之处是利用一个二维数组去保存了两个字母之间的关系，所以看起来代码特别清晰。

之后会对动态规划再细说的，这次就说这么多。