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

Перевести код из С++ в Си - C++

Войти
Регистрация
Восстановить пароль
Другие темы раздела
C++ Как читать из файла http://www.cyberforum.ru/cpp-beginners/thread10925.html
Pomogite pojalujsta! Ia ne znaiu kak chitat iz faila. Naprimer kak chitat chto nmapisanno v word faile i uvidet etot rezultat kak parametr.
C++ Двоичный код символа мне нужно из символа, после считывания из файла получить его двоичный код!!! подскажите плиз как это сделать!!! http://www.cyberforum.ru/cpp-beginners/thread10898.html
Помогите с функцией rand C++
Привет всем помогите решить простую задачу Для каждого из следующих наборов целых чисел напишите единственный оператор, который будет печатать случайно выбранное число из набора: a) 2, 4, 6, 8, 10 b) 3, 5, 7, 9, 11 С) 6, 10, 14, 18, 22 т.е. как задавать шаг выбора случайного числа
Объясните, в чем суть директив #ifndef/#define/#endif C++
Не понял в чем суть директив: #ifndef .. #define.. #endifОбъясните пожалуйста..это как - то связано с хедерами, но вот как ??
C++ Секреты инкремента http://www.cyberforum.ru/cpp-beginners/thread10705.html
Пожалуйста!!! помогите! Умоляю!:help: i=5; j=(++i) + (++i): Чему равно j?; пожалуйста обьясните как это решается, ну или просто ответ скаите???:help::huh:
C++ Реализация программы с помощью классов Здравствуйте, коллеги-программеры! Обращаюсь к Вам с таким вопросом. У меня есть написанная мной неплохая демка, с использованием OpenGL, которая случайным образом генерирует звёзды (точки в 3D-пространстве) и движет их на зрителя. Демка на подобие того скринсейвера, который в Виндовс называется «сквозь вселенную». Программа работает, всё нормально. Но я хотел бы реализовать эту прогу с... подробнее

Показать сообщение отдельно
Splendid
0 / 0 / 1
Регистрация: 27.03.2008
Сообщений: 75

Перевести код из С++ в Си - C++

06.08.2008, 14:50. Просмотров 2017. Ответов 11
Метки (Все метки)

Есть код на С++, надо чистый С. Пыталась исправлять ошибки, но чем больше я это делаю, тем больше ошибок он выдает...поэтому привожу изначальный не исправленный на С++. Пожалуйста, помогите разобраться!
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
#include <algorithm>
#include <cmath>
#include <iostream>
#include <map>
#include <vector>
 
using namespace std;
 
 
//! Возвращает true, если n - точный квадрат простого числа
int perfect_square (const n)
{
int sq = ceil (sqrt ((double)n));
    return sq*sq == n;
}
 
 
//! Возвращает количество бит в числе (т.е. минимальное количество бит, которыми можно представить данное число)
 
unsigned bits_in_number (int n)
{
    if (n == 0)
  return 1;
    unsigned result = 0;
    while (n)
    {
  n /= 2;
  ++result;
    }
    return result;
}
 
//! Возвращает значение k-го бита числа (биты нумеруются с нуля)
 
int test_bit (const n, unsigned k)
{
    return (n & ((1) << k)) != 0;
}
 
//! Умножает a *= b (mod n)
 
void mulmod (int & a, int b, const & n)
{
    // наивная версия, годится только для длинной арифметики
    a *= b;
    a %= n;
}
 
 
//! Переводит число n в форму q*2^p
void transform_num (int n, int & p, int & q)
{
    int p_res = 0;
    while (n % 2 == 0)
    {
  ++p_res;
  n /= 2;
    }
    p = p_res;
    q = n;
}
 
 
int gcd(int a, int b)
{
    int c;
    while (b)
    {
        c = a % b;
        a = b;
        b = c;
    }
    return a;
}
 
int jacobi(long d, long n)
 
{
  int s;
  long d1, b = d, e = 0, m, n1;
 
 
  if (d == 0) return 0;
  if (d == 1) return 1;
 
  while ((b & 1) == 0)
    b >>= 1, e++;
  d1 = b;
  m = n % 8;
 
  if (!(e & 1)) s = 1;
 
  else if (m == 1 || m == 7) s = + 1;
 
  else if (m == 3 || m == 5) s = - 1;
 
  if (n % 4 == 3 && d1 % 4 == 3) s = - s;
 
  if (d1 != 1) n1 = n % d1; else n1 = 1;
 
  return s * jacobi(n1, d1);
 
}
 
//! Усиленный алгоритм Лукаса-Селфриджа проверки n на простоту. Используется усиленный алгоритм Лукаса с параметрами Селфриджа. Работает только с знаковыми типами!!! Второй параметр unused не используется, он только дает тип
 
int lucas_selfridge (const & n)
{
 
    // сначала проверяем тривиальные случаи
    if (n == 2)
  return true;
    if (n < 2 || (n % 2 == 0))
  return false;
 
    // проверяем, что n не является точным квадратом, иначе алгоритм даст ошибку
    if (perfect_square (n))
  return false;
 
    // алгоритм Селфриджа: находим первое число d такое, что:
    // jacobi(d,n)=-1 и оно принадлежит ряду { 5,-7,9,-11,13,... }
    int dd;
    for (int d_abs = 5, d_sign = 1;; d_sign = -d_sign, ++++d_abs)
    {
  dd = d_abs * d_sign;
  int g = gcd (n, d_abs);
  if (1 < g && g < n)
    // нашли делитель - d_abs
    return false;
  if (jacobi (int(dd), n) == -1)
    break;
    }
 
    // параметры Селфриджа
    int
  p = 1,
  q = (p*p - dd) / 4;
    
    // разлагаем n+1 = d*2^s
    int n_1 = n;
    ++n_1;
    int s, d;
    transform_num (n_1, s, d);
 
    // алгоритм Лукаса
    int
  u = 1,
  v = p,
  u2m = 1,
  v2m = p,
  qm = q,
  qm2 = q*2,
  qkd = q;
    for (unsigned bit = 1, bits = bits_in_number(d); bit < bits; bit++)
    {
  mulmod (u2m, v2m, n);
  mulmod (v2m, v2m, n);
  while (v2m < qm2)
    v2m += n;
  v2m -= qm2;
  mulmod (qm, qm, n);
  qm2 = qm;
  qm2 *= 2;
  if (test_bit (d, bit))
  {
    int t1, t2;
    t1 = u2m;
    mulmod (t1, v, n);
    t2 = v2m;
    mulmod (t2, u, n);
    
    int t3, t4;
    t3 = v2m;
    mulmod (t3, v, n);
    t4 = u2m;
    mulmod (t4, u, n);
    mulmod (t4, (int)dd, n);
 
    u = t1 + t2;
    if (!(u % 2 == 0))
    u += n;
    u /= 2;
    u %= n;
 
    v = t3 + t4;
    if (!(v % 2 == 0))
    v += n;
    v /= 2;
    v %= n;
    mulmod (qkd, qm, n);
  }
    }
 
    // точно простое (или псевдо-простое)
    if (u == 0 || v == 0)
  return true;
 
    // вычисляем оставшиеся члены
    int qkd2 = qkd;
    qkd2 *= 2;
    for (int r = 1; r < s; ++r)
    {
  mulmod (v, v, n);
  v -= qkd2;
  if (v < 0) v += n;
  if (v < 0) v += n;
  if (v >= n) v -= n;
  if (v >= n) v -= n;
  if (v == 0)
    return true;
  if (r < s-1)
  {
    mulmod (qkd, qkd, n);
    qkd2 = qkd;
    qkd2 *= 2;
  }
    }
 
    return false;
 
}
 
int main()
{
   
        int num=0;
        printf( "Enter any number -> \n");
        scanf("%u", &num);
        printf("%s\n",
            (
                lucas_selfridge (num) ?
                "It is possibly a prime number" :
                "It is composite number"
            )
        );
    
    return 0;
}
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
 
Текущее время: 23:30. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru