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.

  1 #include <stack>
  2 #include <iostream>
  3 #include <stdio.h>
  4 #define R19 0
  5 #define R20 1
  6 #define R21 2
  7 #define INF 1000000000
  8 using namespace std;
  9
 10 char ch[28][22][22];
 11 char mat[1210][22];
 12 int dif[28][22][1210];
 13 int cost[1210][4]; // cost[i][j]表示从mat的第i行开始匹配j行得到的最小差距
 14 int optimalCh[1210][4]; // optimalCh[i][j]表示从mat的第i行开始匹配j行最优的匹配字符
 15 int d[1210]={0}; // d[i]表示匹配到第i行所得到的最小差距
 16 int step[1210]; // step[i]表示第i行开始的最优匹配的字符行数
 17 // 得到行数中不同的个数
 18 /*
 19 num为第几个字符
 20 l1为num字符的行数
 21 l2为mat的行数
 22 */
 23 int getDif(int num,int l1,int l2)
 24 {
 25     int sum=0;
 26     for(int i=0;i<20;i++)
 27     {
 28         if(ch[num][l1][i]!=mat[l2][i])
 29             sum++;
 30     }
 31     return sum;
 32 }
 33
 34 // num为字符，l1为mat开始匹配的行数，del为删除了哪一行
 35 int getCost(int num,int l1,int del)
 36 {
 37     int sum=0;
 38     for(int i=0,p=0;p<19;i++,p++)
 39     {
 40         if(i==del)
 41             i++;
 42         sum+=dif[num][i][l1+p];
 43     }
 44     return sum;
 45 }
 46
 47 int getCost2(int num,int l1,int del)
 48 {
 49     int sum=0;
 50     for(int i=0,p=0;i<20;i++,p++)
 51     {
 52         if(p==del)
 53             p++;
 54         sum+=dif[num][i][p+l1];
 55     }
 56     return sum;
 57 }
 58
 59
 60 int main()
 61 {
 62     freopen("font.in","r",stdin);
 63     int n; // 行数
 64     cin>>n;
 65     for(int i=0;i<n/20;i++)
 66     {
 67         for(int j=0;j<20;j++)
 68         {
 69             scanf("%s",ch[i][j]);
 70
 71         }
 72     }
 73
 74     freopen("charrec.in","r",stdin);
 75     freopen("charrec.out","w",stdout);
 76
 77     cin>>n;
 78
 79     for(int i=0;i<n;i++)
 80         scanf("%s",mat[i]);
 81
 82     // 初始化dif数组
 83     for(int i=0;i<27;i++)
 84         for(int j=0;j<20;j++)
 85         {
 86             for(int k=0;k<n;k++)
 87             {
 88                 dif[i][j][k]=getDif(i,j,k);
 89             }
 90         }
 91
 92     // 初始化cost
 93     for(int i=0;i<n;i++) // 从第i行开始匹配
 94     {
 95         cost[i][R19]=cost[i][R20]=cost[i][R21]=INF;
 96         for(int j=0;j<27;j++) // 第j个字符
 97         {
 98             if(i+19<n+1)
 99             {
100                 for(int k=0;k<20;k++) // 删掉行数
101                 {
102                     if(cost[i][R19]>getCost(j,i,k))
103                     {
104                         cost[i][R19]=getCost(j,i,k);
105                         optimalCh[i][R19]=j;
106                     }
107                 }
108             }
109             if(i+20<n+1)
110             {
111                 if(cost[i][R20]>getCost(j,i,-1))
112                 {
113                     cost[i][R20]=getCost(j,i,-1);
114                     optimalCh[i][R20]=j;
115                 }
116             }
117             if(i+21<n+1)
118             {
119                 for(int k=0;k<21;k++) // 多出的行数
120                 {
121                     if(cost[i][R21]>getCost2(j,i,k))
122                     {
123                         cost[i][R21]=getCost2(j,i,k);
124                         optimalCh[i][R21]=j;
125                     }
126                 }
127             }
128         }
129     }
130
131     fill_n(d,n+5,INF);
132     // 初始值
133     d[18]=cost[0][R19];step[18]=19;
134     d[19]=cost[0][R20];step[19]=20;
135     d[20]=cost[0][R21];step[20]=21;
136     // DP
137     for(int i=21;i<n;i++) // 匹配的最后一行为i
138     {
139         if(d[i-19]+cost[i-18][R19]<d[i])
140         {
141             d[i]=d[i-19]+cost[i-18][R19];
142             step[i]=19;
143         }
144         if(d[i-20]+cost[i-19][R20]<d[i])
145         {
146             d[i]=d[i-20]+cost[i-19][R20];
147             step[i]=20;
148         }
149         if(d[i-21]+cost[i-20][R21]<d[i])
150         {
151             d[i]=d[i-21]+cost[i-20][R21];
152             step[i]=21;
153         }
154     }
155
156     stack<int> S;
157
158     for(int p=n-1;p!=-1;p=p-step[p])
159     {
160         S.push(optimalCh[p-step[p]+1][step[p]-19]);
161     }
162
163     while(!S.empty())
164     {
165         int tmp=S.top();
166         S.pop();
167         printf("%c",tmp==0?‘ ‘:tmp-1+‘a‘);
168     }
169     cout<<endl;
170     return 0;
171 }