开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > Windows编程 > DNA > 查看源码
ShowDlg.cpp
资源名称：DNA4.rar　[点击查看]
上传用户：kanghaisi
上传日期：2017-04-20
资源大小：5727k
文件大小：10k
源码类别：

DNA

开发平台：

Visual C++

ShowDlg.cpp：源码内容
  1. // ShowDlg.cpp : implementation file
  2. //
  4. #include "stdafx.h"
  5. #include "DNA.h"
  6. #include "ShowDlg.h"
  8. #ifdef _DEBUG
  9. #define new DEBUG_NEW
  10. #undef THIS_FILE
  11. static char THIS_FILE[] = __FILE__;
  13. #define maxNode 20//最大顶点数
  14. #define codeLen  20//边的编码长度
  16. static  char code[maxNode][codeLen+1]={//边的编码
  17.     "CAACCCAAAACCTGGTAGAG",
  18. "ATATCGCGGGTTCAACGTGC",
  19. "ACGTTGACATGCAGGATCGA",
  20. "TATAGTACTGATCTGTAGTC",
  21. "AACCTGGTACCAGACTTGAC",
  22. "TGGTTTGGACTGGTCAAGTT",
  23. "TATAGCGCATGCAGGATCGA"
  24. };
  26. static int AdjMatrix[maxNode][maxNode];
  28. int arc;
  30. char s[maxNode*maxNode][maxNode*codeLen+1];//保存所有生成路径的编码
  31. char newstr[maxNode*codeLen+1];
  33. char result[maxNode][maxNode+1];//保存正确结果路径
  34. char resultold[maxNode][maxNode+1];//传统方法保存正确结果 路径
  36. int answer=0;//最终路径数
  37. int answerold=0;//传统方法的最终路径数
  39. bool flagold[maxNode];
  40. int  travalold[maxNode];
  41. //bool flag[maxNode];//标志每个顶点是否出现过
  42. //int  traval[maxNode];//保存路径
  44. #endif
  46. /////////////////////////////////////////////////////////////////////////////
  47. // CShowDlg dialog
  50. CShowDlg::CShowDlg(CWnd* pParent /*=NULL*/)
  51. : CDialog(CShowDlg::IDD, pParent)
  52. {
  53. //{{AFX_DATA_INIT(CShowDlg)
  54. m_result = _T("");
  55. m_graph = _T("");
  56. m_resultold = _T("");
  57. m_head = 0;
  58. m_tail = 0;
  59. m_vertex = 0;
  60. m_arc = 0;
  61. m_multiresult = _T("");
  62. //}}AFX_DATA_INIT
  63. }
  66. void CShowDlg::DoDataExchange(CDataExchange* pDX)
  67. {
  68. CDialog::DoDataExchange(pDX);
  69. //{{AFX_DATA_MAP(CShowDlg)
  70. DDX_Control(pDX, IDC_MULTITHREADRUN, m_multithreadrun);
  71. DDX_Control(pDX, IDC_RUNOLD, m_runold);
  72. DDX_Control(pDX, IDC_FINISH, m_finish);
  73. DDX_Control(pDX, IDC_ADD, m_add);
  74. DDX_Control(pDX, IDC_RUNOK, m_runok);
  75. DDX_Text(pDX, IDC_RESULT, m_result);
  76. DDX_Text(pDX, IDC_GRAPH, m_graph);
  77. DDX_Text(pDX, IDC_RESULTOLD, m_resultold);
  78. DDX_Text(pDX, IDC_HEAD, m_head);
  79. DDX_Text(pDX, IDC_TAIL, m_tail);
  80. DDX_Text(pDX, IDC_VERTEX, m_vertex);
  81. DDX_Text(pDX, IDC_ARC, m_arc);
  82. DDX_Text(pDX, IDC_MULTIRESULT, m_multiresult);
  83. //}}AFX_DATA_MAP
  84. }
  87. BEGIN_MESSAGE_MAP(CShowDlg, CDialog)
  88. //{{AFX_MSG_MAP(CShowDlg)
  89. ON_BN_CLICKED(IDC_RUNOK, OnRunok)
  90. ON_BN_CLICKED(IDC_ADD, OnAdd)
  91. ON_BN_CLICKED(IDC_FINISH, OnFinish)
  92. ON_BN_CLICKED(IDC_RUNOLD, OnRunold)
  93. //}}AFX_MSG_MAP
  94. END_MESSAGE_MAP()
  96. /////////////////////////////////////////////////////////////////////////////
  97. // CShowDlg message handlers
  101. void CShowDlg::init()
  102. {
  103. char fileName[1][10];//文件名,保存形成的路径
  105. int u,v;  
  106. char str[2*codeLen+1];
  108. FILE *fpTest,*fp2;
  110. strcpy(fileName[0],"0.txt");//将0.txt赋值给fileName[0]
  111. if(!(fpTest=fopen("Test.txt","r"))) 
  112. {
  113. MessageBox("不能打开Test.txt文件！");
  114. return ;
  115. }   
  116.     if(!(fp2=fopen(fileName[0],"w"))) 
  117. {
  118. MessageBox("不能打开0.txt文件！");
  119. return ;
  120. } 
  121.     
  122. //将test.txt文件中的点的个数和边的条数分别赋值给nodeN和edgeN
  123.    // fscanf(fpTest,"%d%d",&nodeN,&edgeN);
  125.     while(arc--)
  126.     {
  127.         fscanf(fpTest,"%d%d",&u,&v); //将边的信息依次读入u,v中
  129. char temp[codeLen+1];
  130. int k;
  132. //对边进行编码
  133. if(u!=0)//u不是起点
  134. {
  135. //顶点编码的后部分
  136. for(k=0;k<codeLen/2;k++)
  137. temp[k]=code[u][k+codeLen/2];
  138. temp[k]='';
  139. strcpy(str,temp);
  140. }
  141. else
  142.             strcpy(str,code[u]);//起点的全部编码
  144. if(v!=m_vertex-1)//v不是终点
  145. {
  146. //编码的前部分
  147. for(k=0;k<codeLen/2;k++)
  148. temp[k]=code[v][k];
  149. temp[k]='';
  150. strcat(str,temp);
  151. }
  152. else
  153. strcat(str,code[v]);
  154.         fprintf(fp2,"%sn",str);//将str的内容写入fp2中      
  155.     } 
  156.     fclose(fpTest);
  157.     fclose(fp2); 
  159.     return ;
  160. }
  162. int CShowDlg::getid(char *str, int n)
  163. {
  164. int i;
  166. for(i=0;i<n;i++)
  167. if(strcmp(str,code[i])==0)
  168. {
  169. return i;
  170. }
  172. return -1;
  173. }
  175. bool CShowDlg::str_compare(char *str, int n)
  176. {
  177. bool flag[maxNode];//标志每个顶点是否出现过
  178.     int  traval[maxNode];//保存路径
  179. char s[codeLen+1];
  180. int i,k,j,len;
  182. for(i=0;i<codeLen;i++)
  183. s[i]=str[i];    
  184. s[i]='';
  186. //如果str第一个顶点不是顶点0，则舍弃
  187. if(strcmp(s,code[0])!=0)
  188. return false;
  190. len=strlen(str);
  191. k=0; 
  193. for(i=len-codeLen;i<len;i++)
  194. s[k++]=str[i] ;
  195. s[k]='';
  197. //如果str最后一个顶点不是顶点nodeN-1，则舍弃
  198. if(strcmp(s,code[n-1])!=0)
  199. return false;
  201. //将所有的flag均设置为false
  202. memset(flag,false,sizeof(flag));
  203. flag[0]=flag[n-1]=true;
  204. k=0;
  205. j=0;
  206. traval[j++]=0;
  208. for(i=codeLen;i<len-codeLen;i++)
  209. {
  210. s[k++]=str[i];
  211. if(k==codeLen)
  212. {
  213. s[k]='';
  214. k=0;   
  215. int t=getid(s,n); 
  216. flag[t]=true;
  217. traval[j++]=t;  
  218. } 
  220. }
  222. for(i=0;i<n;i++)
  223. if(!flag[i])
  224. return false;
  226. traval[n-1]=n-1;
  228. {   
  229. for(i=0;i<n;i++)
  230.              result[answer][i]=traval[i];
  231. answer++;
  232. }               
  234. return true;     
  235. }
  238. int CShowDlg::head(char *str, char *strCode)//判断编码的头部是否匹配
  239. {
  241.     int i;
  242.     char temp[11],temps[11];
  244.     for(i=0;i<codeLen/2;i++)
  245. temp[i]=str[i];
  246. temp[i]='';
  248.     for(i=0;i<codeLen/2;i++)
  249. temps[i]=strCode[i+codeLen/2];
  250. temps[i]='';
  252.     return strcmp(temp,temps)==0;     
  253. }
  255. int CShowDlg::tail(char *str, char *strCode)
  256. {
  257.     int i,k,len;
  258. char temp[11],temps[11];
  259. len=strlen(str);
  260. for(i=len-codeLen/2,k=0;i<len;i++)
  261. temp[k++]=str[i];
  262. temp[k]='';
  264. len=strlen(strCode);
  265. for(i=0,k=0;i<codeLen/2;i++)
  266. temps[k++]=strCode[i];
  267. temps[k]='';
  268. return strcmp(temp,temps)==0;   
  269. }
  273. void CShowDlg::connect(int n)
  274. {
  275.     FILE *fp1,*fpi,*fpj;
  276. int i,j,m;
  278. char fileName[maxNode][15];
  279. char select[codeLen+1];//加入的点 
  281. for(i=0;i<n;i++)
  282. {
  283. fileName[i][0]=i+'0';
  284. fileName[i][1]='';
  285. strcat(fileName[i],".txt");                   
  286. }
  288. char str[maxNode*codeLen+1];
  289. bool flaghead,flagtail;
  291. fp1=fopen(fileName[0],"r");
  292. i=0;
  293. //读入边
  294. while(!feof(fp1))
  295. {
  296. fscanf(fp1,"%s",s[i]);
  297. i++;
  298. }
  299. m=i-1;
  300.     fclose(fp1);
  302. for(i=1;i<n-1;i++)
  303. {
  304. strcpy(select,code[i]);//将第i个顶点的编码赋值给select
  305. fpi=fopen(fileName[i-1],"r");
  306. fpj=fopen(fileName[i],"w");
  307. while(!feof(fpi))
  308. {
  309. fscanf(fpi,"%s",str);//将长度为i-1的串读入str中
  310. if(feof(fpi))break;
  312. flaghead=flagtail=false;
  314. if(head(str,select))//str的头部与第i个顶点的编码匹配
  315. flaghead=true;
  317. if(tail(str,select))//str的尾部与第i个顶点的编码匹配
  318. flagtail=true;
  320. for(j=0;j<m;j++)
  321. {
  322. if(head(s[j],select)&&flagtail)
  323. {
  324. strcpy(newstr,str);
  325. strcat(newstr,s[j]);
  326. fprintf(fpj,"%sn",newstr);                           
  327. }
  329. if(tail(s[j],select)&&flaghead)
  330. {
  333. strcpy(newstr,s[j]);
  334. strcat(newstr,str);
  335. fprintf(fpj,"%sn",newstr);                           
  336. }              
  337. }              
  338. }
  339. fclose(fpi);
  340. fclose(fpj);              
  341. }
  343. return ;  
  344. }
  346. void CShowDlg::gettravel(int n)
  347. {
  348. char str[maxNode*codeLen+1];
  349. FILE *fp;
  350. int i;
  351. char fileName[maxNode][10];
  353. //将文件名分别赋值为i.txt(i=0,1,...,n-1)
  354. for(i=0;i<n;i++)
  355. {
  356. fileName[i][0]=i+'0';
  357. fileName[i][1]='';
  358. strcat(fileName[i],".txt");                   
  359. }
  361. fp=fopen(fileName[n-2],"r");
  363. while(!feof(fp))
  364. {
  365. fscanf(fp,"%s",str);
  367. if(feof(fp))break;
  369. if(strlen(str)!=n*codeLen)
  370. continue;
  372. str_compare(str,n);
  374. }
  375. fclose(fp);
  376. return ;
  377. }
  380. void CShowDlg::OnRunok() 
  381. {
  382. // TODO: Add your control notification handler code here
  383. int i,j,n;
  384. n=m_vertex;
  385.      init(); 
  386. connect(m_vertex);
  387. gettravel(m_vertex);
  388. m_result="";
  389. CString strTemp;
  390. for(i=0;i<answer;i++)
  391. {
  392. for(j=0;j<n;j++)
  393. {
  394. strTemp.Format("%d",result[i][j]);
  395. m_result+=strTemp;
  396. m_result+=" ";
  397. }
  398.             m_result+="rn";
  399.              i++;
  400. }
  402. UpdateData(FALSE);
  403. }
  405. void CShowDlg::OnAdd() 
  406. {
  407. // TODO: Add your control notification handler code here
  408. UpdateData(TRUE);
  410. FILE *fpTest;
  411.      if(!(fpTest=fopen("Test.txt","a"))) 
  412. {
  413. MessageBox("不能打开Test.txt文件！");
  414. }  
  416.  fprintf(fpTest,"%d %dn",m_head,m_tail);
  417.      fclose(fpTest);
  418.  
  419.     // fscanf(fpi,"%s",str);//将长度为i-1的串读入str中
  420. UpdateData(FALSE);
  421.     }
  423. void CShowDlg::OnFinish() 
  424. {
  425. // TODO: Add your control notification handler code here
  426. UpdateData(TRUE);
  427. int u,v,n;
  428. FILE *fpTest;
  429.      if(!(fpTest=fopen("Test.txt","r"))) 
  430. {
  431. MessageBox("不能打开Test.txt文件！");
  432. }  
  433. n=m_vertex;
  434. arc=m_arc;
  435. for(u=0;u<m_vertex;u++)
  436. for(v=0;v<m_vertex;v++)
  437. AdjMatrix[u][v]=0;
  438.       while(!feof(fpTest)) 
  439.   {
  440.   fscanf(fpTest,"%d%d",&u,&v);
  441.   AdjMatrix[u][v]=1;
  442.   }     
  443.   fclose(fpTest);
  444.    m_graph="";
  445. CString strTemp;
  446. for(u=0;u<m_vertex;u++)
  447. {
  448. for(v=0;v<m_vertex;v++)
  449. {
  450. strTemp.Format("%d",AdjMatrix[u][v]);
  451. m_graph+=strTemp;
  452. m_graph+=" ";
  453. }
  454.             m_graph+="rn";         
  455. }
  456.  UpdateData(FALSE);
  457. }
  459. void CShowDlg::OnRunold() 
  460. {
  461. // TODO: Add your control notification handler code here
  462.     int i,j;
  463.     solve();
  464. m_resultold="";
  465. CString strTemp;
  466. for(i=0;i<answerold;i++)
  467. {
  468. for(j=0;j<m_vertex;j++)
  469. {
  470. strTemp.Format("%d",resultold[i][j]);
  471. m_resultold+=strTemp;
  472. m_resultold+=" ";
  473. }
  474.             m_resultold+="rn";             
  475. }
  477. UpdateData(FALSE);
  478.     return ;  
  480. }
  482. void CShowDlg::Tradtion(int t,int d)
  483. {
  484.     int i; 
  485.   travalold[d]=t;   
  486.    if(d==m_vertex-1)
  487.    {
  488.       if(travalold[m_vertex-1]==m_vertex-1)
  489.       {
  490.       for(i=0;i<m_vertex;i++)
  491.       resultold[answerold][i]=travalold[i];    
  492.       }  
  493.   answerold++;
  494.       return ;     
  495.    }
  496.   
  497.    for(i=1;i<m_vertex;i++)
  498.    if(!flagold[i]&&AdjMatrix[t][i])
  499.    {
  500.       flagold[i]=true;
  501.       Tradtion(i,d+1);
  502.       flagold[i]=false;            
  503.    }
  504.    return ;
  505. }
  507. void CShowDlg::solve()
  508. {
  509.    int i;
  510.    memset(flagold,false,sizeof(flagold)); 
  511.    flagold[0]=true;
  512.    travalold[0]=0;
  513.    for(i=1;i<m_vertex;i++)
  514.    if(AdjMatrix[0][i])
  515.    {
  516.    flagold[i]=true;
  517.    Tradtion(i,1);
  518.    flagold[i]=false;
  519.    }
  520.    return ;    
  521. }