Форум программистов, компьютерный форум 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-пространстве) и движет их на зрителя. Демка на подобие того скринсейвера, который в Виндовс называется «сквозь вселенную». Программа работает, всё нормально. Но я хотел бы реализовать эту прогу с... подробнее

Показать сообщение отдельно
XuTPbIu_MuHTAu
Эксперт С++
2224 / 739 / 10
Регистрация: 27.05.2008
Сообщений: 1,508
06.08.2008, 15:20     Перевести код из С++ в Си
Поправил все ошибки. Через 5 минут выложу код

Добавлено через 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
#include <stdio.h>
#include <math.h>
 
 
 
//! Возвращает true, если n - точный квадрат простого числа
int perfect_square (const n)
{
int    sq = ceil (sqrt ((double)n));
    return sq*sq == n;
}
 
 
//! Возвращает количество бит в числе (т.е. минимальное количество бит, которыми можно представить данное число)
 
unsigned int bits_in_number (int n)
{
    
    unsigned int result = 0;
    if (n == 0)  return 1;
    while (n)
    {
  n /= 2;
  ++result;
    }
    return  result;
}
 
//! Возвращает значение k-го бита числа (биты нумеруются с нуля)
 
int test_bit (const n, unsigned int k)
{
    return (n & ((1) << k)) != 0;
}
 
//! Умножает a *= b (mod n)
 
void mulmod (int * a, int b, const int 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 int   n)
{
    int dd,d_abs,d_sign,g,p,q,n_1,s,d,u,v,u2m,v2m,qm,qm2,qkd,t1,t2,t3,t4,qkd2,r;
    unsigned int bit,bits;
    // сначала проверяем тривиальные случаи
    if (n == 2)
  return 1;
    if (n < 2 || (n % 2 == 0))
  return 0;
 
    // проверяем, что n не является точным квадратом, иначе алгоритм даст ошибку
    if (perfect_square (n))
  return 0;
 
    // алгоритм Селфриджа: находим первое число d такое, что:
    // jacobi(d,n)=-1 и оно принадлежит ряду { 5,-7,9,-11,13,... }
    
    for (d_abs = 5, d_sign = 1;; d_sign = -d_sign, d_abs+=2)
    {
  dd = d_abs * d_sign;
  g = gcd (n, d_abs);
  if (1 < g && g < n)
      // нашли делитель - d_abs
      return 0;
  if (jacobi ((dd), n) == -1)
      break;
    }
 
    // параметры Селфриджа
  p = 1,
  q = (p*p - dd) / 4;
    
    // разлагаем n+1 = d*2^s
    n_1 = n;
    ++n_1;
    //int s, d;
    transform_num (n_1, &s, &d);
 
    // алгоритм Лукаса
  u = 1,
  v = p,
  u2m = 1,
  v2m = p,
  qm = q,
  qm2 = q*2,
  qkd = q;
    for (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))
  {
      t1, t2;
      t1 = u2m;
      mulmod (&t1, v, n);
      t2 = v2m;
      mulmod (&t2, u, n);
      
      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 1;
 
    // вычисляем оставшиеся члены
    qkd2 = qkd;
    qkd2 *= 2;
    for (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 1;
  if (r < s-1)
  {
      mulmod (&qkd, qkd, n);
      qkd2 = qkd;
      qkd2 *= 2;
  }
    }
 
    return 0;
 
}
 
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;
}
Убрал ссылки.Инициализацию переменных перенес в начало. Убрал заголовочные файлы С++. Убрал конструкцию "++++"
Вроде всё. Спасибо,памятника не надо
 
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru