您的位置:首页 > 其它

7月31日,实验室可以上网。开始写BLOG

2006-07-31 10:57 309 查看
很久没有写BLOG了,原因很简单,就是因为不能上网。整整六天网络极度不稳定。想想我好像已经4天没有AC了。不过最近几天我一直在看算法书。是沙特出的。 发现其中的伪代码写得的确很好。关于动态规划那章我至少看了3个版本。沙特那个版本我看了两遍。写了一个求最大公共子串的函数LCS。没得上网效率不怎么样,看来AC真的能够给人很大的动力。觉得没过多久又要交经验总结了。

一个LCS函数,为了求DNA Squence Aligment写的。

我的程序代码:

#include "iostream"

using namespace std;

int L[100][100];
int n;
char DNA1[100],DNA2[200];
void processlcs()
{
int i,j;
for(i=0;i<=n;i++)
L[i][0]=0;
for(j=0;j<=n;j++)
L[0][j]=0;
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
{
if(DNA1[i-1]==DNA2[j-1])
L[i][j]=L[i-1][j-1]+1;
else
{
if(L[i][j-1]>L[i-1][j])
L[i][j]=L[i][j-1];
else
L[i][j]=L[i-1][j];
}
}
return;
}
int LCS(int i,int j)
{
return L[i+1][j+1];
}

void finds(int i,int j,int n1,int n2)
{

}
int main()
{
int i,j;
cin>>DNA1>>DNA2;
i=strlen(DNA1);
j=strlen(DNA2);
if(i>j)
n=i;
else
n=j;
processlcs();
finds(i-1,j-1,i,j);
cout<<LCS(i-1,j-1)<<endl;
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: