- Notifications
You must be signed in to change notification settings - Fork 1.6k
/
Copy pathMatchstickTriangleSequence.cs
70 lines (68 loc) · 2.3 KB
/
MatchstickTriangleSequence.cs
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
usingSystem.Collections.Generic;
usingSystem.Numerics;
namespaceAlgorithms.Sequences;
/// <summary>
/// <para>
/// Sequence of number of triangles in triangular matchstick arrangement of side n for n>=0.
/// </para>
/// <para>
/// M. E. Larsen, The eternal triangle – a history of a counting problem, College Math. J., 20 (1989), 370-392.
/// https://web.math.ku.dk/~mel/mel.pdf.
/// </para>
/// <para>
/// OEIS: http://oeis.org/A002717.
/// </para>
/// </summary>
publicclassMatchstickTriangleSequence:ISequence
{
/// <summary>
/// <para>
/// Gets number of triangles contained in an triangular arrangement of matchsticks of side length n.
/// </para>
/// <para>
/// This also counts the subset of smaller triangles contained within the arrangement.
/// </para>
/// <para>
/// Based on the PDF referenced above, the sequence is derived from step 8, using the resulting equation
/// of f(n) = (n(n+2)(2n+1) -(delta)(n)) / 8. Using BigInteger values, we can effectively remove
/// (delta)(n) from the previous by using integer division instead.
/// </para>
/// <para>
/// Examples follow.
/// <pre>
/// .
/// / \ This contains 1 triangle of size 1.
/// .---.
///
/// .
/// / \ This contains 4 triangles of size 1.
/// .---. This contains 1 triangle of size 2.
/// / \ / \ This contains 5 triangles total.
/// .---.---.
///
/// .
/// / \ This contains 9 triangles of size 1.
/// .---. This contains 3 triangles of size 2.
/// / \ / \ This contains 1 triangles of size 3.
/// .---.---.
/// / \ / \ / \ This contains 13 triangles total.
/// .---.---.---.
/// </pre>
/// </para>
/// </summary>
publicIEnumerable<BigInteger>Sequence
{
get
{
varindex=BigInteger.Zero;
vareight=newBigInteger(8);
while(true)
{
vartemp=index*(index+2)*(index*2+1);
varresult=BigInteger.Divide(temp,eight);
yieldreturnresult;
index++;
}
}
}
}