关于破解维吉尼亚密码的问题

mathe

统计结果表示那个数据有问题，自然无法解开。
我认为要么那个数据有误，要么就是加密算法有问题，或者数据不是来源于英文。

lmax

我想问一下  这个 $chi^2$  叫什么，怎么念，英文应该叫什么？？？

BTW, 我想问一下，相减后，如果不除以Pn会怎么样呢？好像也可以呀？

就是直接 fabs (standard_freq - cipher_freq)会如何？

mathe

$chi^2$是希腊字母，应该就按这个希腊字母发音，所以就叫Chi Square。
不除以$p_n4$是也可以看成一个准则，但是不同部分权重不同，没有除以$p_n$后的好，而且最重要的是$chi^2$是一个很重要的分布，所以结果很好处理

mathe

关于$chi^2$分布，可以查看http://mathworld.wolfram.com/Chi-SquaredDistribution.html

荼靡破晓

购买了竟然不能阅读！

荼靡破晓

要是有java版的就好了。。

yy31742

金币是够的，怎么TM不能下载，草！

yy31742

看不到代码，看不到代码。。。。。。。

yy31742

靠，我购买了附件怎么不能下载呢

mathe

1. #include <stdio.h>
   
2. #include <string.h>
   
3. #include <stdlib.h>
   
4. #include <math.h>
   
5. 
   
6. const char *raw = "BZWIEZSEHPVCIEBSYHAWWYLIIGWHMILZLMHTCSXWWCWXWWJLEIWOLMBSWPXWWXSVTZPWEXVESXWWXPMHLPRXDLSMWSJPEQXZLHATOPVIQAYHMCYDLIPNPWSUUZVRDMEMRIZPJMTDOALTFDYHSWYPCBQDLIPXCSWTSYHWIGZHYEJTKLIOSMPTQZYVHZPEZTKREXWWCIHGGFRHBAYIECVMSATVOSACLZMXWADFVDLSIVHKLMHIGSMQSGJSYXFEIRSLZVIXYYSZTJFWAXDWCSJSNXYPDWCVJDPYWPFOXLTQSEXTVSMQPDWXLTEZVIQWNEYHWZJLXKOVIPELRHLZLXLTZLHWPAO";
   
7. 
   
8. #define KEY_LEN 12
   
9. #define NUMC 26
   
10. char freq[KEY_LEN][NUMC];
    
11. int count[KEY_LEN];
    
12. double prob[NUMC]={0.0788,0.0156,0.0268,0.0389,0.1268,0.0256,0.0187,0.0573,0.0707,0.0010,0.0060,0.0394,0.0244,0.0706,0.0776,0.0186,0.0009,0.0594,0.0634,0.0978,0.0280,0.0102,0.0214,0.0016,0.0202,0.0006};
    
13. int main()
    
14. {
    
15.     int len=strlen(raw);
    
16.     int key_len = 5;
    
17.     int i,j,k;
    
18.     double entropy=0.0;
    
19.     double var=0.0;
    
20.     for(i=0;i<NUMC;i++){
    
21.             entropy-=prob[i]*log(prob[i]);
    
22.             var+=prob[i]*log(prob[i])*log(prob[i]);
    
23.     }
    
24.     for(i=0;i<len;i++){
    
25.             char c=raw[i];
    
26.             if('a'<=c&&c<='z'){
    
27.                     freq[i%KEY_LEN][c-'a']++;
    
28.             }else if('A'<=c&&c<='Z'){
    
29.                     freq[i%KEY_LEN][c-'A']++;
    
30.             }else{
    
31.                     fprintf(stderr, "Invalid char %c encounted\n",c);
    
32.             }
    
33.             count[i%KEY_LEN]++;
    
34.     }
    
35.     for(i=0;i<KEY_LEN;i++){
    
36.             double best_entropy = 0;
    
37.             int best_j=-1;
    
38.             for(j=0;j<NUMC;j++){
    
39.                     double real_entropy=0.0;
    
40.                     for(k=0;k<NUMC;k++){
    
41.                             real_entropy-=freq[i][(j+k)%NUMC]*log(prob[k]);
    
42.                     }
    
43.                     if(best_j<0||real_entropy<best_entropy){
    
44.                             best_j=j;
    
45.                             best_entropy=real_entropy;
    
46.                     }
    
47.             }
    
48.             double exp_entropy = count[i]*entropy;
    
49.             double exp_var = count[i]*var;
    
50.             exp_var = sqrt(exp_var-exp_entropy*exp_entropy/count[i]);
    
51.             if(best_entropy>exp_entropy+4*exp_var){
    
52.                     fprintf(stderr, "\nNo valid solution found\n");
    
53.                     exit(-1);
    
54.             }else{
    
55.                     printf("%c",'A'+best_j);
    
56.             }
    
57.     }
    
58.     printf("\n");
    
59.     return 0;
    
60. }
    

复制代码