Estimate the number of discrepancy 2 sequences

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.IO;

namespace EDP
{
    class Program
    {
        static double[] answers = { 1, 2, 4, 6, 12, 18, 28, 44, 88, 100, 152, 240, 370, 556, 882, 750, 1500, 2250, 2784, 4284, 6438, 6062, 9526, 14856, 22944, 26164, 39528, 35122, 54800, 80940, 81326, 122422, 244844, 234934, 356154, 309068, 388042, 589796, 900000, 813466, 1212450, 1837030, 1882194, 2921946, 4544342, 2274560, 3542738, 5495686, 8436986, 9597362, 11352364, 10876536, 16040144, 23626898, 24060696, 22332908, 34665316, 32813078, 48774494, 77333978, 79086932, 118815026, 181723488, 101003862, 202007724, 171112060, 175332604, 266556200, 388397590, 379477372, 355489350, 544231036, 685850834, 1035438526, 1600139990, 967443184, 1435672238, 1335563798, 1305891128, 1985108178, 2984972718, 3449035210, 2245417744, 5401519132, 5272055840, 4470731784, 7013199284, 6544305958, 10138992314, 15284191798, 9839916650, 8579138578, 13353765596, 12170187396, 18839898968 };
        static void Run()
        {
            double logLower = 0;
            double logUpper = 0;
            double epsilon = 0.000001;
            for (int idx = 1; idx <= 94; idx++)
            {
                int HAPs = 0;
                for (int i = 1; i <= idx; i++)
                {
                    if (i != idx && idx % i == 0 && (((idx / i) - 1) % 2) == 0) HAPs++;
                }
                double p1 = Math.Pow(0.75, HAPs) - Math.Pow(0.5, HAPs);
                double p2 = Math.Pow(0.5, HAPs);
                double multiplier = p1 + p1 + 2 * p2;
                double lowerEV = Math.Log(2 * Math.Floor(Math.Exp(logLower + Math.Log(multiplier)) / 2 + epsilon));
                double upperEV = Math.Log(2 * Math.Ceiling(Math.Exp(logUpper + Math.Log(multiplier)) / 2 - epsilon));
                logLower = lowerEV;
                logUpper = upperEV;
                Console.WriteLine(idx + ":\t" + HAPs + "\t" + Math.Exp(lowerEV) + "\t" + Math.Exp(upperEV) + "\t" + answers[idx] / Math.Exp(upperEV));
            }
            Console.ReadLine();
            return;
        }
        static void Main(string[] args)
        {
            Run();
            return;
        }
    }
}