Estimate the number of discrepancy 2 sequences - Revision history

Guy Srinivasan at 21:20, 20 January 2010

2010-01-20T21:20:59Z

← Older revision		Revision as of 14:20, 20 January 2010
Line 10:		Line 10:
	class Program		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 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, 2245417744, 3449035210, 5401519132, 5272055840, 4470731784, 7013199284, 6544305958, 10138992314, 15284191798, 9839916650, 8579138578, 13353765596, 12170187396, 18839898968, 17144676512, 26322363104, 39457407576 };
	static void Run()		static void Run()
	{		{
Line 17:		Line 17:
	double log = 0;		double log = 0;
	double epsilon = 0.000001;		double epsilon = 0.000001;
	for (int idx = 1; idx <= ~~1125~~; idx++)		for (int idx = 1; idx <= answers.Length - 1; idx++)
	{		{
	int HAPs = 0;		int HAPs = 0;

Guy Srinivasan at 21:17, 20 January 2010

2010-01-20T21:17:28Z

Show changes

Guy Srinivasan: New page:

 using System; using System.Collections.Generic; using System.Linq; using System.Text; using System.IO;  namespace EDP {     class Program     {         static double[] answers = { 1,...

2010-01-20T03:27:10Z

New page: <pre> using System; using System.Collections.Generic; using System.Linq; using System.Text; using System.IO; namespace EDP { class Program { static double[] answers = { 1,...

New page

<pre>
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;
}
}
}
</pre>