// Layout.cpp : Defines the entry point for the console application.
/*
 作者：成曉旭
 時間：2001年10月11日(11:35:38-12:35:00)
 內容：完成插花問題的“動態規劃法”算法及注解
*/
#include "stdafx.h"
#include "string.h"
#define MAX(A,B) ((A) > (B) ? (A):(B))
//--------------------鮮花問題--------------------
#define F 100
#define V 100
/*
 插花問題描述:
   將f束鮮花插入v個花瓶中,使達到最徍的視覺效果,
  問題相關約定及插花要求:
   鮮花被編號為1--f,花瓶被編號為1--v,花瓶按從小到
  大順序排列,一只花瓶只能插一支花,鮮花i插入花瓶j中的
  視覺效果效果值已知,編號小的鮮花所放入的花瓶編號也小 
 問題求解思路:
   花瓶j(1<=j<=v)中插入鮮花的可能編號為[1..j](編號
  小的鮮花所放入的花瓶編號也小);
   設數組p[i][j]表示鮮花i插入花瓶j的好看程度,數組
  q[i][j]表示[1..i]束鮮花插入[1..j]個花瓶所能得到的最大
  好看程度,初始化q[0][0] = 0;q[0][j]=0(1<=j<=v),則q[f][v]
  是問題的解.
   特別地,j束鮮花插入到前面的j只花瓶中,所得到的好看
  程度是q[j][j] = p[1][1]+p[2][2]+...+[j][j].現將插花過
  程按花瓶排列順序劃分成不同階段,則在第j階段,第i束鮮花
  若放入第j號花瓶,最大好看程度是q[i-1][j-1]+p[i][j];第i束鮮
  花若放入前j-1個花瓶中的某一個,所得的好看程度是q[i][j-1],
  那么在第j階段,插入第i束鮮花所能得到的最大好看程度為:
  q[i][j] = MAX(q[i-1][j-1]+p[i][j],q[i][j-1]),要使q[i][j]
  最大,應使q[i-1][j-1]和q[i][j-1]也最大
*/
//初始化函數
void Initialize(int *f,int *v,int p[][V])
{
 int i,j;
 printf("輸入鮮花數量及花瓶個數:");
 scanf("%d%d",f,v);
 printf("順序輸入各鮮花插入各花瓶的好看程度:\n");
 for(i=1;i<=*f;i++)
  for(j=1;j<=*v;j++)
   p[i][j] = i+j;
   //scanf("%d",&p[i][j]);
}
//鮮花問題處理函數
int Sove(int p[][V],int f,int v,int *way)
{
 int i,j,newv,q[F][V];
 q[0][0] = 0;
 /*設置v個花瓶分別被插入v束鮮花時各號花瓶對應的(初始)最大好看程度*/
 for(j=1;j<=v;j++)
 {
  q[0][j] = 0;
  /*設置第j束鮮花放入第j號花瓶中的最大好看程度*/
  q[j][j] = q[j-1][j-1]+p[j][j];
 }
 for(j=1;j<=v;j++)
  for(i=1;i<j;i++)
   q[i][j] = MAX(q[i-1][j-1]+p[i][j],q[i][j-1]);
 newv = v;
 for(i=f;i>0;i--)
 {
  while(q[i-1][newv-1]+p[i][newv] < q[i][newv])
   newv--;
  //確定鮮花i插在花瓶newv中,并準備考慮前一只花瓶　
  way[i] = newv--;
 }
 return(q[f][v]);
}
//--------------------鮮花問題--------------------
//--------------------最長子串問題--------------------
#define N 100
char a[N],b[N],str[N];
//計算兩個序列最長公共子序列的長度
int Get_LongSubStr_Len(char *a,char *b,int c[][N])
{
 int m=strlen(a),n=strlen(b),//兩個序列的長度
  i,j;//循環變量
 for(i=0;i<=m;i++) c[i][0] = 0;
 for(i=1;i<=n;i++) c[0][i] = 0;
 for(i=1;i<=m;i++)
  for(j=1;j<=n;j++)
   if(a[i-1]==b[j-1])
    c[i][j] = c[i-1][j-1]+1;
   else
    c[i][j] = MAX(c[i-1][j],c[i][j-1]);
   /*
    if(c[i-1][j]>=c[i][j-1])
     c[i][j] = c[i-1][j];
    else
     c[i][j] = c[i][j-1];
     */
 return(c[m][n]);
}
//構造最長公共子序列
char *Build_LongSubStr(char s[],char *a,char *b)
{
 int i=strlen(a),j=strlen(b),
  k,c[N][N];
 k = Get_LongSubStr_Len(a,b,c);
 s[k] = '\0';
 while(k>0)
 {
  if(c[i][j]==c[i-1][j])
   i--;
  else
  {
   if(c[i][j]==c[i][j-1])
    j--;
   else
   {
    s[--k]=a[i-1];
    i--;
    j--;
   }
  }
 }
 return(s);
}
//--------------------最長子串問題--------------------
int main(int argc, char* argv[])
{
 int i,f,v,p[F][V],way[F];
 //-----------------------------------
 /*
 Initialize(&f,&v,p);
 printf("最大好看程度為%d\n",Sove(p,f,v,way));
 printf("插有鮮花的花瓶是:\n");
 for(i=1;i<=f;i++)
  printf("%4d",way[i]);
 */
 //-----------------------------------
 printf("輸入兩個字符串(長度<%d):\n",N);
 scanf("%s%s",a,b);
 printf("兩個串的最長公共子序列是:%s\n",Build_LongSubStr(str,a,b));
 //-----------------------------------
 printf("\n\n應用程序正在運行......\n");
 return 0;
}