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

Длинная арифметика - C++

Восстановить пароль Регистрация
 
eocron
Кактус
 Аватар для eocron
66 / 66 / 6
Регистрация: 23.05.2012
Сообщений: 343
04.04.2013, 16:22     Длинная арифметика #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
          
class BigInt
{
        static const int max_size = 100000;//количество десятичных разрядов, которые должно вмещать
        static const int osn = 10000;//Основание - {0,1,....,9999}
        
        int len_osn;//длинная в чарах нашего основания
        int dig_size;//истинный размер в интах = max_size/len_osn
        
        int *digits;//массив разрядов по основанию osn
public:
          int amount;//размер нашего числа по osn
          ....
};
          friend BigInt operator - (const BigInt &a, const BigInt &b)
          {
                 BigInt res = a;
                 int r = 0;
                 for (int i = 0;i<res.amount;i++)
                 {
                        res.digits[i] -= b.digits[i] + r;
                        if (res.digits[i]<0)
                        {
                                res.digits[i]+=osn;
                                res.digits[i+1]--;
                        }   
                 }
                 int pos = a.amount;
                 while (pos>0 && !res.digits[pos])
                 {
                        pos--;
                 }
                 res.amount = pos+1;
 
                 return res;
          }
Ну, думаю, весь код здесь ни к чему. Проблема в том, что когда я провожу операцию a-b, массив a (именно массив digits в a) портится (появляется мусор в конце, причем совершенно рандомный при каждом запуске), все остальные характеристики остаются верными. Не пойму почему. Помогите.

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
          //Assignment    
          BigInt& operator = (const BigInt& other)
          {
                 if(this != &other)
                 {
                         memset(digits,0,sizeof(int)*dig_size);
                         amount = other.amount;
                         for(int i=0;i < other.amount;i++)
                         {
                                 digits[i] = other.digits[i];
                         }
                 }
                 return *this;
          }
Вот еще функция копирования
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
04.04.2013, 16:22     Длинная арифметика
Посмотрите здесь:

Длинная арифметика C++
C++ Длинная арифметика
C++ Длинная арифметика
C++ Длинная арифметика
C++ Длинная арифметика
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
skiff26
2 / 2 / 0
Регистрация: 15.11.2011
Сообщений: 47
04.04.2013, 16:32     Длинная арифметика #2
Приведи тестовый пример на котором от рабатываешь класс, ну и если не сложно дай архив класса для ознакомления, т.к. ты изобретаешь колесо.
iifat
2179 / 1332 / 96
Регистрация: 05.06.2011
Сообщений: 3,690
04.04.2013, 16:33     Длинная арифметика #3
Интересно, конструктор копирования есть? Если нет, вычитание портит первый аргумент. Также интересен бы деструктор: если он освобождает память под digits, то он ещё и освободит память вычитаемого и туда и правда может любой мусор залететь...
eocron
Кактус
 Аватар для eocron
66 / 66 / 6
Регистрация: 23.05.2012
Сообщений: 343
04.04.2013, 16:36  [ТС]     Длинная арифметика #4
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
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
#include <iostream>
#include <string>
#include <math.h>
 
using namespace std;
 
class BigInt
{
        static const int max_size = 100000;// òðåáóåìûé Г°Г*çìåð Гў 10-îé ÑÑ
        static const int osn = 10000;//îñГ*îâГ*Г*ГЁГҐ ÑÑ
        
        int len_osn;//äëèГ*Г*Г* îñГ*îâГ*Г*ГЁГї ÑÑ
        int dig_size;//ГЁГ±ГІГЁГ*Г*ûé Г°Г*çìåð Г¬Г*Г±Г±ГЁГўГ*
        
        int *digits;
public:
          int amount;//êîëè÷åñòâî Г§Г*Г*ÿòûõ Г°Г*çðÿäîâ Г°Г*çìåðГ* len_osn
          //int to BigInt
          BigInt(int a):amount(0)
          {
                   start_initialize();
                   
                   int* curr = digits;
                   do 
                   {
                        *curr++ = a % osn;
                        a /= osn;
                        amount++;
                   }while(a != 0);
          }
          //Simple construct
          BigInt():amount(0)
          {
                   start_initialize();
          }
          ~BigInt()
          {
                   free(digits);
          }
          //Assignment    
          BigInt& operator = (const BigInt& other)
          {
                 if(this != &other)
                 {
                         memset(digits,0,sizeof(int)*dig_size);
                         amount = other.amount;
                         for(int i=0;i < other.amount;i++)
                         {
                                 digits[i] = other.digits[i];
                         }
                 }
                 return *this;
          }
//###===###===###===###===###===###===###===###===###===###===###===###===###
//###                             >> or <<                                ###     
//###===###===###===###===###===###===###===###===###===###===###===###===###
          friend ostream &operator<<(ostream& output, BigInt& a)
          {
                  a.output(output);
                  return output;
          }
          
          friend istream &operator>>(istream& input, BigInt& a)
          {
                  a.input(input);
                  return input;
          }
//###===###===###===###===###===###===###===###===###===###===###===###===###
//###                             + or -                                  ###     
//###===###===###===###===###===###===###===###===###===###===###===###===###
          friend BigInt operator + (const BigInt &a, const BigInt &b)
          {
                  BigInt res;
                  res.amount = max(a.amount,b.amount);
                  int r = 0;
                  for(int i=0;i<res.amount | r;i++)
                  {
                        res.digits[i] = a.digits[i] + b.digits[i] + r;
                        if (res.digits[i]>=osn)
                        {
                              res.digits[i]-=osn;
                              r = 1;
                        }
                        else
                        {
                              r = 0;
                        }
                  }
                  if(res.digits[res.amount])
                  {
                     res.amount++; 
                  }
                  return res;
          }
          
          friend BigInt operator - (const BigInt &a, const BigInt &b)
          {
                 BigInt res = a;
                 int r = 0;
                 for (int i = 0;i<res.amount;i++)
                 {
                        res.digits[i] -= b.digits[i] + r;
                        if (res.digits[i]<0)
                        {
                                res.digits[i]+=osn;
                                res.digits[i+1]--;
                        }   
                 }
                 int pos = a.amount;
                 while (pos>0 && !res.digits[pos])
                 {
                        pos--;
                 }
                 res.amount = pos+1;
                 cout<<"++++++++++++++++++++++++++"<<endl;
                 cout<<res<<endl;
                 return res;
          }
//###===###===###===###===###===###===###===###===###===###===###===###===###
//###                             * or /                                  ###     
//###===###===###===###===###===###===###===###===###===###===###===###===###
          friend BigInt operator * (const BigInt &a, const BigInt &b)
          {
                  BigInt res;
                  for(int i=0;i<a.amount;i++)
                  {
                        int r = 0;
                        for (int j=0;j<b.amount | r;j++)
                        {
                              res.digits[i+j] += a.digits[i] * b.digits[j] + r;
                              //cout<<a.digits[i]<<'*'<<b.digits[j]<<'-'<<r<<'='<<res.digits[i+j]<<endl;
                              r = res.digits[i+j] / osn;
                              res.digits[i+j] -= r*osn;
                        }
                  }
                  int pos = a.amount + b.amount;
                  while (pos>0 && !res.digits[pos])
                  {
                        pos--;
                  }
                  res.amount = pos + 1;
                  cout<<res.amount<<endl;
                  return res;
          }
 
          friend BigInt operator / (const BigInt &a, const BigInt &b)
          {
          }
            
//###===###===###===###===###===###===###===###===###===###===###===###===###
//###                             BOOLSHIT                                ###     
//###===###===###===###===###===###===###===###===###===###===###===###===###            
            friend bool operator == (const BigInt &a, const BigInt &b)
            {
                  if (a.amount!=b.amount)
                  {
                        return false;
                  }
                  for (int i=0;i<a.amount;i++)
                  {
                        if (a.digits[i]!=b.digits[i])
                        {
                               return false;
                        }
                  }
                  return true;
            }
            
            friend bool operator > (const BigInt &a, const BigInt &b)
            {
                  if (a.amount!=b.amount)
                  {
                        return a.amount>b.amount;
                  }
                  for (int i=a.amount-1;i>=0;i--)
                  {
                        if (a.digits[i]!=b.digits[i])
                        {
                              return a.digits[i]>b.digits[i];
                        }
                  }
                  return false;
            }
            
            friend bool operator < (const BigInt &a, const BigInt &b)
            {
                  if (a.amount!=b.amount)
                  {
                       return a.amount<b.amount;
                  }
                  for (int i=a.amount-1;i>=0;i--)
                  {
                        if (a.digits[i]!=b.digits[i])
                        {
                              return a.digits[i]<b.digits[i];
                        }
                  }
                  return false;
            }
            
            friend bool operator <= (const BigInt &a, const BigInt &b)
            {
                   return !(a>b);
            }
            
            friend bool operator >= (const BigInt &a, const BigInt &b)
            {
                   return !(a<b);
            }
private: 
          //Some Junk for >>
          void input(istream& input)
          {
                  memset(digits,0,sizeof(int)*dig_size);
                  amount = 0;
                  string str;
                  int pos = 0;
                  input>>str;
                  
                  for (int i=str.size()-1;i>=0;i-=len_osn)
                  {
                      int start = i - len_osn + 1;
                      if(start<0)
                      {
                             start = 0;
                      }
                      string dig = str.substr(start,i-start+1);
                      digits[pos++] = atoi(dig.c_str());
                  }
                  //Cut off zero's
                  amount = str.size();
                  for(int i=str.size()-1 ; i>0 ; i--)
                  {
                          if(digits[i]!=0)
                          {
                                break;
                          }
                          else
                          {
                                amount--;
                          }
                  }
          }
          //Some Junk for <<
          void output(ostream& output)
          {
                      if(amount==0){output << 0;}
                      cout<<digits[amount - 1];
                      for(int i=amount - 2;i >= 0;i--)
                      {
                              printf("%.*d",len_osn,digits[i]);
                      }
          }
          //Some Junk for /
          void LevelUp()
          {
                 digits[0]=0;
                 for(int i = amount;i>=1;i--)
                 {
                        digits[i] = digits[i-1];
                 }
                 if(digits[amount])
                 {
                        amount++;
                 }
          }
          void Error()
          {
                 cout<<"Error\n";
          }
          void start_initialize()
          {
                   len_osn = (int)log10((float)osn);
                   dig_size = max_size/len_osn;
                   
                   digits=(int *)malloc( sizeof(int)*dig_size );
                   memset(digits,0, sizeof(int)*dig_size );
                   //cout<<"len_osn:"<<len_osn<<' '<<"dig_size:"<<dig_size<<' '<<"osn:"<<osn<<' '<<"max_size:"<<max_size<<endl;
          }
};
int main()
{
    BigInt a,c,b;
    //int d;
    cin>>a;
    cin>>b;
/*    cout<<'\n'<<a<<' '<<b<<endl<<a.amount<<' '<<b.amount<<endl;
    c=a+b;
    cout<<"a+b:"<<c<<endl;
    cout<<'\n'<<a<<' '<<b<<endl<<a.amount<<' '<<b.amount<<endl;*/
    c=a-b;
    cout<<"a-b:"<<c<<endl;
    cout<<'\n'<<a<<' '<<b<<endl<<a.amount<<' '<<b.amount<<endl;
 /*    c=a*b;
    cout<<"a*b:"<<c<<endl;
    cout<<'\n'<<a<<' '<<b<<endl<<a.amount<<' '<<b.amount<<endl;
    c=a/b;
    cout<<"a/b:"<<c<<endl;*/
 
    system("pause");
}
Окай, ловите.

Тест такой:

1736816287368162387162836186381 1000000000000000
skiff26
2 / 2 / 0
Регистрация: 15.11.2011
Сообщений: 47
04.04.2013, 16:41     Длинная арифметика #5
смотрю у тебя есть не явный конструктор, воспользуйся для начала им: а(234), б(456).
eocron
Кактус
 Аватар для eocron
66 / 66 / 6
Регистрация: 23.05.2012
Сообщений: 343
04.04.2013, 16:49  [ТС]     Длинная арифметика #6
Все, понял. Спасибо.
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
05.04.2013, 10:00     Длинная арифметика
Еще ссылки по теме:

Длинная арифметика C++
C++ Длинная арифметика

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

Или воспользуйтесь поиском по форуму:
skiff26
2 / 2 / 0
Регистрация: 15.11.2011
Сообщений: 47
05.04.2013, 10:00     Длинная арифметика #7
При сложении чисел для res вызывается декструктор. Возможно некорректно перегружены операторы или нехватает перегруженного оператора, уточнить надо по хелпу или по какой-либо книжке по С++(T operoator..., T& operator... - разные вещи).
И может не стоит сразу гнаться за производительность используя С-функции работы с памятью;
PS. В С++ использую new и delete.
Yandex
Объявления
05.04.2013, 10:00     Длинная арифметика
Ответ Создать тему
Опции темы

Текущее время: 00:38. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru