Ricostruzione
di polyomini L-convessi
G. Castiglione, A. Restivo, R. Vaglica.
Dipartimento di Matematica e Applicazioni, Università di Palermo
1
Polyomino
insieme finito e connesso di celle (quadrati unitari)
del piano cartesiano definito a meno di una
traslazione
• Polyomino convesso : polyomino le cui righe e colonne
sono connesse
Polyomino convesso
due celle di un insieme discreto si dicono connesse se esiste
un cammino (sequenza di celle adiacenti) che le collega
2
contenuto nell’insieme.
Polyomini L-convessi
polyomini convessi in cui ogni coppia di celle risulta connessa
da un cammino monotono con al più un solo cambiamento di
direzione (L-path)
L-convesso
L-convesso
• Cammino monotono: cammino autoevitante costituito
da passi in due sole direzioni
In un polyomino convesso qualunque coppia di celle
3
risulta connessa da almeno un cammino monotono.
Ricostruzione di insiemi discreti a partire da
informazioni parziali
Proiezioni orizzontali e verticali
(Tomografia discreta)
V2 2 3 8 7 7 3 3
1
3
3
8
8
6
3
3
H
L-cammini
{ 1, 2, ... ,}
4
Ricostruzione di polyomini L-convessi
Proiezioni orizzontali e
verticali
(Tomografia discreta)
unicità
L-cammini
L-cammini
massimali
L-cammini
bordati
unicità
unicità
Ricostruzione banale
(L-convessi)
Algoritmo di ricostruzione
5
L-cammini
• L(P): insieme degli L-cammini contenuti nel polyomino
L-convesso P
Denotiamo con x,y (x,yℤ-{0}) un L-cammino fatto da |x|-1
passi orizzontali, |y|-1 passi verticali orientato come segue
•
se x>0 e y>0
•
se x>0 e y<0
•
se x<0 e y>0
•
se x<0 e y<0
L(P) { x,y / x,y P }
Esempio
L-cammino 3,4
contenuto in un
polyomino L-convesso
6
L-cammini massimali
(L(P), ) è un insieme parzialmente ordinato
• x,y è massimale (in L(P)) se
x',y' L(P) , x,y x',y' x x' y y'
Lmax(P): insieme di tutti gli elementi massimali di L(P)
Osservazione: gli elementi di Lmax(P) possono avere più occorrenze
in P. Tutte queste occorrenze connettono celle appartenenti ai bordi
di P
Relazione tra altezza e larghezza di
P e l’insieme Lmax(P)
h( P) max{ y \ x, y Lmax ( P)}
w( P) max{ x \ x, y Lmax ( P)}
7
Rettangoli massimali
[x,y] Polyomino di forma rettangolare avente larghezza x ed
altezza y.
parzialmente ordinato
R(P) { [x,y] t.c. [x,y] P }
Rmax(P): insieme di tutti gli elementi massimali di R(P)
Osservazione: Rmax(P) è un insieme finito di rettangoli non
confrontabili ovvero [x,y], [x',y']Rmax(P) tali
che
[x,y] [x',y'] [x',y'] [x,y]
Rmax(P) {[x1,y1] , [x2,y2] , … , [xn,yn]}
ordinamento canonico
x1 x2 … xn and y1 y2 … yn
8
Rettangoli non confrontabili in posizione crossing
Dati i rettangoli [x,y], [x',y'] non confrontabili si dice che essi si trovano
in una posizione crossing se [ x, y ] [ x, y] [min( x, x), min( y, y)]
rettangoli in posizione crossing
rettangoli in posizione
non crossing
Teorema: Un polyomino convesso P è un L-convesso se e solo se i suoi
rettangoli massimali hanno a due a due una posizione crossing in P. 9
I polyomini L-convessi sono caratterizzati
dai loro L-cammini massimali?
Rmax(P) {[x1 , y1] , … , [xi , yi] , … , [xn , yn] }
• [x1 , y1] [w(P) , min{y : x,y Lmax(P), xw(P) }]
• [xn , yn] [ min{x : x,y Lmax(P), yh(P) } , h(P) ]
ordine canonico
se Rmax(P) {[x,y]}, ovvero P [x,y], e QL è tale che
Lmax(P) Lmax(Q) Q P
L famiglia dei polyomini L-convessi
10
Corrispondenza tra un polyomino
L-convesso P e la famiglia Lmax(P )
Lemma. Sia Rmax(P) 2. QL tale che Lmax(P) Lmax(Q), la dimensione e
la reciproca posizione del primo e secondo (penultimo ed ultimo) rettangolo
massimale di Q coincideranno con quelle del primo e del secondo (penultimo ed
ultimo rispettivamente) rettangolo massimale di P.
Teorema. Sia P L. Se Rmax(P) 3 allora P è univocamente determinato da
Lmax(P) .
Lmax(P) P
nel caso in cui P è costituito
al più da 3 rettangoli
massimali
11
Controesempio
Questo esempio considera due polyomini L-convessi
distinti, con più di tre rettangoli massimali, aventi lo
stesso insieme di L-cammini massimali riportato in
tabella.
11, 2
11,2
11, 2
11,2
10,3
10, 3
10,3 10, 3
9,4
9, 4
9,4
9, 4
8,5
7, 5
7,5
6, 5
7,6
6, 6
6,6
5, 6
6,7
4, 7
5,7
4, 7
3,8
3, 8
3,8
3, 8
2,9
2, 9
2,9
2, 9
4 rettangoli massimali
P1
P2
2 occorrenze
12
Multiset ?
Questo esempio considera due polyomini L-convessi differenti aventi lo stesso
multiset di L-cammini massimali.
P1
P2
non è massimale
multiset
12, 2 , 11,3 , 10, 4 , 8, 6 , 7, 7 , 7, 7 , 6,8 , ... , 2, 12
13
L- cammini bordati
2
3
8
8
8
6
3
6
7
3
3
SE (8,2) , (6,3) , (3,8),
1
E
N
(8,1) , (6,1) , (3,7),
;
Polyomino L-convesso
14
Problemi affrontati
• Consistenza
• Ricostruzione
• Unicità
15
L-cammini bordati
Sia P un polyomino convesso.
Un L-cammino è bordato in P se
parte da una cella del bordo
procede ortogonalmente rispetto allo stesso bordo
quando incontra il bordo opposto gira in senso antiorario
#
quindi procede dritto fino al bordo opposto
In particolare è detto :
• SE bordato se parte dal bordo superione
• EN bordato se parte dal bordo sinistro
• NW bordato se parte dal bordo inferiore
• WS bordato se parte dal bordo destro
#
#
# #
#
#
# #
# 16
Definizione di un L-cammino bordato
Sia
π (i0 , j0 ), , (is , j s ), , (im , jm ) P
un L-cammino che cambia direzione nella cella (is , js ) .
π è detto
• SE bordato se è in direzione SE e (i0 , j0 1), (is , js 1), (im 1, jm ) P
• EN bordato se è in direzione EN e (i0 1, j0 ), (is 1, js ), (im , jm 1) P
• NW bordato se è in direzione NW e (i0 , j0 1), (is , js 1), (im 1, jm ) P
• WS bordato se è in direzione WS e (i0 1, j0 ), (is 1, js ), (im , jm 1) P
(i, j ) denota la cella
i 1, i j 1, j
17
Definizione di una funzione di valutazione
per un L-cammino
La size di un L-cammino
è la funzione
π (i0 , j0 ), , (is , js ), , (im , jm )
s : L( P) N \ {0} N \ {0}
definita da
dove
SE =
s (π ) (h, k )
h im i0 1
k j m j0 1 .
tutte le coppie di interi positivi che rappresentano
la size di un L-cammino SE bordato
Analogamente E N , W S, N W
card (SE) = card (N W) = (P)
card (E N) = card (W S) = h(P)
18
Esempio
SE (1,1) , (8,2), (7,2),(6,3),(8,3),(7,2),(6,1), (3,2), (3,1)
##
#
###
##
#
#
#
#
##
#
##
#
#
#
#
#
E N (1,1) , (3,1), (3,2),(8,1),(8,2),(6,3),(3,6), (3,7)
19
Struttura dell’algoritmo
Prima fase determina gli elementi di Rmax(P)
[x1,y1]
[x2,y2]
x1 x2 x3 x4 and y1 y2 y3 y4
[x3,y3]
[x4,y4]
20
Struttura dell’algoritmo
1,1 1,1
2 , 2
3 , 3 4 , 4
Seconda fase determina la mutua posizione dei
rettangoli massimali a partire…
Ω = (ω1, ω2, ω3, ω4)
ascisse dei SW corners
Σ = (σ1, σ2, σ3, σ4)
ordinate dei SW corners
21
Prima fase
LEMMA. Sia P un polyomino L-convesso. Gli elementi di
Rmax(P) sono univocamente determinati da SE
(o equivalentemente da E N , N W , W S)
yi 1 ki 1
xi
xi 1
Dalla prova di questo lemma si ricava una procedura per
determinare gli elementi di Rmax(P) a partire dall’insieme SE
22
Seconda fase
Due procedure che determinano Ω
OMEGA1 (SE, E
N)
Ω
Ω
OMEGA2 (SE, WS)
incrociato
allineato a sinistra
i 1 i j
i 1 i
allineato a destra
i 1 i xi xi 1
Determinare la reciproca posizione di due rettangoli massimali
significa stabilire quale tipo di intersezione crossing hanno.
23
Seconda fase
Due procedure che determinano Ω
OMEGA1 (SE, E
N)
OMEGA2 (SE, WS)
Ω
Ω
• Scegliendo solo una delle due procedure …
(SE, WS)
(SE, E N)
clock.rotation of π/2
OMEGA2
(SE *, WS *)
… due tipi of sizes sono necessari!!!
Ω
OMEGA2
P
Ω*
2
P*
Σ
24
Seconda fase
OMEGA1 (SE, E
N)
OMEGA2 (SE, WS)
Ω
Ω
Tuttavia il nostro scopo è ricostruire univocamente P da un’unica
coppia di set di sizes.
(SE, E N)
(SE, E N)
clock.rotation of π/2
OMEGA1
(SE *, WS *)
Ω
OMEGA2
Ω*
Σ
Teorema: Ogni polyomino L-convesso è univocamente determinato da (SE, E N).
25