LightOJ 1110 An Easy LCS--DP

来源:互联网 发布:delphi源码 编辑:程序博客网 时间:2024/05/17 08:07
/*找出最长的字典序最小的公共子串这是我同学想的一方法在动态规划的时候给每个状态都求出   最长的字典序最小的公共子串动态规划是二维的,同时还要存一个公共子串 所以需要一个三维的字符数组*/#include<stdio.h>#include<string.h>int dp[110][110],ls,lss;//dp动态规划用,记录长度  ls、lss分别代表s、ss的长度char s[110],ss[110],map[110][110][110];//map 就是记录各个状态的公共子串的三维数组void pro(){int i,j,ii,jj;//ii、jj指示  最长的字典序最小的  那个状态j=ls;if(lss>j)j=lss;for(i=0;i<=j;++i)//初始化{dp[0][i]=dp[i][0]=0;map[i][0][0]=map[0][i][0]=0;}for(i=1;i<=lss;++i)//开始dp{for(j=1;j<=ls;++j){if(ss[i]==s[j])//可以增长   dp[i-1][j-1]以后不会用了  所以可以修改成增长以后的情况  方便一会儿复制{dp[i-1][j-1]++;map[i-1][j-1][dp[i-1][j-1]-1]=ss[i];map[i-1][j-1][dp[i-1][j-1]]=0;}ii=i-1;//假设dp[i-1][j-1]是  最长的字典序最小的jj=j-1;if(dp[i][j-1]>dp[ii][jj]||(dp[i][j-1]==dp[ii][jj]&&strcmp(map[i][j-1],map[ii][jj])<0))//若dp[i][j-1]是  最长的字典序最小的  则更新{ii=i;jj=j-1;}if(dp[i-1][j]>dp[ii][jj]||(dp[i-1][j]==dp[ii][jj]&&strcmp(map[i-1][j],map[ii][jj])<0))//若dp[i-1][j]是  最长的字典序最小的  则更新{ii=i-1;jj=j;}dp[i][j]=dp[ii][jj];//把  最长的字典序最小的  复制到当前strcpy(map[i][j],map[ii][jj]);}}}int main()//主函数就不用解释了吧{int t,i;ss[0]=s[0]=' ';scanf("%d",&t);for(i=1;i<=t;++i){scanf("%s%s",s+1,ss+1);ls=strlen(s+1);lss=strlen(ss+1);pro();if(dp[lss][ls]==0)printf("Case %d: :(\n",i);else printf("Case %d: %s\n",i,map[lss][ls]);}return 0;}

原创粉丝点击