Генерация ПСП

@S.l.a.y.e.R · Регистрация: 02.05.2009

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

Пожалуста, помогите. Задача состоит в том чтобы написать на С++ генератор псевдослучайных последовательностей по формуле:

Y[i+1]=F*Y[i][/B],

где F - матрица(любая невырожденная), Y - вектора.
Все вычисления должны быть по модулю 5 (элементы матрицы и векторов соответственно тоже - 0 1 2 3 4 ).
Нужно написать программу, которая реализует это умножение, причем циклическое, т.е. вектор-результат еще раз умножится на ту же матрицу и получится следующий вектор и т. д. Также нужно чтобы вычислялся период т.е. когда вектора-результаты начнут повторяться. И еще прорамма должна позволять выводить результат задаваемого умножения т.е. например результат после 15-го умножения (1 результат) или результат с 10-го до 16-го умножения (7 результатов).
Я пока только написал достаточно дилетантскую программу, которая только считает результат 1-го умножения. Вот она:

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
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
#include<stdio.h>
#include<stdlib.h>
#include<conio.h>
 
int main(void) {
int c1;
int i, j, k;
randomize();
printf("vvedite rang matrici(2, 3 ili 4) - n: \n");
scanf("%d",&c1);
    switch(c1) {
 
case 2:
int S2=2;
int A2[2][2];
int B2[2][1];
int C2[2][1];
for(j=0; j<S2; j++){
for(k=0; k<1; k++){
B2[j][k]=random(5);}}
 
for(i=0; i<S2; i++){
for(j=0; j<S2; j++){
A2[i][j]=random(5);}}
printf("matica F: \n");
for(i=0; i<S2; i++){
for(j=0; j<S2; j++){
printf("%3d ",A2[i][j]);}
   printf("\n");
   printf("\n");}
 
    printf("\n");
printf("vector Y[i]: \n");
for(j=0; j<S2; j++){
for(k=0; k<1; k++){
printf("%3d ",B2[j][k]);}
   printf("\n");
   printf("\n");}
 
   printf("\n");
 
for(i=0; i<S2; i++){
for(k=0; k<1; k++){
C2[i][k]=0;
for(j=0; j<S2; j++){
C2[i][k]=(C2[i][j]+(A2[i][j]*B2[j][k])%5)%5;}}}
printf("vector Y[i+1]=F*Y[i]: \n");
for(i=0; i<S2; i++){
for(k=0; k<1; k++){
printf("%3d ",C2[i][k]);}
printf("\n");
printf("\n");}
break;
 
case 3:
int S3=3;
int A3[3][3];
int B3[3][1];
int C3[3][1];
for(j=0; j<S3; j++){
for(k=0; k<1; k++){
B3[j][k]=random(5);}}
 
for(i=0; i<S3; i++){
for(j=0; j<S3; j++){
A3[i][j]=random(5);}}
printf("matica F: \n");
for(i=0; i<S3; i++){
for(j=0; j<S3; j++){
printf("%3d ",A3[i][j]);}
   printf("\n");
   printf("\n");}
 
    printf("\n");
printf("vector Y[i]: \n");
for(j=0; j<S3; j++){
for(k=0; k<1; k++){
printf("%3d ",B3[j][k]);}
   printf("\n");
   printf("\n");}
 
   printf("\n");
 
for(i=0; i<S3; i++){
for(k=0; k<1; k++){
C3[i][k]=0;
for(j=0; j<S3; j++){
C3[i][k]=(C3[i][j]+(A3[i][j]*B3[j][k])%5)%5;}}}
printf("vector Y[i+1]=F*Y[i]: \n");
for(i=0; i<S3; i++){
for(k=0; k<1; k++){
printf("%3d ",C3[i][k]);}
printf("\n");
printf("\n");}
break;
    
case 4:
int S4=4;
int A4[4][4];
int B4[4][1];
int C4[4][1];
for(j=0; j<S4; j++){
for(k=0; k<1; k++){
B4[j][k]=random(5);}}
 
for(i=0; i<S4; i++){
for(j=0; j<S4; j++){
A4[i][j]=random(5);}}
printf("matica F: \n");
for(i=0; i<S4; i++){
for(j=0; j<S4; j++){
printf("%3d ",A4[i][j]);}
   printf("\n");
   printf("\n");}
 
    printf("\n");
printf("vector Y[i]: \n");
for(j=0; j<S4; j++){
for(k=0; k<1; k++){
printf("%3d ",B4[j][k]);}
   printf("\n");
   printf("\n");}
 
   printf("\n");
 
for(i=0; i<S4; i++){
for(k=0; k<1; k++){
C4[i][k]=0;
for(j=0; j<S4; j++){
C4[i][k]=(C4[i][j]+(A4[i][j]*B4[j][k])%5)%5;}}}
printf("vector Y[i+1]=F*Y[i]: \n");
for(i=0; i<S4; i++){
for(k=0; k<1; k++){
printf("%3d ",C4[i][k]);}
printf("\n");
printf("\n");}
break;
    //return 0;
}   return 0;}

Программа достаточно кривая, но она у меня работает. Я не знаю как задать чтоб умножало дальше и выводило только требуемые результаты(задаются с клавиатуры по запросу программы). Ранг матрицы - 2х2, 3х3, 4х4(также задаётся с клавиатуры по запросу программы). Больше не нужно.
Буду рад любой помощи. Заранее спасибо.

@CheshireCat · 02.05.2009, 22:09

Хмм.... она у тебя компилируется? А какой компилятор?

@S.l.a.y.e.R · 03.05.2009, 14:15 **[ТС]**

компилятор С++

Добавлено через 1 минуту 30 секунд
компилирует без ошибок и предупреждений

@CheshireCat · 03.05.2009, 14:19

Вот поэтому я и спрашиваю, какой компилятор. Ибо GCC 4.4.0, например, наотрез отказывается это компилировать:
=== Build finished: 12 errors, 0 warnings ===

@S.l.a.y.e.R · 03.05.2009, 14:56 **[ТС]**

а где посмотреть какой именно у меня компилятор?

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .	Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .
Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .	SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .	SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip Сканируйте QR-код на мобильном и вы увидите, что появится джойстик для управления главным героем. . . .

@CheshireCat 2924 / 1274 / 114 Регистрация: 27.05.2008 Сообщений: 3,465
	02.05.2009, 22:09
	Хмм.... она у тебя компилируется? А какой компилятор? 0

@S.l.a.y.e.R 0 / 0 / 0 Регистрация: 02.05.2009 Сообщений: 10
	03.05.2009, 14:15 [ТС]
	компилятор С++ Добавлено через 1 минуту 30 секунд компилирует без ошибок и предупреждений 0

@CheshireCat 2924 / 1274 / 114 Регистрация: 27.05.2008 Сообщений: 3,465
	03.05.2009, 14:19
	Вот поэтому я и спрашиваю, какой компилятор. Ибо GCC 4.4.0, например, наотрез отказывается это компилировать: === Build finished: 12 errors, 0 warnings === 0

@S.l.a.y.e.R 0 / 0 / 0 Регистрация: 02.05.2009 Сообщений: 10
	03.05.2009, 14:56 [ТС]
	а где посмотреть какой именно у меня компилятор? 0