Форум StopLinux

Что то компилер с ума спрыгивает . Тупой вопрос, 7^21 это скока? Просто у двух компиляторов си мнения капитально расходятся с виндовым калькулятором, а на глаз прав именно калькулятор . Разница в младших разрядах, а как раз младшие то разряды хоть как нужны .

И виндовый кальк так же говорит. А Си почему то уверен что 558545864083284030, несмотря на long double
В результате у моего быдлокода нет шансов

#include <stdio.h>
#include <string.h>
#include <conio.h>
#include <math.h>


long double A,B,C,D;
long double F_n(long double n)
{
long double f=0;
                  //D=3*(A-C);
                 while (n<=B)
                        {
                        n=n+A;
                        f=f+D;
                        }
                 f=f+(n-C);
return f;
}

int main(void)
{
 long double f,n=0,n1=0,s=0,s1=0;//A,B,C;
  // A=pow(21,7);
  // B=pow(7,21);
  // C=pow(12,7);

   A=50;
   B=2000;
   C=40;
   D=3*(A-C);
    
    //f=F_n(n); 

    // while (n<=B)
    //  {
    //    s=s+F_n(n);
    //     n++;
    //  }
      while (n1<B)
        {
        n1=n1+A;
        s1=s1+A*((F_n(n1)-A)+((A+1)/2));
        }s1=s1+F_n(0);

   //printf(" \n %Lf %Lf %Lf",A,B,C);
   //printf(" \n %Lf %Lf %Lf",f,s,s1);
      printf(" \n %Lf",s1);
   getch();
   return 0;
}

Хотя у него и так нет шансов, явно нужен другой подход

Да пожалуйста, хоть 7^221, на
5843416508321765137412223462176261055753559544095817424625430781485355014524409616639950700057743469282269362223748969437289370453147445512994056218802090543396973006957507863824356124007

Калькуль то как раз выдает нормально. Просто для таких чисел встроеных в С типов данных маловато будет . У long double 15 знаков точности. Ну и учитывая что при степени B = 17 считается 30 секунд для 21 это  около 3 дней . Явно не то решение.

Это не решение. Время вычисления получится несколько дней, а нужно только несколько младших разрядов. Должна быть какая то бяка позволяющая это сделать проще

Маладэц, а теперь посчитай сколько итераций там. Можешь воспользоваться калькулятором.

Не, видно не совсем внимательно читал.
Итак,
a=1
ищем первое простое равное 2.
ищем максимальную степень этого простого не превышающую 20, т.е 4
a=a*2^4
ищем второе простое, т.е 3.
ищем максимальную степень 3 не превышающую 20, нашли.
a=a*3^2
следующее простое - 5, 5^2>20 поэтому устанавливаем флаг отключающий возведение простого в степень и сравнение с 20
a=a*5
a=a*7
Итд.

Раскопаем труп тему, но это надо запостить!

Haskell:

main = do print $ foldl lcm 1 [1..20]

$ time ./hello.exe
232792560

real    0m0.069s
user    0m0.031s
sys     0m0.015s