foreach循环

pony2001mx

大家好：
我想输出所有长20个字母的字符串，只能为a，b，c和d，其中a、b、c的总数小于等于7。以下foreach循环脚本哪里有误？谢谢！

1. #!/usr/bin/perl
   
2. 
   
3. @a = (a,b,c,d);
   
4. $x='';
   
5. 
   
6. for $i1 (@a){
   
7. $x .= $i1;
   
8. for $i2 (@a){
   
9. $X .= $i2;
   
10. for $i3 (@a){
    
11. $x .= $i3;
    
12. for $i4 (@a){
    
13. $x .= $i4;
    
14. for $i5 (@a){
    
15. $x .= $i5;
    
16. for $i6 (@a){
    
17. $x .= $i6;
    
18. for $i7 (@a){
    
19. $x .= $i7;
    
20. for $i8 (@a){
    
21. $x .= $i8; $no = &abc($x); next if $no > 7;
    
22. for $i9 (@a){
    
23. $x .= $i9; $no = &abc($x); next if $no > 7;
    
24. for $i10 (@a){
    
25. $x .= $i10; $no = &abc($x); next if $no > 7;
    
26. for $i11 (@a){
    
27. $x .= $i11; $no = &abc($x); next if $no > 7;
    
28. for $i12 (@a){
    
29. $x .= $i12; $no = &abc($x); next if $no > 7;
    
30. for $i13 (@a){
    
31. $x .= $i13; $no = &abc($x); next if $no > 7;
    
32. for $i14 (@a){
    
33. $x .= $i14; $no = &abc($x); next if $no > 7;
    
34. for $i15 (@a){
    
35. $x .= $i15; $no = &abc($x); next if $no > 7;
    
36. for $i16 (@a){
    
37. $x .= $i16; $no = &abc($x); next if $no > 7;
    
38. for $i17 (@a){
    
39. $x .= $i17; $no = &abc($x); next if $no > 7;
    
40. for $i18 (@a){
    
41. $x .= $i18; $no = &abc($x); next if $no > 7;
    
42. for $i19 (@a){
    
43. $x .= $i19; $no = &abc($x); next if $no > 7;
    
44. for $i20 (@a){
    
45. $x .= $i20; $no = &abc($x); next if $no > 7;
    
46. print "$x\n";
    
47. }}}}}}}}}}}}}}}}}}}}    #20 brackets
    
48. 
    
49. sub abc
    
50. {
    
51. $string = shift;
    
52. @string = split//, $string;
    
53. for $elem (@string){
    
54. if ($elem eq 'a') {$ha{$elem} +=1}
    
55. elsif ($elem eq 'b') {$hb{$elem} += 1}
    
56. elsif ($elem eq 'c') {$hc{$elem} += 1}
    
57. elsif ($elem eq 'd') {$hd{$elem} += 1}
    
58. }
    
59. 
    
60. return $ha{'a'}+$hb{'b'}+$hc{'c'};
    
61. }

复制代码

tadegenban

写这种程序是为了好玩吗？

pony2001mx

不是的，其实我上面只举了个例子，我主要是做基因分析的。我的perl水平非常有限，所以请高手指点。谢谢！

fender0107401

这代码好暗黑啊。。。

xiumu2280

你还是写出需求 重新写吧··· 这个真心没法看···

104359176

哈哈，典型的把自己绕进去的算法。就是排列组合。

stanley_tam

代码重复的略多哦。
把a b c d分别当作4进制里的0，1，2，3。 那么问题就变成了从4进制数0一直数到'3' x 20。不知道对不对。。。{:3_188:}
可以用模块Math::Base::Convert，它把4进制数弄成了 用字符表示的dna进制  => a c t g。可以把dna进制替换成abcd
a、b、c的总数小于等于7，那就是d的数量大于等于13个咯
没时间跑完，你试下。{:3_204:}

1. #!perl
   
2. use Modern::perl;
   
3. use FileHandle;
   
4. use Math::Base::Convert qw{ cnv };
   
5. 
   
6. sub count_d;
   
7. sub main;
   
8. 
   
9. main;
   
10. 
    
11. sub main {
    
12.     my $number = 0;
    
13.     my $out_fh = FileHandle->new('out.txt', 'w');
    
14. 
    
15.     while (1) {
    
16.         $number++;
    
17.         my $converted = cnv($number, 10 => ['a', 'b', 'c', 'd']); # convert to string like aaaaaaaaaaaaaabdbdcd
    
18. 
    
19.         my $len = length $converted;
    
20.         if ( $len < 20) {
    
21.             # body...
    
22.             $converted = 'a' x (20-$len) . $converted;
    
23.         }
    
24.         elsif ($len > 20){
    
25.             last
    
26.         }
    
27. 
    
28.         my $d_num = count_d($converted);
    
29.         if ($d_num >= 13) {
    
30.             $out_fh->print("$converted\n");
    
31.         }
    
32.     }
    
33. 
    
34.     close $out_fh;
    
35. }
    
36. 
    
37. sub count_d {
    
38.     my $str   = shift;
    
39.     my $g_num = 0;
    
40. 
    
41.     for my $i ( unpack("U*", $str)){
    
42.         $g_num++ if $i == 100;  # ord('d') => 100
    
43.     }
    
44. 
    
45.     return $g_num;
    
46. }
    
47. 
    
48. __END__

复制代码

afukada

回复 1# pony2001mx


先對結果做預估

a+b+c=0:

N_0=C(20,0) * 3^0 = 1 * 1 = 1

a+b+c=1:

N_1=C(20,1) * 3^1 = 20 * 3 = 60

...

所以可以得到當a+b+c=n:

N_n=C(20,n) * 3^n

把結果列表出來:

a+b+c	N_n
0	1
1	60
2	1710
3	30780
4	392445
5	3767472
6	28256040
7	169536240
Total	201984748

可以檢查一下 (希望沒有算錯)

huang6894

回复 3# pony2001mx


    这样写，你确定你跑生物大数据没把pc跑死麽。。。

huang6894

以下代码，我把电脑跑死了~
4^20 = 1 099 511 627 776

1. my @out;
   
2. sub X {
   
3.     @_ ? map { my $h = $_; map [ $h, @$_ ], X(@_) } @{ +shift } : [];
   
4.     }
   
5.     my @A = ([ A,T,C,G ]) x 20;
   
6.     push @out,"@$_" for X @A;
   
7.     for my $i (@out){
   
8.             my $count = 0;
   
9.             ++$count while($i =~ m/G/g);
   
10.             print "$i\n" if $count > 13;
    
11.     }
    

复制代码