我刚刚在UVA's Online Judge遇到了这个小问题,并且认为它可能是一个很好的候选人进行一些代码高尔夫。
问题:
您将设计一个程序,以帮助建筑师根据城市中建筑物的位置绘制城市的天际线。为了使问题易于处理,所有建筑物都是矩形的,并且它们共用一个共同的底部(它们内置的城市非常平坦)。这个城市也被视为二维的。建筑物由有序三重(Li,Hi,Ri)指定,其中 Li 和 Ri 分别是建筑物的左右坐标我和你好是建筑物的高度。
在下图中,建筑物左侧显示三元组
(1,11,5), (2,6,7), (3,13,9), (12,7,16), (14,3,25), (19,18,22), (23,13,29), (24,4,28)
,右侧显示的天际线由序列表示:
1, 11, 3, 13, 9, 0, 12, 7, 16, 3, 19, 18, 22, 3, 23, 13, 29, 0
输出应包含描述天际线的矢量,如上例所示。在天际线矢量(v1,v2,v3,... vn)中, vi 使得i是偶数表示水平线(高度)。 vi 使得i是奇数表示垂直线(x坐标)。天际线矢量应该表示所采取的“路径”,例如,从最小x坐标开始并在定义天际线的所有线上水平和垂直行进的bug。因此,天际线矢量中的最后一个条目将为0.坐标必须用空格分隔。
如果我不计算提供(测试)建筑物的声明并包括所有空格和制表符,我的解决方案在Python中的长度为 223 。
以下是精简版:
B=[[1,11,5],[2,6,7],[3,13,9],[12,7,16],[14,3,25],[19,18,22],[23,13,29],[24,4,28]]
# Solution.
R=range
v=[0 for e in R(max([y[2] for y in B])+1)]
for b in B:
for x in R(b[0], b[2]):
if b[1]>v[x]:
v[x]=b[1]
p=1
k=0
for x in R(len(v)):
V=v[x]
if p and V==0:
continue
elif V!=k:
p=0
print "%s %s" % (str(x), str(V)),
k=V
我认为我没有犯过任何错误,但如果是这样的话 - 随意批评我。
我没有太多的声誉,所以我只需要支付100美元的奖金 - 我很好奇,如果有人可以尝试解决这个问题,那么就算是80个字符。 cobbal 发布的解决方案是101 characters long,目前是最好的解决方案。
我认为,对于这类问题,80个字符是病态限制。 cobbal ,他的46个角色解决方案让我感到惊讶 - 尽管我必须承认,在我部分理解他所写的内容之前,我花了一些时间阅读他的解释。答案 0 :(得分:25)
我刚刚开始学习J,所以这是我第一次尝试打高尔夫球:
103 62 49
46个字符
b =: 8 3 $ 1 11 5 2 6 7 3 13 9 12 7 16 14 3 25 19 18 22 23 13 29 24 4 28
,(,.{&s)I.s~:_1|.s=:0,~>./(1&{*{.<:[:i.{:)"1 b
1 11 3 13 9 0 12 7 16 3 19 18 22 3 23 13 29 0
虽然我确信有人真正了解这种语言可能会缩短这一点
代码说明:
NB. list numbers up to right bound of the building ([: i. {:) 14 3 25 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 NB. compare to left bound of building and then multiply by height (1&{ * {. <: [: i. {:) 14 3 25 0 0 0 0 0 0 0 0 0 0 0 0 0 0 3 3 3 3 3 3 3 3 3 3 3 NB. apply to each row of b, note how shorter entries are padded with 0s (1&{ * {. <: [: i. {:)"1 b 0 11 11 11 11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 6 6 6 6 6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ... NB. collapse to find max, add a 0 to the end for rotate later, assign to s ] s =: 0 ,~ >./ (1&{ * {. <: [: i. {:)"1 b 0 11 11 13 13 13 13 13 13 0 0 0 7 7 7 7 3 3 3 18 18 18 3 13 13 13 13 13 13 0 NB. rotate s right 1 and then compare to s to find where height changes s ~: _1 |. s 0 1 0 1 0 0 0 0 0 1 0 0 1 0 0 0 1 0 0 1 0 0 1 1 0 0 0 0 0 1 NB. find indices of all differences I. s ~: _1 |. s 1 3 9 12 16 19 22 23 29 NB. pair each index with the height of the building there (,. {&s) I. s ~: _1 |. s 1 11 3 13 9 0 ... NB. and finally flatten the list , (,. {&s) I. s ~: _1 |. s 1 11 3 13 9 0 12 7 16 3 19 18 22 3 23 13 29 0
答案 1 :(得分:17)
Python,89个字符,也使用Triptych的5001作弊:
B=[[1,11,5],[2,6,7],[3,13,9],[12,7,16],[14,3,25],[19,18,22],[23,13,29],[24,4,28]]
x=o=0
while x<5001:
n=max([H for L,H,R in B if L<=x<R]+[0])
if n-o:print x,n,
o=n;x+=1
5001
替换max(map(max,B))+1
以允许(几乎)任意大城市留下102个字符。
修订历史记录:
答案 2 :(得分:10)
Python:115个字符
与OP一样,我不包括数据声明,但我在计算空白。
D = [(1,11,5), (2,6,7), (3,13,9), (12,7,16),
(14,3,25), (19,18,22), (23,13,29), (24,4,28)]
P=[max([0]+[h for s,h,e in D if s<=x<e])for x in range(5001)]
for i,x in enumerate(P[1:]):
if x!=P[i]:print i+1,x,
请注意,我使用OP提供的链接作为问题的确切定义。例如,假设没有建立超过5000的坐标,并且所有坐标都是正整数,我会作弊。在我看来,原帖不会受到足够的限制,因为这很有趣。
编辑:感谢John Pirie关于将列表构造折叠到打印循环中的提示。我怎么想念那个?!
修改:在决定使用原始问题中给出的完全定义后,我将range(1+max(zip(*D)[2]))
更改为range(5001)
。第一个版本将处理任意正整数的建筑物(假设它们都适合内存)。
编辑:意识到我过度复杂了。检查我的修订版。
顺便说一句 - 我有一种预感,有一种更优雅,也可能更短的方式来做到这一点。有人打败了我!答案 3 :(得分:9)
一个176字节的WinXP .COM可执行文件:
vQoAgMYQjsKO2jPAM / + 5AIDzq7QLzSE8 / 751AXQDvoQB6DkAi / noNACL2egvACn5A / 87HXYCiR2D xwLi9eviM8mZ9 / VSQQvAdfeI + rQCzSG3LFqAwjC0As0h4vbD / 9Y8CnP6D7bI / 9Y8CnPwtACR9 + UD yOvxtAvNITz / dRO0CM0hLDDDtAHNITwadfW + kAHDM / YZ / 7cgrTn4dA + L + I1E / tHo6Jr / i8folf8L 9nXozSA =
Base64编码,我used this site to encode it。解码为.com文件。程序读取stdin直到EOF,从控制台读取时为Ctrl-Z,然后将结果输出到stdout。
编辑:源代码:
mov bp,10
add dh,10h
mov es,dx
mov ds,dx
xor ax,ax
xor di,di
mov cx,8000h
rep stosw
mov ah,0bh
int 21h
cmp al,255
mov si,offset l9
je l1
mov si,offset l11
l1:
call l7
mov di,cx
call l7
mov bx,cx
call l7
sub cx,di
add di,di
l2:
cmp bx,[di]
jbe l3
mov [di],bx
l3:
add di,2
loop l2
jmp l1
l4:
xor cx,cx
l5:
cwd
div bp
push dx
inc cx
or ax,ax
jnz l5
mov dl,bh
mov ah,2
int 21h
mov bh,44
l6:
pop dx
add dl,48
mov ah,2
int 21h
loop l6
ret
l7:
call si
cmp al,10
jae l7
db 0fh, 0b6h, 0c8h
l8:
call si
cmp al,10
jae ret
mov ah,0
xchg cx,ax
mul bp
add cx,ax
jmp l8
l9:
mov ah,0bh
int 21h
cmp al,255
jne l12
mov ah,8
int 21h
l10:
sub al,48
ret
l11:
mov ah,1
int 21h
cmp al,26
jne l10
mov si,offset l12
ret
l12:
xor si,si
xor di,di
mov bh,32
l13:
lodsw
cmp ax,di
je l14
mov di,ax
lea ax,[si-2]
shr ax,1
call l4
mov ax,di
call l4
l14:
or si,si
jne l13
int 20h
像我一样使用A86进行编译。
答案 4 :(得分:5)
Python有133个字符,内存和时间效率,对数据输入没有限制
D = [(1,11,5), (2,6,7), (3,13,9), (12,7,16), (14,3,25), (19,18,22), (23,13,29), (24,4,28)]
l,T=0,zip(*D)
for x,h in map(lambda x:(x,max([y for a,y,b in D if a<=x<b]or[0])),sorted(T[0]+T[2])):
if h!=l: print x,h,
l=h
解释:
lambda x:(x,max([y for a,y,b in D if a<=x<b]or[0])
返回位置x处的位置和高度。
现在遍历由sorted(zip(*D)[0]+zip(*D)[2])
编译的已排序坐标列表,并在高度发生变化时输出。
第二个版本效率不如上面那个并且有一个坐标限制,但只使用 115个字符:
for x in range(100):
E=[max([y for a,y,b in D if a<=(x-i)<b]+[0])for i in(0,1)]
if E[0]-E[1]:print x,E[0],
答案 5 :(得分:5)
98 J个字符,默认定义(无变量名!):
,@(([:(1:,{:"1)}.~:}:)#])@((],[:>./(1&{@[*(]>:{.@[)*]<{:@[)"1)"2 0(([:<./{."1)}.[:i.@>:[:>./{:"1))
行动中:
$ jconsole s =: ,@(([:(1:,{:"1)}.~:}:)#])@((],[:>./(1&{@[*(]>:{.@[)*]<{:@[)"1)"2 0(([:<./{."1)}.[:i.@>:[:>./{:"1)) |:b =: 8 3 $ 1 11 5 2 6 7 3 13 9 12 7 16 14 3 25 19 18 22 23 13 29 24 4 28 1 2 3 12 14 19 23 24 11 6 13 7 3 18 13 4 5 7 9 16 25 22 29 28 s b 1 11 3 13 9 0 12 7 16 3 19 18 22 3 23 13 29 0
假设最左边的坐标始终为1,只有86个字符。
,@(([:(1:,{:"1)}.~:}:)#])@((],[:>./(1&{@[*(]>:{.@[)*]<{:@[)"1)"2 0([:>:[:i.[:>./{:"1))
只有76,另外假设最右边的坐标最多为99。
,@(([:(1:,{:"1)}.~:}:)#])@((],[:>./(1&{@[*(]>:{.@[)*]<{:@[)"1)"2 0&(>:i.99))
借用cobbal的一些技巧,我可以把它变成68。
[:,@({:"1#>:@i.@#,.{."1)[:1&|.[:(,.(~:_1&|.))[:>./(1&{*{.<:[:i.{:)"1
隐性定义无法与使用s=:…
竞争消除冗余。
每当我用 J 回答问题时,我会花时间解释发生了什么,因为我认为其他人可能会喜欢看外语的运作。
NB. The first, second, and last element of a vector ({. 0{b), (1 { 0{b), ({: 0{b) 1 11 5 NB. Count from 0 to (last element of vector)-1 i. {: 0{b 0 1 2 3 4 NB. Booleans: first element of vector less than or equal to (above)? ({. <: [:i.{:) 0{b 0 1 1 1 1 NB. Multiply by second element of vector (1&{ * {.<:[:i.{:) 0{b 0 11 11 11 11 NB. Stack up results for each vector, then find maximum by column >./ (1&{*{.<:[:i.{:) " 1 b 0 11 11 13 13 13 13 13 13 0 0 0 7 7 7 7 3 3 3 18 18 18 3 13 13 13 13 13 13 NB. Identify leaders and make table |: (,. (~: _1 & |.)) >./(1&{*{.<:[:i.{:)"1 b 0 11 11 13 13 13 13 13 13 0 0 0 7 7 7 7 3 3 3 18 18 18 3 13 13 13 13 13 13 1 1 0 1 0 0 0 0 0 1 0 0 1 0 0 0 1 0 0 1 0 0 1 1 0 0 0 0 0 NB. Rotate left |: 1 |. (,.(~:_1&|.))>./(1&{*{.<:[:i.{:)"1 b 11 11 13 13 13 13 13 13 0 0 0 7 7 7 7 3 3 3 18 18 18 3 13 13 13 13 13 13 0 1 0 1 0 0 0 0 0 1 0 0 1 0 0 0 1 0 0 1 0 0 1 1 0 0 0 0 0 1 NB. 1-based index and first element, when last element is true |: ({:"1 # >: @ i. @ # ,. {."1) 1&|.(,.(~:_1&|.))>./(1&{*{.<:[:i.{:)"1 b 1 3 9 12 16 19 22 23 29 11 13 0 7 3 18 3 13 0 NB. Flatten , ({:"1#>:@i.@#,.{."1)1&|.(,.(~:_1&|.))>./(1&{*{.<:[:i.{:)"1 b 1 11 3 13 9 0 12 7 16 3 19 18 22 3 23 13 29 0 NB. Rearrange for tacit verb ([:,@({:"1#>:@i.@#,.{."1)[:1&|.[:(,.(~:_1&|.))[:>./(1&{*{.<:[:i.{:)"1) b 1 11 3 13 9 0 12 7 16 3 19 18 22 3 23 13 29 0
答案 6 :(得分:5)
2 C#答案 - 方式太久了,但我希望看到更好?
LINQ方法(135个字符,不包括数组行):
var a=new[]{new[]{1,11,5},new[]{2,6,7},new[]{3,13,9},new[]{12,7,16},new[]{14,3,25},new[]{19,18,22},new[]{23,13,29},new[]{24,4,28}};
int l=0,y,x=-1;while(++x<5001){var b=a.Where(c=>c[0]<=x&&c[2]>x);if((y=b.Any()?b.Max(c=>c[1]):0)!=l)Console.Write(x+", "+(l=y)+", ");}
或非LINQ答案(179 185 字符,不包括数组行):
var a={1,11,5,2,6,7,3,13,9,12,7,16,13,3,25,19,18,22,23,13,29,24,4,28};
var b=new int[5000];int i=-1,j,z;while(++i<a.Length)for(j=a[i*3];j<a[i*3+2];j++)if((z=a[i*3+1])>b[j])b[j]=z;i=-1;z=0;while(++i<5000)if(b[i]!=z)Console.Write(i+", "+(z=b[i])+", ");
答案 7 :(得分:3)
代码被压缩(几行到代码),这对锦标赛有利(时间是最稀缺的资源),看起来是正确的(我不知道python,但我认为我理解代码)。
您的解决方案基本上在缓冲区中绘制城市天际线,然后以所需格式输出缓冲区的内容。
您从问题中忽略的额外信息是,最多将有5000个建筑物,水平位置将小于10.000。这意味着在您的情况下内存似乎不是问题(假设32位架构的天际线为40kb,建筑物描述为45kb - 可选,您可以在读取循环中绘制天际线)。该算法在建筑物数量上呈线性,因此速度很快。
由于内存限制较严格,你可以选择一次通过算法,但我相信在这种情况下它会执行得更慢并且实现起来要复杂得多(更多的时间,更多的CPU时间)
现在您应该考虑以给定格式读取输入并将该数据用于计算而不是预先存储的数据数组。
BTW,现在在ACM竞赛中python是一种有效的语言吗?答案 8 :(得分:3)
(快速,我的意思是不到两个小时)
(不包括“ #/
”以改进突出显示)
use 5.010;
$/=undef;
@s=map{[split',',$_]}grep{$_}split/\)\s*(?:$|,\s*\()|^\s*\(/,<>; #/
for$s(@s){($l,$y,$r)=@$s;
for$x($l..$r){$c=$p[$x];$p[$x]=$c>$y?$c:$y;}}
for($x=1;$x<=@p;$x++){$y=$p[$x]||0;
if(!defined$z){$l=$x;$z=$y;
}elsif($y!=$z){push@n,[$l,$z,$x-1];$z=$y;$l=$x;}}
push@n,[$l,$z];
say join', ',map{($_->[0],$_->[1])}@n;
#! /usr/bin/env perl
use strict;
use warnings;
use 5.010;
use YAML;
use List::Util 'max';
my $file;
{
local $/ = undef;
$file = <>;
}
my @sections = map { [split ',', $_] } grep {$_} split m'
\)\s* (?:$|,\s*\() |
^ \s* \(
'x, $file;
#my $max_x = max map{ $_->[2] } @sections;
#say $max_x;
my @points;
for my $reg( @sections ){
my($l,$y,$r) = @$reg;
for my $x ( $l .. $r ){
my $c = $points[$x] || 0;
$points[$x] = max $c, $y;
}
}
my @new;
my($l,$last_y);
for( my $x=1; $x <= @points; $x++ ){
my $y = $points[$x] || 0;
# start
if( ! defined $last_y ){
$l = $x;
$last_y = $y;
next;
}
if( $y != $last_y ){
push @new, [$l,$last_y,$x-1];
$last_y = $y;
$l = $x;
next;
}
}
push @new, [$l,$last_y];
say Dump \@sections, \@points, \@new;
say join ', ', map { ($_->[0],$_->[1]) } @new;
(不包括“ #/
”以改进突出显示)
#! /usr/bin/env perl
use strict;
use warnings;
use YAML;
use 5.010;
$/=undef;
my@s=map{[split',',$_]}grep{$_}split/\)\s*(?:$|,\s*\()|^\s*\(/,<>; #/
my@p;
{
no strict; no warnings 'uninitialized';
for$s(@s){
($l,$y,$r)=@$s;
for$x($l..$r){
$c=$p[$x];
$p[$x]=$c>$y?$c:$y;
}
}
}
# $last_y => $z
my @n;
{
no strict;
#my($l,$z);
for($x=1;$x<=@p;$x++){
$y=$p[$x]||0;
if(!defined$z){
$l=$x;
$z=$y;
}elsif($y!=$z){
push@n,[$l,$z,$x-1];
$z=$y;
$l=$x;
}
}
push@n,[$l,$z];
}
say Dump \@s, \@p, \@n;
say join', ',map{($_->[0],$_->[1])}@n;
我使用YAML来确保获得正确的数据,并且不同版本的工作方式相同。
答案 9 :(得分:3)
假设输入:
b=[(1,11,5),(2,6,7),(3,13,9),(12,7,16),(14,3,25),(19,18,22),(23,13,29),(24,4,28)]
Haskell:105个字符
h x=maximum$0:[y|(l,y,r)<-b,l<=x,x<r]
main=putStr$unwords[show x++" "++show(h x)|x<-[1..9999],h x/=h(x-1)]
字符串格式似乎是Haskell落后于Python解决方案的地方。必须使用额外的5个字符来编写'main ='也无济于事,但也许它不应该被包括在内,如果他们的代码必须演示完整的程序,C#/ Java解决方案将是巨大的:)
Haskell:76个字符 (没有字符串格式&amp; no main)
h x=maximum$0:[y|(l,y,r)<-b,l<=x,x<r]
print[(x,h x)|x<-[1..9999],h x/=h(x-1)]
回顾original problem它需要你从文件中读取输入,所以我认为看到添加了多少个字符会很有趣。
Haskell:149个字符(完整解决方案)
main=interact f
f i=unwords[show x++" "++show(h x)|x<-[1..9999],h x/=h(x-1)] where
h x=maximum$0:[y|[l,y,r]<-b,l<=x,x<r]
b=map(map read.words)$lines i
以下是完整解决方案的内容,尽可能使用更具描述性的变量名称和类型签名。
main :: IO ()
main = interact skyline
skyline :: String -> String
skyline input =
unwords [show x ++ " " ++ show (heightAt x) |
x <- [1..9999], heightAt x /= heightAt (x-1)]
where heightAt :: Int -> Int
heightAt x = maximum $ 0 : [h | [l,h,r] <- buildings, l <= x, x < r]
buildings :: [[Int]]
buildings = map (map read . words) $ lines input
答案 10 :(得分:2)
重读UVA规则,我们不限于最大X为5000,而是5000 建筑物。允许使用最多(包括)9999的X和Y值。
此外,显然只有C,C ++,C#和Java才是官方认可的语言,所以我确实使用Java。这些数字只是空格分隔,但可以放回逗号(以两个总字符数为代价)。总共153个字符(不包括数组行):
int[][]b=new int[][]{{1,11,5},{2,6,7},{3,13,9},{12,7,16},{14,3,25},{19,18,22},{23,13,29},{24,4,28}};
int[]y=new int[10000];int i;for(int[]o:b)for(i=o[0];i<o[2];y[i]=Math.max(y[i++],o[1]));for(i=0;i<9999;)if(y[i++]!=y[i])System.out.print(i+" "+y[i]+" ");
逻辑很简单。使流程变得有点不稳定的唯一因素是变量重用和后增量的非标准放置。生成:
1 11 3 13 9 0 12 7 16 3 19 18 22 3 23 13 29 0
答案 11 :(得分:2)
C
int main(int arc, char **argv) {
int B[][3]={{1,11,5},{2,6,7},{3,13,9},{12,7,16},{14,3,25},{19,18,22},{23,13,29},{24,4,28}},o=0,y=0,x=0,blen=8,bs=0,b;
for (;x<9001;x++,o=y,y=0) {
for (b=bs;b<blen;b++) {
if (x >= B[b][0] && x < B[b][2] && B[b][1] > y) y=B[b][1];
if (x > B[b][2]) bs = b;
}
if (y-o) printf("%d %d ", x, y);
}
}
答案 12 :(得分:2)
红宝石,80个字符
B=[[1,11,5],[2,6,7],[3,13,9],[12,7,16],[14,3,25],[19,18,22],[23,13,29],[24,4,28]]
o=0;1.upto(5001){|x|y=(B.map{|b|x>=b[0]&&x<b[2]&&b[1]||0}.max);o!=(o=y)&&p(x,y)}
答案 13 :(得分:2)
除了挑战。
结果设置是否正确? 在第22位,最高点是18,而23是13,所以3不是最高点。
我也尝试制作一个php版本,它给了我一个不同的最终向量。它没有针对速度进行优化。
<?php
$buildings = array(
array(1,11,5),
array(2,6,7),
array(3,13,9),
array(12,7,16),
array(14,3,25),
array(19,18,22),
array(23,13,29),
array(24,4,28)
);
$preSkyline = array();
for( $i = 0; $i<= 30; $i++){
foreach( $buildings as $k => $building){
if( $i >= $building[0] && $i<= $building[2]){
$preSkyline[$i][$k] = $building[1];
} else{
$preSkyline[$i][$k] = 0;
}
}
}
foreach( $preSkyline as $s =>$a){
$skyline[$s] = max( $a );
}
$finalSkyline = array();
foreach( $skyline as $k => $v){
if( $v !== $skyline[ $k-1]){
$finalSkyline[$k] = $v;
}
}
echo "<pre>";
print_r( $finalSkyline );
?>
返回:
Array
(
[0] => 11
[2] => 13
[9] => 0
[11] => 7
[16] => 3
[18] => 18
[22] => 13
[29] => 0
)
是拐点和最大高度。
答案 14 :(得分:2)
这是我在Perl中的尝试。 137个字符,其中33个用于查找天际线的末尾。
@a = ([1,11,5],[2,6,7],[3,13,9],[12,7,16],[14,3,25],[19,18,22],[23,13,29],[24,4,28]);
($x)=sort{$b<=>$a}map{$$_[2]}@a;
for(;$c<=$x;$c++){$n=0;map{$n=$$_[1]if$c>=$$_[0]&&$c<$$_[2]&&$n<$$_[1]}@a;print"$c $n "if$n!=$h;$h=$n;}
答案 15 :(得分:1)
#include <stdio.h>
#define MAX_B 5000
static unsigned max_y[MAX_B];
main() {
signed i, j;
int max_x = 0, last_new = 0, curr_ht = 0;
for (;!feof(stdin);) {
int l, h, r;
fscanf(stdin, "%d %d %d\n", &l, &h, &r);
if (r > max_x)
max_x = r;
for (i = l; i <= r; i++)
if (max_y[i] < h)
max_y[i] = h;
}
max_x += 2;
for (i = 0; i < max_x; i++) {
j = max_y[i] - last_new;
curr_ht += j;
last_new = max_y[i];
if (j > 0)
printf("%d %d ", i, curr_ht);
if (j < 0)
printf("%d %d ", i - 1, curr_ht);
}
printf("\n");
}
非常简单的C解决方案...... 540个字符。
答案 16 :(得分:1)
即使这是一篇老帖子,我还是认为我会用 137 字符分享我的gnu octave实现:
function[p]=sl(b)s=zeros(max(b)(:,2),max(b(:,3)));for i=b's(1:i(2),i(1):i(3)-1)=1;end;t=sum(s);u=find(shift(t,1)~=t);p=[u;t(u)](:)';end;
将任意大小的三元组矩阵作为b