-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathpart2.pl
87 lines (74 loc) · 2.01 KB
/
part2.pl
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
#!/usr/bin/perl
@velocities = (
[0, 0, 0],
[0, 0, 0],
[0, 0, 0],
[0, 0, 0],
);
@positions = ();
foreach my $line (<>) {
my ($x, $y, $z) = ($line =~ /<x=(-?\d+), y=(-?\d+), z=(-?\d+)>/);
push(@positions, [$x, $y, $z]);
}
@initial_velocities = map { [@$_] } @velocities;
@initial_positions = map { [@$_] } @positions;
@cycle_length = (-1, -1, -1);
for (my $n = 0; $n < 100000000; ++$n) {
for (my $i = 0; $i < 4; ++$i) {
for (my $j = 0; $j < 3; ++$j) {
my $delta = 0;
for (my $k = 0; $k < 4; ++$k) {
if ($positions[$i][$j] < $positions[$k][$j]) {
++$delta;
} elsif ($positions[$i][$j] > $positions[$k][$j]) {
--$delta;
}
}
$velocities[$i][$j] += $delta;
}
}
for (my $i = 0; $i < 4; ++$i) {
for (my $j = 0; $j < 3; ++$j) {
$positions[$i][$j] += $velocities[$i][$j];
}
}
my $checked = 0;
for (my $k = 0; $k < 3; ++$k) {
if ($cycle_length[$k] == -1) {
$checked = 1;
my $ok = 1;
for (my $i = 0; $i < 4; ++$i) {
if ($positions[$i][$k] != $initial_positions[$i][$k]) {
$ok = 0;
} else {
}
if ($velocities[$i][$k] != $initial_velocities[$i][$k]) {
$ok = 0;
}
}
if ($ok) {
$cycle_length[$k] = $n;
}
}
}
if ($n % 10000 == 0) {
print $n, " ", $cycle_length[0], " ", $cycle_length[1] , " ", $cycle_length[2], "\n";
}
last if $checked == 0;
}
for (my $i = 0; $i < 3; ++$i) {
print ++$cycle_length[$i], " "
}
print "\n";
sub gcd {
my ($a, $b) = @_;
return $b if ($a == 0);
return $a if ($b == 0);
return gcd($b, $a % $b);
}
$ret = 1;
for (my $i = 0; $i < 3; ++$i) {
my $g = gcd($ret, $cycle_length[$i]);
$ret = $ret / $g * $cycle_length[$i];
}
print "$ret\n";