怎么获得最大值和最大值所在日期之前的最小值

alunfirst

如下数据：想得到10月31至12月9之间a的最大值和b的最大值，和最大值所在日期之前的最小值

2016/10/31	1	a	39.9
2016/11/1	1	a	41.72
2016/11/2	1	a	43.06
2016/11/3	1	a	42.7
2016/11/4	1	a	46.97
2016/11/7	1	a	45.33
2016/11/8	1	a	2
2016/11/9	1	a	38.1
2016/11/10	1	a	45.3
2016/11/11	1	a	44.93
2016/11/14	1	a	45.6
2016/11/15	1	a	45.7
2016/11/16	1	a	47.49
2016/11/17	1	a	46.88
2016/11/18	1	a	46.12
2016/11/21	1	a	100
2016/11/22	1	a	46.55
2016/11/23	1	a	46.08
2016/11/24	1	a	45
2016/11/25	1	a	43.92
2016/11/28	1	a	43.4
2016/11/29	1	a	41.8
2016/11/30	1	a	42.18
2016/12/1	1	a	41.82
2016/12/2	1	a	41.35
2016/12/5	1	a	41.71
2016/12/6	1	a	42.32
2016/12/7	1	a	43.56
2016/12/8	1	a	42.98
2016/12/9	1	a	1
2016/10/31	2	b	7.36
2016/11/1	2	b	100
2016/11/2	2	b	7.32
2016/11/3	2	b	7.38
2016/11/4	2	b	7.4
2016/11/7	2	b	7.53
2016/11/8	2	b	7.49
2016/11/9	2	b	7.43
2016/11/10	2	b	7.53
2016/11/11	2	b	7.56
2016/11/14	2	b	7.53
2016/11/15	2	b	7.61
2016/11/16	2	b	7.57
2016/11/17	2	b	7.55
2016/11/18	2	b	7.56
2016/11/21	2	b	7.54
2016/11/22	2	b	7.59
2016/11/23	2	b	7.61
2016/11/24	2	b	7.55
2016/11/25	2	b	7.56
2016/11/28	2	b	7.63
2016/11/29	2	b	7.48
2016/11/30	2	b	7.42
2016/12/1	2	b	7.45
2016/12/2	2	b	7.24
2016/12/5	2	b	7.2
2016/12/6	2	b	7.19
2016/12/7	2	b	7.26
2016/12/8	2	b	7.17
2016/12/9	2	b	7.19

得到数据：

2016/11/21

1

a

45.93

2016/11/8

1

a

2

2016/11/1

2

b

100

2016/10/31

2

b

7.36

zy86416779

1. [root@CentOS_6 tmp]# head temp.txt
   
2. 2016/10/31
   
3. 
   
4. 1
   
5.         a
   
6. 39.9
   
7. 2016/11/1
   
8. 
   
9. 1
   
10.         a
    
11. 41.72
    
12. [root@CentOS_6 tmp]# sed -i 's/\t//g' temp.txt
    
13. [root@CentOS_6 tmp]# head temp.txt            
    
14. 2016/10/31
    
15. 
    
16. 1
    
17. a
    
18. 39.9
    
19. 2016/11/1
    
20. 
    
21. 1
    
22. a
    
23. 41.72
    
24. [root@CentOS_6 tmp]# sed -i "/^$/d" temp.txt
    
25. [root@CentOS_6 tmp]# head temp.txt         
    
26. 2016/10/31
    
27. 1
    
28. a
    
29. 39.9
    
30. 2016/11/1
    
31. 1
    
32. a
    
33. 41.72
    
34. 2016/11/2
    
35. 1
    
36. [root@CentOS_6 tmp]# cat test.awk
    
37. {
    
38.   printf $0
    
39.   if(NR%4==0){
    
40.     printf "\n"
    
41.     next
    
42.   }
    
43.   printf "\t"
    
44. }
    
45. [root@CentOS_6 tmp]# awk -f test.awk temp.txt > after_temp.txt
    
46. [root@CentOS_6 tmp]# heed after_te
    
47. -bash: heed: command not found
    
48. [root@CentOS_6 tmp]# heed after_temp.txt
    
49. -bash: heed: command not found
    
50. [root@CentOS_6 tmp]# head after_temp.txt  
    
51. 2016/10/31      1       a       39.9
    
52. 2016/11/1       1       a       41.72
    
53. 2016/11/2       1       a       43.06
    
54. 2016/11/3       1       a       42.7
    
55. 2016/11/4       1       a       46.97
    
56. 2016/11/7       1       a       45.33
    
57. 2016/11/8       1       a       2
    
58. 2016/11/9       1       a       38.1
    
59. 2016/11/10      1       a       45.3
    
60. 2016/11/11      1       a       44.93
    
61. [root@CentOS_6 tmp]# cat after_test.awk
    
62. BEGIN{i=1;j=1}
    
63. {
    
64.   if(i==1){
    
65.     a = split($0, array)
    
66.         i=0
    
67.         next
    
68.   }
    
69.   if($3==array[3]){
    
70.     if($4<=array[4])
    
71.       a = split($0, array)
    
72.   }
    
73.   else if(j==1){
    
74.     b = split($0, brray)
    
75.         j=0
    
76.         next
    
77.   }
    
78.   else{
    
79.     if($3==brray[3]){
    
80.       if($4>=brray[4])
    
81.       b = split($0, brray)
    
82.     }
    
83.   }
    
84. }
    
85. END{
    
86.   for (i = 1; i <= a; ++i){
    
87.     printf array[i]
    
88.     printf "\t"
    
89.   }
    
90.   printf "\n"
    
91.   for (i = 1; i <= b; ++i){
    
92.     printf brray[i]
    
93.     printf "\t"
    
94.   }
    
95.   printf "\n"
    
96. }
    
97. [root@CentOS_6 tmp]# awk -f after_test.awk after_temp.txt  
    
98. 2016/12/9       1       a       1
    
99. 2016/11/1       2       b       100

复制代码

jason680

回复 1# alunfirst

$ awk '{k=$3;if(t[k]>$4){t[k]=$4;u[k]=$0};if(M[k]==""){m[k]=M[k]=t[k]=$4;u[k]=a[k]=b[k]=$0;}if(M[k]<$4){if(m[k]>t[k]){m[k]=t[k];b[k]=u[k]};t[k]=M[k]=$4;a[$3]=$0}}END{for(n in a){print a[n],b[n]}}' file
2016/11/21     1     a     100      2016/11/8     1     a     2
2016/11/1     2     b     100      2016/10/31     2     b     7.36     

alunfirst

jason680 发表于 2016-12-23 10:59
回复 1# alunfirst

$ awk '{k=$3;if(t[k]>$4){t[k]=$4;u[k]=$0};if(M[k]==""){m[k]=M[k]=t[k]=$4;u[k]= ...

好复杂,谢谢高手。

elu_ligao

回复 3# jason680

讲解下吧， 没看明白啊， 写个伪代码也行， 谢谢， 学习下

sunzhiguolu

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. 
   
5. my ($da, $db) = ('2016/10/31', '2016/12/9');
   
6. my @aData = ();
   
7. 
   
8. while (<DATA>){
   
9.     next if ($_ lt $da or $_ gt $db);
   
10.     my $pos = /a/ ? 0 : 1;
    
11.     push (@{$aData[$pos]}, [split]);
    
12. }
    
13. 
    
14. foreach (@aData){
    
15.     my @aT = sort {$a->[-1] <=> $b->[-1]} @$_;
    
16.     my $raMax = $aT[-1];
    
17.     @aT = sort {$a->[-1] <=> $b->[-1]} grep {$_->[0] =~ s/\b(?=\d\b)/0/gr lt $raMax->[0] =~ s/\b(?=\d\b)/0/gr} @aT;
    
18.     print join ("\t", @{$raMax}, @{$aT[0]}), "\n";
    
19. }
    
20. 
    
21. __DATA__
    
22. 2016/10/31  1   a   39.9
    
23. 2016/11/1   1   a   41.72
    
24. 2016/11/2   1   a   43.06
    
25. 2016/11/3   1   a   42.7
    
26. 2016/11/4   1   a   46.97
    
27. 2016/11/7   1   a   45.33
    
28. 2016/11/8   1   a   2
    
29. 2016/11/9   1   a   38.1
    
30. 2016/11/10  1   a   45.3
    
31. 2016/11/11  1   a   44.93
    
32. 2016/11/14  1   a   45.6
    
33. 2016/11/15  1   a   45.7
    
34. 2016/11/16  1   a   47.49
    
35. 2016/11/17  1   a   46.88
    
36. 2016/11/18  1   a   46.12
    
37. 2016/11/21  1   a   100
    
38. 2016/11/22  1   a   46.55
    
39. 2016/11/23  1   a   46.08
    
40. 2016/11/24  1   a   45
    
41. 2016/11/25  1   a   43.92
    
42. 2016/11/28  1   a   43.4
    
43. 2016/11/29  1   a   41.8
    
44. 2016/11/30  1   a   42.18
    
45. 2016/12/1   1   a   41.82
    
46. 2016/12/2   1   a   41.35
    
47. 2016/12/5   1   a   41.71
    
48. 2016/12/6   1   a   42.32
    
49. 2016/12/7   1   a   43.56
    
50. 2016/12/8   1   a   42.98
    
51. 2016/12/9   1   a   1
    
52. 2016/10/31  2   b   7.36
    
53. 2016/11/1   2   b   100
    
54. 2016/11/2   2   b   7.32
    
55. 2016/11/3   2   b   7.38
    
56. 2016/11/4   2   b   7.4
    
57. 2016/11/7   2   b   7.53
    
58. 2016/11/8   2   b   7.49
    
59. 2016/11/9   2   b   7.43
    
60. 2016/11/10  2   b   7.53
    
61. 2016/11/11  2   b   7.56
    
62. 2016/11/14  2   b   7.53
    
63. 2016/11/15  2   b   7.61
    
64. 2016/11/16  2   b   7.57
    
65. 2016/11/17  2   b   7.55
    
66. 2016/11/18  2   b   7.56
    
67. 2016/11/21  2   b   7.54
    
68. 2016/11/22  2   b   7.59
    
69. 2016/11/23  2   b   7.61
    
70. 2016/11/24  2   b   7.55
    
71. 2016/11/25  2   b   7.56
    
72. 2016/11/28  2   b   7.63
    
73. 2016/11/29  2   b   7.48
    
74. 2016/11/30  2   b   7.42
    
75. 2016/12/1   2   b   7.45
    
76. 2016/12/2   2   b   7.24
    
77. 2016/12/5   2   b   7.2
    
78. 2016/12/6   2   b   7.19
    
79. 2016/12/7   2   b   7.26
    
80. 2016/12/8   2   b   7.17
    
81. 2016/12/9   2   b   7.19
    

复制代码

sunzhiguolu

perl abc.pl
------------------------------------------------------------------------------------
2016/11/21      1       a       100     2016/11/8       1       a       2
2016/11/1       2       b       100     2016/10/31      2       b       7.36

elu_ligao

3楼代码明了了，学着写了一个

1. awk '!M[$3]{M[$3]=m[$3]=$4;PM[$3]=Pm[$3]=p[$3]=$0;next}{if($4>M[$3]){M[$3]=$4;PM[$3]=$0;Pm[$3]=p[$3]}else if($4<M[$3]){if($4<m[$3]){m[$3]=$4;p[$3]=$0}}}END{for(i in M)print PM[i]"\t"Pm[i]}'

复制代码

moperyblue

日期是不是已经排序?

1. awk '{x=mktime(gensub("/",FS,"g",$1)" 0 0 0");k=$3}NR==FNR{if($4>a[k]){a[k]=$4;max[k]=$0;b[k]=x};next}x<b[k]{if(!c[k]||$4<c[k]){c[k]=$4;min[k]=$0}}END{for(i in a)print max[i]"\t"min[i]}' file file

复制代码