Метод Лагранжа

@Andy_Coldfield · Регистрация: 25.11.2017

Студворк — интернет-сервис помощи студентам

Помогите, не знаю как подступиться к этой задаче. Суть следующая, пользователь вводит произвольный многочлен. И вот его надо привести к канонической форме и проверить ранг с помощью метода Лагранжа.
Вот к примеру есть следующий многочлен:
f(x₁

Добавлено через 6 минут
f(x₁,x₂,x₃) = -x²₁ - 5x²₂ - 7x²₃ + 2x₁x₂-2x₁x₃+6x₂x₃

@Andy_Coldfield · 09.02.2020, 01:19 **[ТС]**

assemberist, не поможете разобраться с двумя кусками кода...

Кликните здесь для просмотра всего текста

C
1
2
3
4
5
6
7
8
9
10
11
12
13
int gf_single_divide(int a, int b)
{
  int sum_j;
 
  gf_modar_setup();
  if (b == 0) return -1;
  if (a == 0) return 0;
  sum_j = B_TO_J[a] - B_TO_J[b];
#ifdef W_16
  if (sum_j < 0) sum_j += Modar_nwm1;
#endif
  return (int) J_TO_B[sum_j];
}

Кликните здесь для просмотра всего текста

C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
void gf_modar_setup()
{
  int j, b, t;
  if (gf_already_setup) return;
 
  B_TO_J = (int *) malloc(sizeof(int)*Modar_nw);
  if (B_TO_J == NULL) {
    perror("gf_modar_setup, malloc B_TO_J");
    exit(1);
  }
   /* When the word size is 8 bits, make three copies of the table so that
      you don't have to do the extra addition or subtraction in the
      multiplication/division routines */
 
#ifdef W_8
  J_TO_B = (int *) malloc(sizeof(int)*Modar_nw*3);
#elif W_16
  J_TO_B = (int *) malloc(sizeof(int)*Modar_nw);
#endif
  if (J_TO_B == NULL) {
    perror("gf_modar_setup, malloc J_TO_B");
    exit(1);
  }
  for (j = 0; j < Modar_nw; j++) {
    B_TO_J[j] = Modar_nwm1;
    J_TO_B[j] = 0;
  } 
 
  b = 1;
  for (j = 0; j < Modar_nwm1; j++) {
    if (B_TO_J[b] != Modar_nwm1) {
      fprintf(stderr, "Error: j=%d, b=%d, B->J[b]=%d, J->B[j]=%d (0%o)\n",
              j, b, B_TO_J[b], J_TO_B[j], (b << 1) ^ Modar_poly);
      exit(1);
    }
    B_TO_J[b] = j;
    J_TO_B[j] = b;
    b = b << 1;
    if (b & Modar_nw) b = (b ^ Modar_poly) & Modar_nwm1;
  }
/*   for (j = 0; j < Modar_nw; j++) { */
/*     printf("%3d b->j %3d        %3d j->b %3d\n",  */
/*             j, B_TO_J[j], j, J_TO_B[j]; */
/*   } */
 
#ifdef W_8
    for (j = 0; j < Modar_nwm1; j++) {
      J_TO_B[j+Modar_nwm1] = J_TO_B[j];
      J_TO_B[j+2*Modar_nwm1] = J_TO_B[j];
    }
    J_TO_B += Modar_nwm1;
#endif
  gf_already_setup = 1;
 
}

Добавлено через 28 секунд
assemberist, Мне второй не очень понятен....

@assemberist · 09.02.2020, 09:29

Сообщение от Andy_Coldfield

Мне второй не очень понятен

Что uменно не понятно?
Что-то uз сuнтаксuса?
Потому что я не могу сказать как он работает, пока не увuжу заченuя u тuпы всех переменных.

@Andy_Coldfield · 09.02.2020, 13:49 **[ТС]**

assemberist,

Кликните здесь для просмотра всего текста

C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
#define prim_poly_32 020000007
#define prim_poly_16 0210013
#define prim_poly_8 0435
 
static int gf_already_setup = 0;
 
#ifdef W_8
  static int Modar_w = 8;
  static int Modar_nw = 256;
  static int Modar_nwm1 = 255;
  static int Modar_poly = prim_poly_8;
#elif W_16
  static int Modar_w = 16;
  static int Modar_nw = 65536;
  static int Modar_nwm1 = 65535;
  static int Modar_poly = prim_poly_16;
#endif
 
static int *B_TO_J;
static int *J_TO_B;

Добавлено через 19 минут
assemberist, меня больше смущают обороты ifdef ... elsedef..

@assemberist · 09.02.2020, 16:02

Сообщение от Andy_Coldfield

ifdef ... elsedef..

Это дuректuвы препроцессора. Онu позволяют пuсать макроссы, компuлuровать по условuю u делать много другuх полезных вещей. В твоем случае онu определяют размерность данных: 1 uлu 2 байта.
Правда тогда нужно лuбо пропuсать #define для нужной опцuu лuбо указать компuлятору, чтобы он это сделал сам прu компuляцuu.

@Avazart · 11.02.2020, 20:30

Сообщение от Andy_Coldfield

меня больше смущают обороты ifdef ... elsedef..

И правильно что смущают, нужно к минимуму сводить.

@Andy_Coldfield · 11.02.2020, 23:54 **[ТС]**

Avazart, что значит сводить к минимуму?

@Avazart · 12.02.2020, 01:29

Очевидно что

C++
1
2
3
  static int Modar_nw = 256;
  static int Modar_nwm1 = 255;
  static int Modar_poly = prim_poly_8

Зависят от Modar_w т.е. стоило бы выразить одно через другое и ifdef ... elsedef был бы меньшим или вообще можно было бы викинуть.
Так же вероятно можно было бы подумать о применении шаблонов.

@Andy_Coldfield · 26.02.2020, 17:11 **[ТС]**

assemberist, а не подскажите алгоритм как привести матрицу к ступенчатому виду...

@assemberist · 26.02.2020, 23:02

Сообщение от Andy_Coldfield

а не подскажите алгоритм как привести матрицу к ступенчатому виду

Лень ввестu "привести матрицу к ступенчатому виду с++" в яндекс, зато не лень доставать assemberistа.
http://www.intuit.ru/studies/c... cture/1978

Новые блоги и статьи Все статьи Все блоги /
сукцессия 29. Переход от одних деревьев на другие делать более или менее вероятностным? anaschu 12.07.2026 Насколько смена типов микоризы — исключительное событие в двухвековой сукцессии? Оценка вероятности в пространстве параметров В текущей версии модели успешно реализован ключевой механизм. . .	сукцессия 27. Думаю, как переделывать уже написанную статью с планами на сукцессию. anaschu 12.07.2026 Анализ соответствия модели требованиям Реализованные компоненты: Механизм закисления почвы через протонную помпу Конкуренция между типами микориз pH как триггер сукцессии C/ P соотношение. . .	Сукцессия 26. Мат модель создана. anaschu 12.07.2026 Модель смены растительных сукцессий посредством управления грибами работает внутри небольшой ячейки почвы, восстанавливающейся после пожара, где ненадолго бывшее царство хвойных снова захватили. . .	Решил проблему с ошибкой пагинации сообщений с сервера на алгоритме обхода дерева "Эстафета хвоста". Hrethgir 12.07.2026 Проблема была в том, что удалялась именно новая кнопка, а не старая. Ни один ИИ не обнаружил это, а сам я смог только когда с работой стало попроще и когда заставил работать будущее автономное. . .
сукцессия 25. Хронология ошибок anaschu 12.07.2026 # От 50-тонного гриба до устойчивого леса: хроника ошибок при построении модели вековой сукцессии микоризы ## О чём эта статья В процессе построения ОДУ-модели (система дифференциальных. . .	сукцессия 24. Промежуточное общее описание модели anaschu 12.07.2026 Хендофф: модель АМ→ЭКМ сукцессии микоризы (ризосфера, 50 лет) Содержание проекта Симуляция вековой (50 лет) экологической сукцессии в почве леса Основные участники: АМ-гриб, ЭКМ-гриб,. . .	сукцессия 23. Более физиологичная физиология, более экологичная экология, более диффурные диффуры. anaschu 12.07.2026 Что реально нашли и починили за эти 5 часов Правило Линдемана (КПД конверсии сахара в тело, kEff) — раньше 100% полученного углерода шло прямо в биомассу гриба; теперь только kEff=0. 5 (после. . .	сукцессия 22. От артефактов к физиологии: калибровка агентной модели грибной сукцессии для воспроизведения сезонной динамики и pH-плато anaschu 11.07.2026 Аннотация В данной работе представлена калибровка агентной модели динамики грибных сообществ (fungal-succession), направленная на устранение нефизичных артефактов (коллапс биомассы, мгновенное. . .

@Andy_Coldfield 21 / 19 / 6 Регистрация: 25.11.2017 Сообщений: 708

	Метод Лагранжа 04.02.2020, 17:02. Показов 3844. Ответов 28 Метки нет (Все метки) Помогите, не знаю как подступиться к этой задаче. Суть следующая, пользователь вводит произвольный многочлен. И вот его надо привести к канонической форме и проверить ранг с помощью метода Лагранжа. Вот к примеру есть следующий многочлен: f(x₁ Добавлено через 6 минут f(x₁,x₂,x₃) = -x²₁ - 5x²₂ - 7x²₃ + 2x₁x₂-2x₁x₃+6x₂x₃ 0

@Andy_Coldfield 21 / 19 / 6 Регистрация: 25.11.2017 Сообщений: 708
	11.02.2020, 23:54 [ТС]
	Avazart, что значит сводить к минимуму? 0

@Andy_Coldfield 21 / 19 / 6 Регистрация: 25.11.2017 Сообщений: 708
	26.02.2020, 17:11 [ТС]
	assemberist, а не подскажите алгоритм как привести матрицу к ступенчатому виду... 0