View difference between Paste ID: 0xpNJLjJ and sJqyXHBy

-f[1][1]=1.000000, f[1][1]=1.000000, f[2][1]=0.500000, f[1][2]=0.500000
+#include<iostream>
-f[1][2]=0.500000, f[2][1]=0.500000, f[2][2]=1.250000, f[2][2]=1.250000
+#include<stdio.h>
-f[1][3]=0.333333, f[3][1]=0.333333, f[2][3]=2.833333, f[3][2]=2.833333
+#include<math.h>
-f[1][4]=0.250000, f[4][1]=0.250000, f[2][4]=5.125000, f[4][2]=5.125000
+using namespace std;
-f[1][5]=0.200000, f[5][1]=0.200000, f[2][5]=8.100000, f[5][2]=8.100000
+
-f[1][6]=0.166667, f[6][1]=0.166667, f[2][6]=11.750000, f[6][2]=11.750000
+unsigned long long int gcd(unsigned long long int a, unsigned long long int b)
-f[1][7]=0.142857, f[7][1]=0.142857, f[2][7]=16.071429, f[7][2]=16.071429
+{
-f[1][8]=0.125000, f[8][1]=0.125000, f[2][8]=21.062500, f[8][2]=21.062500
+     unsigned long long int c;
-f[1][9]=0.111111, f[9][1]=0.111111, f[2][9]=26.722222, f[9][2]=26.722222
+     while ( b != 0 )
-f[1][10]=0.100000, f[10][1]=0.100000, f[2][10]=33.050000, f[10][2]=33.050000
+        {
-f[1][11]=0.090909, f[11][1]=0.090909, f[2][11]=40.045455, f[11][2]=40.045455
+            c = b;
-f[1][12]=0.083333, f[12][1]=0.083333, f[2][12]=47.708333, f[12][2]=47.708333
+            b = a % b;
-f[1][13]=0.076923, f[13][1]=0.076923, f[2][13]=56.038462, f[13][2]=56.038462
+            a = c;
-f[1][14]=0.071429, f[14][1]=0.071429, f[2][14]=65.035714, f[14][2]=65.035714
+        }
-f[1][15]=0.066667, f[15][1]=0.066667, f[2][15]=74.700000, f[15][2]=74.700000
+
-f[1][16]=0.062500, f[16][1]=0.062500, f[2][16]=85.031250, f[16][2]=85.031250
+        return a;
-f[1][17]=0.058824, f[17][1]=0.058824, f[2][17]=96.029412, f[17][2]=96.029412
+}
-f[1][18]=0.055556, f[18][1]=0.055556, f[2][18]=107.694444, f[18][2]=107.694444
+
-f[1][19]=0.052632, f[19][1]=0.052632, f[2][19]=120.026316, f[19][2]=120.026316
+struct fraction
-f[1][20]=0.050000, f[20][1]=0.050000, f[2][20]=133.025000, f[20][2]=133.025000
+{
-f[1][21]=0.047619, f[21][1]=0.047619, f[2][21]=146.690476, f[21][2]=146.690476
+       unsigned long long int num,denom;
-f[1][22]=0.045455, f[22][1]=0.045455, f[2][22]=161.022727, f[22][2]=161.022727
+};
-f[1][23]=0.043478, f[23][1]=0.043478, f[2][23]=176.021739, f[23][2]=176.021739
+
-f[1][24]=0.041667, f[24][1]=0.041667, f[2][24]=191.687500, f[24][2]=191.687500
+void reducelowest(struct fraction& a)
-f[1][25]=0.040000, f[25][1]=0.040000, f[2][25]=208.020000, f[25][2]=208.020000
+{
-f[1][26]=0.038462, f[26][1]=0.038462, f[2][26]=225.019231, f[26][2]=225.019231
+     unsigned long long int m;
-f[1][27]=0.037037, f[27][1]=0.037037, f[2][27]=242.685185, f[27][2]=242.685185
+     m=gcd(a.num,a.denom);
-f[1][28]=0.035714, f[28][1]=0.035714, f[2][28]=261.017857, f[28][2]=261.017857
+     /*unsigned long long int i;
-f[1][29]=0.034483, f[29][1]=0.034483, f[2][29]=280.017241, f[29][2]=280.017241
+     if(a.num>a.denom)
-f[1][30]=0.033333, f[30][1]=0.033333, f[2][30]=299.683333, f[30][2]=299.683333
+         m=a.denom;
-f[1][31]=0.032258, f[31][1]=0.032258, f[2][31]=320.016129, f[31][2]=320.016129
+     else
-f[1][32]=0.031250, f[32][1]=0.031250, f[2][32]=341.015625, f[32][2]=341.015625
+         m=a.num;
-f[1][33]=0.030303, f[33][1]=0.030303, f[2][33]=362.681818, f[33][2]=362.681818
+
-f[1][34]=0.029412, f[34][1]=0.029412, f[2][34]=385.014706, f[34][2]=385.014706
+     for(i=m;i>=1;i--)
-f[1][35]=0.028571, f[35][1]=0.028571, f[2][35]=408.014286, f[35][2]=408.014286
+     {
-f[1][36]=0.027778, f[36][1]=0.027778, f[2][36]=431.680556, f[36][2]=431.680556
+         if(a.num%i==0&&a.denom%i==0)
-f[1][37]=0.027027, f[37][1]=0.027027, f[2][37]=456.013514, f[37][2]=456.013514
+             break;
-f[1][38]=0.026316, f[38][1]=0.026316, f[2][38]=481.013158, f[38][2]=481.013158
+      }*/
-f[1][39]=0.025641, f[39][1]=0.025641, f[2][39]=506.679487, f[39][2]=506.679487
+     a.num=a.num/m;
-f[1][40]=0.025000, f[40][1]=0.025000, f[2][40]=533.012500, f[40][2]=533.012500
+     a.denom=a.denom/m;
-f[1][41]=0.024390, f[41][1]=0.024390, f[2][41]=560.012195, f[41][2]=560.012195
+ }
-f[1][42]=0.023810, f[42][1]=0.023810, f[2][42]=587.678571, f[42][2]=587.678571
+
-f[1][43]=0.023256, f[43][1]=0.023256, f[2][43]=616.011628, f[43][2]=616.011628
+int main()
-f[1][44]=0.022727, f[44][1]=0.022727, f[2][44]=645.011364, f[44][2]=645.011364
+{
-f[1][45]=0.022222, f[45][1]=0.022222, f[2][45]=674.677778, f[45][2]=674.677778
+    FILE*fp=fopen("omoout.txt","w");
-f[1][46]=0.021739, f[46][1]=0.021739, f[2][46]=705.010870, f[46][2]=705.010870
+    struct fraction f[104][104];
-f[1][47]=0.021277, f[47][1]=0.021277, f[2][47]=736.010638, f[47][2]=736.010638
+    for(int i=1;i<=51;i++)
-f[1][48]=0.020833, f[48][1]=0.020833, f[2][48]=767.677083, f[48][2]=767.677083
+    {
-f[1][49]=0.020408, f[49][1]=0.020408, f[2][49]=800.010204, f[49][2]=800.010204
+            f[1][i].num=1; f[1][i].denom=i;
-f[1][50]=0.020000, f[50][1]=0.020000, f[2][50]=833.010000, f[50][2]=833.010000
+            f[i][1]=f[1][i];
-f[1][51]=0.019608, f[51][1]=0.019608, f[2][51]=866.676471, f[51][2]=866.676471
+            f[2][i].num=(2*i*i*i-2*i+3); f[2][i].denom=6*i;
-f[1][52]=0.019231, f[52][1]=0.019231, f[2][52]=901.009615, f[52][2]=901.009615
+            reducelowest(f[2][i]);
-f[1][53]=0.018868, f[53][1]=0.018868, f[2][53]=936.009434, f[53][2]=936.009434
+            f[i][2].num=f[2][i].num;f[i][2].denom=f[2][i].denom;
-f[1][54]=0.018519, f[54][1]=0.018519, f[2][54]=971.675926, f[54][2]=971.675926
+
-f[1][55]=0.018182, f[55][1]=0.018182, f[2][55]=1008.009091, f[55][2]=1008.009091
+     }
-f[1][56]=0.017857, f[56][1]=0.017857, f[2][56]=1045.008929, f[56][2]=1045.008929
+
-f[1][57]=0.017544, f[57][1]=0.017544, f[2][57]=1082.675439, f[57][2]=1082.675439
+
-f[1][58]=0.017241, f[58][1]=0.017241, f[2][58]=1121.008621, f[58][2]=1121.008621
+    for(int i=3;i<=50;i++)
-f[1][59]=0.016949, f[59][1]=0.016949, f[2][59]=1160.008475, f[59][2]=1160.008475
+    {
-f[1][60]=0.016667, f[60][1]=0.016667, f[2][60]=1199.675000, f[60][2]=1199.675000
+            for(int j=2;j<=49;j++)
-f[1][61]=0.016393, f[61][1]=0.016393, f[2][61]=1240.008197, f[61][2]=1240.008197
+            {
-f[1][62]=0.016129, f[62][1]=0.016129, f[2][62]=1281.008065, f[62][2]=1281.008065
+                    if(f[i][j+1].denom==0)
-f[1][63]=0.015873, f[63][1]=0.015873, f[2][63]=1322.674603, f[63][2]=1322.674603
+                    {
-f[1][64]=0.015625, f[64][1]=0.015625, f[2][64]=1365.007813, f[64][2]=1365.007813
+                      f[i][j+1].num=((f[i-1][j].denom)*((f[i-1][j+1].denom)*(f[i][j].denom)+(f[i-1][j+1].num)*(f[i][j].num)));
-f[1][65]=0.015385, f[65][1]=0.015385, f[2][65]=1408.007692, f[65][2]=1408.007692
+                      f[i][j+1].denom=(f[i-1][j].num)*(f[i-1][j+1].denom)*(f[i][j].denom);
-f[1][66]=0.015152, f[66][1]=0.015152, f[2][66]=1451.674242, f[66][2]=1451.674242
+                      reducelowest(f[i][j+1]);
-f[1][67]=0.014925, f[67][1]=0.014925, f[2][67]=1496.007463, f[67][2]=1496.007463
+                    }
-f[1][68]=0.014706, f[68][1]=0.014706, f[2][68]=1541.007353, f[68][2]=1541.007353
+            }
-f[1][69]=0.014493, f[69][1]=0.014493, f[2][69]=1586.673913, f[69][2]=1586.673913
+     }
-f[1][70]=0.014286, f[70][1]=0.014286, f[2][70]=1633.007143, f[70][2]=1633.007143
+    for(int i=1;i<=50;i++)
-f[1][71]=0.014085, f[71][1]=0.014085, f[2][71]=1680.007042, f[71][2]=1680.007042
+    {
-f[1][72]=0.013889, f[72][1]=0.013889, f[2][72]=1727.673611, f[72][2]=1727.673611
+            //for(int j=1;j<=50;j++)
-f[1][73]=0.013699, f[73][1]=0.013699, f[2][73]=1776.006849, f[73][2]=1776.006849
+            //{
-f[1][74]=0.013514, f[74][1]=0.013514, f[2][74]=1825.006757, f[74][2]=1825.006757
+                    fprintf(fp,"f(%d,%d)=%lld/%lld    %f\n",i,i,f[i][i].num,f[i][i].denom,pow((f[i][i].num-1),0.5));
-f[1][75]=0.013333, f[75][1]=0.013333, f[2][75]=1874.673333, f[75][2]=1874.673333
+            //}
-f[1][76]=0.013158, f[76][1]=0.013158, f[2][76]=1925.006579, f[76][2]=1925.006579
+     }
-f[1][77]=0.012987, f[77][1]=0.012987, f[2][77]=1976.006494, f[77][2]=1976.006494
+     fclose(fp);
-f[1][78]=0.012821, f[78][1]=0.012821, f[2][78]=2027.673077, f[78][2]=2027.673077
+
-f[1][79]=0.012658, f[79][1]=0.012658, f[2][79]=2080.006329, f[79][2]=2080.006329
+     return 0;
-f[1][80]=0.012500, f[80][1]=0.012500, f[2][80]=2133.006250, f[80][2]=2133.006250
+}