Example #1
0
        public ContinuedFraction(int n)
        {
            int root = (int)Math.Sqrt(n);
            if(root * root == n) {
                throw new ArgumentException("N can't be a perfect square");
            }

            continuedFraction = new List<long>();

            HashSet<String> store = new HashSet<String>();
            RadicalFraction current = new RadicalFraction(n, 0, 1, 1, 0);
            while(!store.Contains(current.GetKey())) {
                store.Add(current.GetKey());
                long integerPart = current.GetIntegerPart();
                continuedFraction.Add(integerPart);

                current = current.MinusFraction(new RadicalFraction(n, integerPart, 0, 1, 0));
                current = current.GetReciprocal();
                current.Simplify();
            }

            h = new BigInteger[] { 0, 1 };
            k = new BigInteger[] { 1, 0 };

            index = 0;
        }
Example #2
0
        public RadicalFraction MinusFraction(RadicalFraction other)
        {
            if (other.radical != this.radical) {
                throw new ArgumentException("Radicals differ");
            }

            if (other.denominator[0] != this.denominator[0] || other.denominator[1] != this.denominator[1]) {
                other = other.Multiply(new RadicalFraction(radical, denominator[0], denominator[1], denominator[0], denominator[1]));
            }

            long[] newNumerator = new long[2];
            newNumerator[0] = this.numerator[0] - other.numerator[0];
            newNumerator[1] = this.numerator[1] - other.numerator[1];

            long[] newDenominator = new long[2];
            newDenominator[0] = this.denominator[0];
            newDenominator[1] = this.denominator[1];

            return new RadicalFraction(radical, newNumerator, newDenominator);
        }
Example #3
0
        public RadicalFraction Multiply(RadicalFraction other)
        {
            if (other.radical != this.radical) {
                throw new ArgumentException("Radicals differ");
            }

            long[] newNumerator = new long[2];
            newNumerator[0] = this.numerator[0] * other.numerator[0] + this.numerator[1] * other.numerator[1] * radical;
            newNumerator[1] = this.numerator[0] * other.numerator[1] + this.numerator[1] * other.numerator[0];

            long[] newDenominator = new long[2];
            newDenominator[0] = this.denominator[0] * other.denominator[0] + this.denominator[1] * other.denominator[1] * radical;
            newDenominator[1] = this.denominator[0] * other.denominator[1] + this.denominator[1] * other.denominator[0];

            RadicalFraction result = new RadicalFraction(radical, newNumerator, newDenominator);

            return result;
        }
Example #4
0
        public RadicalFraction Divide(RadicalFraction other)
        {
            if (other.radical != this.radical) {
                throw new ArgumentException("Radicals differ");
            }

            return this.Multiply(other.GetReciprocal());
        }