USACO5.4-Character Recognition

题目大意是字符串识别
一道细节很繁琐的DP,要用到很多数组
一开始还真看不出是DP,后来参考了别人的代码,然后又按自己的思路重头到尾写了,虽然速度不咋的

Executing...
Test 1: TEST OK [0.008 secs, 6504 KB]
Test 2: TEST OK [0.008 secs, 6504 KB]
Test 3: TEST OK [0.019 secs, 6504 KB]
Test 4: TEST OK [0.035 secs, 6504 KB]
Test 5: TEST OK [0.084 secs, 6504 KB]
Test 6: TEST OK [0.116 secs, 6504 KB]
Test 7: TEST OK [0.167 secs, 6504 KB]
Test 8: TEST OK [0.192 secs, 6504 KB]
Test 9: TEST OK [0.178 secs, 6504 KB]

All tests OK.
YOUR PROGRAM ('charrec') WORKED FIRST TIME! That's fantastic
-- and a rare thing. Please accept these special automated
congratulations.

 #include <stack>
 #include <iostream>
 #include <stdio.h>
 #define R19 0
 #define R20 1
 #define R21 2
 #define INF 1000000000
 using namespace std;

 char ch[][][];
 char mat[][];
 int dif[][][];
 int cost[][]; // cost[i][j]表示从mat的第i行开始匹配j行得到的最小差距
 int optimalCh[][]; // optimalCh[i][j]表示从mat的第i行开始匹配j行最优的匹配字符
 int d[]={}; // d[i]表示匹配到第i行所得到的最小差距
 int step[]; // step[i]表示第i行开始的最优匹配的字符行数
 // 得到行数中不同的个数
 /*
 num为第几个字符
 l1为num字符的行数
 l2为mat的行数
 */
 int getDif(int num,int l1,int l2)
 {
     int sum=;
     for(int i=;i<;i++)
     {
         if(ch[num][l1][i]!=mat[l2][i])
             sum++;
     }
     return sum;
 }

 // num为字符，l1为mat开始匹配的行数，del为删除了哪一行
 int getCost(int num,int l1,int del)
 {
     int sum=;
     for(int i=,p=;p<;i++,p++)
     {
         if(i==del)
             i++;
         sum+=dif[num][i][l1+p];
     }
     return sum;
 }

 int getCost2(int num,int l1,int del)
 {
     int sum=;
     for(int i=,p=;i<;i++,p++)
     {
         if(p==del)
             p++;
         sum+=dif[num][i][p+l1];
     }
     return sum;
 }

 int main()
 {
     freopen("font.in","r",stdin);
     int n; // 行数
     cin>>n;
     for(int i=;i<n/;i++)
     {
         for(int j=;j<;j++)
         {
             scanf("%s",ch[i][j]);

         }
     }

     freopen("charrec.in","r",stdin);
     freopen("charrec.out","w",stdout);

     cin>>n;

     for(int i=;i<n;i++)
         scanf("%s",mat[i]);

     // 初始化dif数组
     for(int i=;i<;i++)
         for(int j=;j<;j++)
         {
             for(int k=;k<n;k++)
             {
                 dif[i][j][k]=getDif(i,j,k);
             }
         }

     // 初始化cost
     for(int i=;i<n;i++) // 从第i行开始匹配
     {
         cost[i][R19]=cost[i][R20]=cost[i][R21]=INF;
         for(int j=;j<;j++) // 第j个字符
         {
             if(i+<n+)
             {
                 for(int k=;k<;k++) // 删掉行数
                 {
                     if(cost[i][R19]>getCost(j,i,k))
                     {
                         cost[i][R19]=getCost(j,i,k);
                         optimalCh[i][R19]=j;
                     }
                 }
             }
             if(i+<n+)
             {
                 if(cost[i][R20]>getCost(j,i,-))
                 {
                     cost[i][R20]=getCost(j,i,-);
                     optimalCh[i][R20]=j;
                 }
             }
             if(i+<n+)
             {
                 for(int k=;k<;k++) // 多出的行数
                 {
                     if(cost[i][R21]>getCost2(j,i,k))
                     {
                         cost[i][R21]=getCost2(j,i,k);
                         optimalCh[i][R21]=j;
                     }
                 }
             }
         }
     }

     fill_n(d,n+,INF);
     // 初始值
     d[]=cost[][R19];step[]=;
     d[]=cost[][R20];step[]=;
     d[]=cost[][R21];step[]=;
     // DP
     for(int i=;i<n;i++) // 匹配的最后一行为i
     {
         if(d[i-]+cost[i-][R19]<d[i])
         {
             d[i]=d[i-]+cost[i-][R19];
             step[i]=;
         }
         if(d[i-]+cost[i-][R20]<d[i])
         {
             d[i]=d[i-]+cost[i-][R20];
             step[i]=;
         }
         if(d[i-]+cost[i-][R21]<d[i])
         {
             d[i]=d[i-]+cost[i-][R21];
             step[i]=;
         }
     }

     stack<int> S;

     for(int p=n-;p!=-;p=p-step[p])
     {
         S.push(optimalCh[p-step[p]+][step[p]-]);
     }

     while(!S.empty())
     {
         int tmp=S.top();
         S.pop();
         printf("%c",tmp==?' ':tmp-+'a');
     }
     cout<<endl;
     return ;
 }