Форум программистов, компьютерный форум CyberForum.ru

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
 
Рейтинг: Рейтинг темы: голосов - 10, средняя оценка - 4.90
Larri
0 / 0 / 0
Регистрация: 29.09.2011
Сообщений: 16
#1

При наличие корней в системе уравнений выдает ошибку - C++

29.10.2011, 20:07. Просмотров 1212. Ответов 22
Метки нет (Все метки)

После выполнения программы (она решает систему линейных уравнений) для случая наличия корней появляется ошибка.
Вот код:
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
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
#include <stdio.h>
#include<stdlib.h>
#include<time.h>
#define N 4
typedef struct
{int*a;
bool l;//нулевой ли вектор
}lama;
int* WC(int *w1,int f,int n)//alfa na vector
{
    int i=0;
    int *w;
    w=(int*)malloc(n*sizeof(int));
    for(i=0;i<n;i++)
        w[i]=0;
    for(i=0;i<n;i++)
        w[i]=(w[i]+w1[i]*f)%2;
    return w;
}
int WW(int *w1,int *w2,int n)//skal proizv
{
    int i=0;
    int w=0;
    for(i=0;i<n;i++)
        w=(w+w1[i]*w2[i])%2;
    return w;
}
int* MW(int n, int *a, int *w)
{
    int i=0,j=0;
    int *w1;
    w1=(int*)malloc(n*sizeof(int));
    for (j=0;j<n;j++)
        w1[j]=0;
    for (j=0;j<n;j++)
        for(i=0;i<n;i++)
            w1[j]=(w1[j]+a[j*n+i]*w[i])%2;
    return w1;
}
lama form(int n,int *a,int *w1, int *w0)
{
    int i=0;
    int *w2;
    lama w;
    w.l=false;  
    int q,q1,e=0;
    w.a=(int*)malloc(n*sizeof(int));
    w2=(int*)malloc(n*sizeof(int));
    w2=MW(n,a,w1);
    q=(WW(MW(n,a,w2),w1,n))*(WW(w2,w1,n));  
    if(w0==NULL)
        for(i=0;i<n;i++)
        {
            w.a[i]=(w2[i]+WC(w1,q,n)[i])%2; //в двоичном сложении 1+1=1-1=0, если не учитывать перенос разряда
            e+=w.a[i];
        }
    else
    {
        q1=(WW(MW(n,a,w2),w0,n))*(WW(MW(n,a,w0),w0,n));
        for(i=0;i<n;i++)
        {
            w.a[i]=(w2[i]+WC(w1,q,n)[i]+WC(w0,q1,n)[i])%2;
            e+=w.a[i];          
        }
    }   
    if(e==0) w.l=true;
    free(w2);
    return w;
}
int*XX(int*a,int*w,int k,int n)//nepr
{
    int j,i,q;
    int*x;
    x=(int*)malloc(n*sizeof(int));  
    for(i=0;i<n;i++)
        x[i]=0;
    for(i=0;i<=k;i++)
    {       
        q=WW((w+i*n),w,n)*WW(MW(n,a,(w+i*n)),(w+i*n),n);
        for(j=0;j<n;j++)
            x[j]=(x[j]+q*w[i*n+j])%2;
    }   
    return x;
}
bool prov(int*a,int*x,int*b,int n)
{   
    int*b1,i,j=0;
    b1=(int*)malloc(n*sizeof(int));
    b1=MW(n,a,x);
    for(i=0;i<n;i++)
        if(b1[i]!=b[i]) j++;
    free(b1);
    if(j==0) return true;
    else return false;
}
int main()
{
    int *a,*b,*w,*x,i,j,k=0;
    lama p; 
    a=(int*)malloc(N*N*sizeof(int));
    b=(int*)malloc(N*sizeof(int));
    w=(int*)malloc(N*N*sizeof(int));
    /*srand(time(0));
    for(i=0;i<N;i++)
    {
        for(j=i;j<N;j++)
        {
            a[i*N+j]=rand()%2;
            if(i!=j) a[j*N+i]=a[i*N+j];
        }
        b[i]=rand()%2;
    }
    for (i=0; i<N; i++)
    {
        for (j=0; j<N; j++)
            printf("\t %i", a[i*N+j]);
        printf ("\n");
        }   */
    for(i=0;i<N;i++)//-пользовательское задание матрицы А
        for(j=0;j<N;j++)
        {
            printf("Enter a[%i,%i]\n",i,j);
            scanf("%i",&a[i*N+j]);
        }
    for(i=0;i<N;i++)
    {
        printf("Enter b[%i]\n",i);
        scanf("%i",&b[i]);
    }
    for (i=0; i<N; i++)
    {
        for (j=0; j<N; j++)
            printf("\t %i", a[i*N+j]);
        printf ("\n");
    }
    for(i=0;i<N;i++)    
        printf("b[%i]=%i\n",i,b[i]);
    for(j=0;j<N;j++)
        w[j]=b[j];
    i=1;
    do
    {
        if(i==1) p=form(N,a,b,NULL);
        else p=form(N,a,(w+(i-1)*N),(w+(i-2)*N));
        for(j=0;j<N;j++)//было j<=N
        {
            w[i*N+j]=p.a[j];//запись вектора W[i] в матрицу омег            
        }
        k++;
        i++;
    }
    while(!p.l);
    if(p.l)
    {
        for (i=0; i<=k; i++)
        {
            for (j=0; j<N; j++)
                printf("\t %i", w[i*N+j]);
            printf ("\n");
        }
        x=XX(a,w,k,N);      
        if(!prov(a,x,b,N)) printf("answer does not exist\n");
        else 
        {
            printf("answer is true\n");
            printf("answer: ");
            for(i=0;i<N;i++)
                printf("%i ",x[i]);
            printf("\n");
        }
    }
    else printf("algorithm does not work\n");
    free(a);
    free(b);
    free(w);
    return 0;
}
Лучшие ответы (1)
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
29.10.2011, 20:07     При наличие корней в системе уравнений выдает ошибку
Посмотрите здесь:

При b=1 выдает ошибку - C++
# include &lt;iostream&gt; # include &lt;math.h&gt; using namespace std; double Func( double s, double b); void main () { double c=0;...

Выдает ошибку при построении - C++
List.h #pragma once #include &lt;iostream&gt; using namespace std; struct A { int key; }; struct listStruct { A a;

Почему выдает ошибку при вводе? - C++
Вот моя программа. Она перестала работать после того как я добавил функцию (invers). Эта функция должна инвертировать второй список. ...

программа выдает ошибку при компиляции - C++
подскажите, пожалуста, в чем ошибка..При компиляции выскакивает ошибка undefined reference to 'Obrabotchik::SetListContinents()' ...

Выдает ошибку при создании модуля - C++
у меня задача создать программу и 3 модуля с расширение .h я создал все три, первые 2 отлично работают, но когда я создал и подключил...

При заполнении двумерного массива выдает ошибку - C++
Недавно начал изучать С++ и застрял на двумерном массиве. решил сделать двумерный массив в виде лабиринта, но компилятор выдает ошибку на...

Почему выдает такую ошибку? c++ при компиляции - C++
Вот код программы, нахождение наибольшего числа с 3 заданных #include &lt;iostream&gt; uisng namespace std; int max(int a,int b) ...

После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
alkagolik
Заблокирован
29.10.2011, 21:25     При наличие корней в системе уравнений выдает ошибку #16
Сообщение было отмечено автором темы, экспертом или модератором как ответ
-=ЮрА=-, просто руками писать надо, а не чем-либо еще
тест твоей программы
deadman@kubenta:~/documents/projects/system/bin/Debug$ valgrind ./system
==2486== Memcheck, a memory error detector
==2486== Copyright (C) 2002-2010, and GNU GPL'd, by Julian Seward et al.
==2486== Using Valgrind-3.6.1-Debian and LibVEX; rerun with -h for copyright info
==2486== Command: ./system
==2486==
f1 :
0 1 2 3 4 5 6 7 8 9
f2 :
==2486== Invalid read of size 4
==2486== at 0x8048595: main (main.c:28)
==2486== Address 0x41c5084 is 4 bytes inside a block of size 5 alloc'd
==2486== at 0x402896C: realloc (vg_replace_malloc.c:525)
==2486== by 0x80484FA: f2 (main.c:14)
==2486== by 0x804856F: main (main.c:25)
==2486==
0 1 0 0 0 ==2486==
==2486== HEAP SUMMARY:
==2486== in use at exit: 0 bytes in 0 blocks
==2486== total heap usage: 2 allocs, 2 frees, 45 bytes allocated
==2486==
==2486== All heap blocks were freed -- no leaks are possible
==2486==
==2486== For counts of detected and suppressed errors, rerun with: -v
==2486== ERROR SUMMARY: 4 errors from 1 contexts (suppressed: 11 from 6)

нормальный код
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
#include <stdio.h>
#include <stdlib.h>
 
int * f(int n)
{
    int * x = (int *)malloc(n*sizeof(int));
        for(int i = 0; i < n; i++)
                x[i] = i;
    return x;
}
 
int * f2(int * x, int m)
{
    int *p = x;
    p = realloc(p, m * sizeof(int));
    return p;
}
 
 
int main ()
{
    int i, * a = f(10);
        printf("f1 : \n");
    for(i = 0; i < 10; i++)
                printf("%d ",a[i]);
        a = f2(a, 5);
        printf("\nf2 : \n");
        for(i = 0; i < 5; i++)
                printf("%d ",a[i]);
    free(a);
    return 0;
}
тест моей программы

deadman@kubenta:~/documents/projects/system/bin/Debug$ valgrind ./system
==2623== Memcheck, a memory error detector
==2623== Copyright (C) 2002-2010, and GNU GPL'd, by Julian Seward et al.
==2623== Using Valgrind-3.6.1-Debian and LibVEX; rerun with -h for copyright info
==2623== Command: ./system
==2623==
f1 :
0 1 2 3 4 5 6 7 8 9
f2 :
0 1 2 3 4 ==2623==
==2623== HEAP SUMMARY:
==2623== in use at exit: 0 bytes in 0 blocks
==2623== total heap usage: 2 allocs, 2 frees, 60 bytes allocated
==2623==
==2623== All heap blocks were freed -- no leaks are possible
==2623==
==2623== For counts of detected and suppressed errors, rerun with: -v
==2623== ERROR SUMMARY: 0 errors from 0 contexts (suppressed: 11 from 6)
Larri
0 / 0 / 0
Регистрация: 29.09.2011
Сообщений: 16
29.10.2011, 21:36  [ТС]     При наличие корней в системе уравнений выдает ошибку #17
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение

Не по теме:

Larri, а у вас нет ничего из теории???Это ускорило бы реализацию алгоритма
PS:Мне нужно матоописание не потому что не хочу править Ваш код, а потому что в силу как я понял сложности метода возможно по самому алгоритму что-то упустили, вот Я не хочу обидеть(мол не так написали) просто хочу 100% как логическую так и ошибку реализации исключить

Не по теме:

Вот ссылка на сам алгоритм, но у меня модификация (по ней я вообще ничего не нашла):http://lib.mexmat.ru/books/4633
Вероятность логической ошибки мала, т к 1. программа вычисляет решение правильно, 2. выдаваемая ошибка идет с сообщением Heap Corruption Detected.
А ошибка реализации только если где-то при работе с памятью...

-=ЮрА=-
Заблокирован
Автор FAQ
31.10.2011, 11:55     При наличие корней в системе уравнений выдает ошибку #18
Larri, для начала решил оттестировать ваш алгоритм на счёт простой системы, как видите на скрине в алгоритме скорее всего логическая ошибка, т.к. всегда попадается b1[i] != b[i]
Цитата Сообщение от Larri Посмотреть сообщение
if(b1[i]!=b[i]) j++;
Ваш код с небоьшими изменениями, в плане его сокращения откомпилировал несколько раз в Release, даже на простой системе Ваш алгоритм не работает. Вообщем у меня многое под сомнением, но вот здесь особенно
Цитата Сообщение от Larri Посмотреть сообщение
w1[j]=(w1[j]+a[j*n+i]*w[i])%2;
a[j*n+i] - допустим j = n - 1 i = n тогда вы пытаетесь обратиться к несуществующему элементу
a[(n - 1)*n + (n - 1)] как я понял по программе в а всего n элементов.
Цитата Сообщение от Larri Посмотреть сообщение
http://lib.mexmat.ru/books/4633
Там книжку читать надо, а быстрый поиск по Ланцош ничего не дал, если можно сделайте принт скрин нужных страниц и сюда тогда смогу проверить что к чему

Цитата Сообщение от Larri Посмотреть сообщение
w=(w+w1[i]*w2[i])%2;
- ну и это мне не понятно, вы двоичное сложение хотите сделать или что этим хотели реализовать???Это не двоичное сложение - а остаток от целочисленного деления на два, поясню
2 + 3 = 5 (10 + 11 = 101)
5 % 2 = 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
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#define N 4
 
typedef struct
{
    int*a;
    bool l;//нулевой ли вектор
}
lama;
 
int* WC(int *w1,int f,int n)//alfa na vector
{
    int i=0;
    int *w=(int*)malloc(n*sizeof(int));
    for(i=0;i<n;i++)
        w[i]=0;
    for(i=0;i<n;i++)
        w[i]=(w[i]+w1[i]*f)%2;
    return w;
}
 
int WW(int *w1,int *w2,int n)//skal proizv
{
    int i=0;
    int w=0;
    for(i=0;i<n;i++)
        w=(w+w1[i]*w2[i])%2;
    return w;
}
 
int* MW(int n, int *a, int *w)
{
    int i,j;
    int *w1=(int*)malloc(n*sizeof(int));
    for (j=0;j<n;j++)
        w1[j]=0;
    for (j=0;j<n;j++)
    for (i=0;i<n;i++)
        w1[j]=(w1[j]+a[j*n+i]*w[i])%2;
    return w1;
}
 
lama form(int n,int *a,int *w1, int *w0)
{
        int i=0;
        int *w2;
        lama w;
        w.l=false;      
        int q,q1,e=0;
        w.a=(int*)malloc(n*sizeof(int));
        w2=(int*)malloc(n*sizeof(int));
        w2=MW(n,a,w1);
        q=(WW(MW(n,a,w2),w1,n))*(WW(w2,w1,n));  
        if(w0==NULL)
        {
                for(i=0;i<n;i++)
                {
                        w.a[i]=(w2[i]+WC(w1,q,n)[i])%2; //в двоичном сложении 1+1=1-1=0, если не учитывать перенос разряда
                        e+=w.a[i];
                }
        }
        else
        {
                q1=(WW(MW(n,a,w2),w0,n))*(WW(MW(n,a,w0),w0,n));
                for(i=0;i<n;i++)
                {
                        w.a[i]=(w2[i]+WC(w1,q,n)[i]+WC(w0,q1,n)[i])%2;
                        e+=w.a[i];                      
                }
        }       
        if(e==0) 
            w.l=true;
        free(w2);
        return w;
}
int*XX(int*a,int*w,int k,int n)//nepr
{
        int j,i,q;
        int*x =(int*)malloc(n*sizeof(int));  
        for(i=0;i<n;i++)
            x[i]=0;
        for(i=0;i<=k;i++)
        {               
                q=WW((w+i*n),w,n)*WW(MW(n,a,(w+i*n)),(w+i*n),n);
                for(j=0;j<n;j++)
                    x[j]=(x[j]+q*w[i*n+j])%2;
        }       
        return x;
}
bool prov(int*a,int*x,int*b,int n)
{       
    int i,j = 0;
    int * b1=MW(n,a,x);
    for(i=0;i<n;i++)
    if(b1[i]!=b[i])
        j++;
    return (j == 0 ? true : false);
}
int main()
{
        int *x,i,j,k=0;
        lama p; 
        int * a=(int*)malloc(N*N*sizeof(int));
        int * b=(int*)malloc(N*sizeof(int));
        int * w=(int*)malloc(N*N*sizeof(int));
        for(i=0;i<N;i++)//-пользовательское задание матрицы А
                for(j=0;j<N;j++)
                {
                        printf("Enter a[%i,%i] : ",i,j);
                        scanf("%i",&a[i*N+j]);
                }
        for(i=0;i<N;i++)
        {
                printf("Enter b[%i] : ",i);
                scanf("%i",&b[i]);
        }
        for (i=0; i<N; i++)
        {
                for (j=0; j<N; j++)
                        printf("\t %i", a[i*N+j]);
                printf ("\n");
        }
        for(i=0;i<N;i++)        
                printf("b[%i]=%i\n",i,b[i]);
        for(j=0;j<N;j++)
                w[j]=b[j];
        i=1;
        do
        {
                if(i==1) 
                    p=form(N,a,b,NULL);
                else 
                    p=form(N,a,(w+(i-1)*N),(w+(i-2)*N));
                for(j=0;j<N;j++)//было j<=N
                {
                        w[i*N+j]=p.a[j];//запись вектора W[i] в матрицу омег                    
                }
                k++;
                i++;
        }
        while(!p.l);
        if(p.l)
        {
                for (i=0; i<=k; i++)
                {
                        for (j=0; j<N; j++)
                                printf("\t %i", w[i*N+j]);
                        printf ("\n");
                }
                x=XX(a,w,k,N); 
                prov(a,x,b,N);
               if(!prov(a,x,b,N)) 
                    printf("answer does not exist\n");
                else 
                {
                        printf("answer is true\n");
                        printf("answer: ");
                        for(i=0;i<N;i++)
                                printf("%i ",x[i]);
                        printf("\n");
                }
        }
        else 
            printf("algorithm does not work\n");
 /*       free(a);
        free(b);
        free(w);*/
        return 0;
}
Миниатюры
При наличие корней в системе уравнений выдает ошибку  
Larri
0 / 0 / 0
Регистрация: 29.09.2011
Сообщений: 16
31.10.2011, 20:10  [ТС]     При наличие корней в системе уравнений выдает ошибку #19
Дело в том, что система на писана для системы над кольцом вычетов по модулю 2. Соответственно, там есть только числа 0 и 1, поэтому приведенный тест, естественно, приведет к неправильному результату. Я проверяла для систем над нужным концом, результат либо верный, либо программа падает. В а n*n элементов, т.к. это квадратная матрица
-=ЮрА=-
Заблокирован
Автор FAQ
01.11.2011, 10:22     При наличие корней в системе уравнений выдает ошибку #20
Цитата Сообщение от Larri Посмотреть сообщение
Дело в том, что система на писана для системы над кольцом вычетов по модулю 2. Соответственно, там есть только числа 0 и 1, поэтому приведенный тест, естественно, приведет к неправильному результату.
- т.е работаете с Системой н-го порядка с 1-цами и нулями, а как же тогда вектор свободных членов, что тоже только 1 и 0 в нём???


Цитата Сообщение от Larri Посмотреть сообщение
Я проверяла для систем над нужным концом, результат либо верный, либо программа падает.
- предлагаю заменить все массивы статическими, вы же для системы 4-го порядка алгоритм затачиваете
Larri
0 / 0 / 0
Регистрация: 29.09.2011
Сообщений: 16
01.11.2011, 15:31  [ТС]     При наличие корней в системе уравнений выдает ошибку #21
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение
- т.е работаете с Системой н-го порядка с 1-цами и нулями, а как же тогда вектор свободных членов, что тоже только 1 и 0 в нём???


- предлагаю заменить все массивы статическими, вы же для системы 4-го порядка алгоритм затачиваете
Да, система н-го поряка с 1-ми и 0-ми, в свободном векторе тоже 0-ли и 1-цы. Заменить массивы статическими не получится, т. к. массивы потом будут разных размерностей.
Программа компелируется, но завершается с ошибкой: CRT detected that the application wrote to memory after end of heap buffer, и как это исправить я не знаю=(
Deviaphan
Делаю внезапно и красиво
Эксперт C++
1286 / 1220 / 50
Регистрация: 22.03.2011
Сообщений: 3,744
01.11.2011, 15:43     При наличие корней в системе уравнений выдает ошибку #22
Т.е. все эти массивы содержат только нули и единицы? Тогда не нужно всех этих ужасов с (+) и (%), просто ксорь (^).
И ищи, где портишь кучу. Только сперва избавься от утечек памяти, я уже показывал где именно.
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
07.11.2011, 11:04     При наличие корней в системе уравнений выдает ошибку
Еще ссылки по теме:

Компилятор выдает ошибку при создании списка - C++
Доброго времени суток всем, не могу понять где ошибка создаю список list самым простым способом #include &lt;list&gt; #include...

Выдает ошибку при передаче массива в функцию - C++
Ругается на 59 строку: 2.cpp(59): error C2664: otn: невозможно преобразовать параметр 4 из &quot;double **&quot; в &quot;double&quot; Не существует...

Компилятор выдает ошибку при использовании фабричного метода - C++
Всем привет!!! новая ошибка... смотрите, теперь у меня есть класс template &lt;typename T&gt; struct Node{ T key; Node *...

Выдает ошибку при компиляции. Не могу понять почему - C++
Помогите с кодом. Я вижу что ошибка уровня самых низов но не могу понять все равно что тут не так. Вот эта ошибка: // error...

Почему выдает ошибку при передаче по ссылке параметров? - C++
#include &lt;iostream&gt; /* run this program using the console pauser or add your own getch, system(&quot;pause&quot;) or input loop */ void...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
Larri
0 / 0 / 0
Регистрация: 29.09.2011
Сообщений: 16
07.11.2011, 11:04  [ТС]     При наличие корней в системе уравнений выдает ошибку #23
Цитата Сообщение от Deviaphan Посмотреть сообщение
Т.е. все эти массивы содержат только нули и единицы? Тогда не нужно всех этих ужасов с (+) и (%), просто ксорь (^).
И ищи, где портишь кучу. Только сперва избавься от утечек памяти, я уже показывал где именно.
Да, массивы содержат только нули и единицы.
Насчет кучи постараюсь.
Yandex
Объявления
07.11.2011, 11:04     При наличие корней в системе уравнений выдает ошибку
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru