Develop and Download Open Source Software

Browse Subversion Repository

Contents of /branches/mty-makai/mty.c

Parent Directory Parent Directory | Revision Log Revision Log


Revision 115 - (show annotations) (download) (as text)
Tue Apr 10 10:30:12 2007 UTC (17 years ago) by notanpe
Original Path: trunk/naniya/mty/mty.c
File MIME type: text/x-csrc
File size: 20982 byte(s)
速度表示を変更。起動からの平均と区間の速度を表示。
1 /***********************************************************************
2 *
3 * file: mty.c
4 *
5 * まあ、待て屋。
6 *
7 * $Id$
8 *
9 */
10
11 #include <assert.h>
12 #include <ctype.h>
13 #include <errno.h>
14 #include <malloc.h>
15 #include <limits.h>
16 #include <stdarg.h>
17 #include <stddef.h>
18 #include <stdio.h>
19 #include <stdlib.h>
20 #include <string.h>
21 #include <time.h>
22
23 #if defined(WIN32)
24
25 #include <windows.h>
26 #include <sys/timeb.h>
27
28 #elif defined(__GNUC__)
29
30 #include <sys/time.h>
31
32 #endif
33
34 #include "config.h"
35 #include "cp932.h"
36 #include "crypt64.h"
37 #include "desconst.h"
38 #include "expr_parse.h"
39 #include "scoreboard.h"
40 #include "synth.h"
41 #include "tr64.h"
42 #include "translate.h"
43 #include "wdict.h"
44
45 #if USE_DT
46 #include "dt4.h"
47 #endif
48
49 /* 鍵文字列 */
50 unsigned char key[8 + 8];
51 unsigned char okey[8 + 8];
52
53 /* 鍵クラス */
54 static
55 struct
56 {
57 unsigned cls;
58 unsigned map[256];
59 } kcls[8 + 8];
60
61
62 /* 拡張鍵クラス */
63 #define KCLS_DT0 64
64 #define KCLS_DT1 128
65 #define KCLS_K2 256
66
67 #if USE_DT
68 /* 鍵キメ用辞書インデクス */
69 struct DT *kd[8 + 8];
70
71 /* 辞書インデクス */
72 struct DT *dtidx[0x100 + 1];
73 #endif
74
75 /* 指定されたクラスと入っているキーから、classify を行う */
76 void
77 key_make_map(int n)
78 {
79 int i, j;
80 unsigned c = kcls[n].map[key[n]];
81
82 #if USE_DT
83 if (3 <= n && n < 7 && kd[n - 3])
84 {
85 /* 辞書のケツの文字。後ろにナニヤラキャラクタが来る */
86 c = kd[n - 3]->c[0];
87 if ((0x81 <= c && c <= 0x9F)
88 || (0xE0 <= c && c <= 0xFC))
89 c = KCLS_K2;
90 else
91 c = (cp932[256 * key[n]]
92 | cp932[256 * (key[n] ^ 0x80)]);
93 #if DEBUG>=1
94 printf("*n=%d, key=%02X, cls=%04X\n",
95 n,
96 key[n],
97 c);
98 #endif
99 }
100 else if (2 <= n && n < 6 && kd[n - 2])
101 {
102 return;
103 }
104 else if (1 <= n && n < 5 && kd[n - 1])
105 {
106 return;
107 }
108 else if (1 <= n && n < 5 && !kd[n - 1]
109 //&& (c & KCLS_K2)
110 && (c & KCLS_DT1))
111 {
112 /* 漢字2文字を拾っていきまつ */
113 #if DEBUG>=1
114 printf("(%d)%02X %02X(%02X:%02X:%02X:%02X)\n",
115 n, key[n - 1], key[n],
116 cp932[(256 * key[n - 1] + key[n])],
117 cp932[(256 * key[n - 1] + key[n]) ^ 0x0080],
118 cp932[(256 * key[n - 1] + key[n]) ^ 0x8000],
119 cp932[(256 * key[n - 1] + key[n]) ^ 0x8080]);
120 #endif
121 if (n != 1 && n != 2
122 && (cp932[(256 * key[n - 1] + key[n]) ^ 0x0080] & KCLS_DT1))
123 key[n] ^= 0x80;
124 else if (n != 2 && n != 3
125 && (cp932[(256 * key[n - 1] + key[n]) ^ 0x8000] & KCLS_DT1))
126 key[n - 1] ^= 0x80;
127 else if (n > 3 && (cp932[(256 * key[n - 1] + key[n]) ^ 0x8080] & KCLS_DT1))
128 key[n - 1] ^= 0x80, key[n] ^= 0x80;
129 if (cp932[256 * key[n - 1] + key[n]] & KCLS_DT1)
130 {
131 for (kd[n - 1] = dtidx[key[n - 1]];
132 kd[n - 1]->c[1] != key[n];
133 kd[n - 1]++)
134 assert(kd[n - 1]->c[0] == key[n - 1]);
135 #if DEBUG>=1
136 printf("(%02X%02X:%02X%02X)%c%c%c%c\n",
137 kd[n - 1]->c[0],
138 kd[n - 1]->c[1],
139 kd[n - 1]->c[2],
140 kd[n - 1]->c[3],
141 kd[n - 1]->c[0],
142 kd[n - 1]->c[1],
143 kd[n - 1]->c[2],
144 kd[n - 1]->c[3]);
145 #endif
146 return;
147 }
148 }
149 else if (n < 4 && (c & KCLS_DT0) && kd[n] == NULL)
150 {
151 /* カタカナ埋め込みいきます */
152 assert(kd[n] == NULL);
153 #if DEBUG>=1
154 printf("n=%d, key=%02X\n", n, key[n]);
155 #endif
156 kd[n] = dtidx[key[n]];
157 if (!kd[n]
158 && !(n == 1 || n == 2)
159 && dtidx[key[n] ^ 0x80])
160 {
161 key[n] ^= 0x80;
162 kd[n] = dtidx[key[n]];
163 }
164 if (kd[n])
165 return;
166 }
167 else
168 {
169 kd[n] = NULL;
170 }
171 #endif
172
173 /* 最後の部分は class map を生成する必要ナシ */
174 if (n >= 6)
175 return;
176
177 for (i = 0; i < 256; i++)
178 {
179 unsigned bm = 0;
180 #if 1
181 if (c & KCLS_K1)
182 {
183 if (cp932[256 * key[n] + i] & KCLS_K1)
184 bm |= KCLS_K2 | (cp932[256 * key[n] + i] & KCLS_DT1);
185 if (cp932[256 * (key[n] ^ 0x80) + i] & KCLS_K1)
186 bm |= KCLS_K2 | (cp932[256 * (key[n] ^ 0x80) + i] & KCLS_DT1);
187 #if 0
188 bm |= ((cp932[256 * key[n] + i] & KCLS_K1)
189 || (cp932[256 * (key[n] ^ 0x80) + i] & KCLS_K1)
190 ? KCLS_K2 : 0);
191 #endif
192 }
193 if (c & (KCLS_AN | KCLS_KA | KCLS_K2))
194 for (j = 0; j < 256; j++)
195 {
196 bm |= cp932[256 * i + j] & (KCLS_AN | KCLS_KA | KCLS_K1
197 | KCLS_DT0);
198 #if 0
199 if (j >= 127 && !(n == 0 || n == 1))
200 break;
201 #endif
202 }
203 kcls[n + 1].map[i] = bm;
204 #endif
205 if (i >= 128 && !(n == 0 || n == 1))
206 kcls[n + 1].map[i - 128] |= kcls[n + 1].map[i];
207 }
208
209 if (n < 6)
210 kcls[n + 1].map[0x00] = kcls[n + 1].map[0x80] = 0;
211 if (n == 6)
212 kcls[7].map[0x00] |= KCLS_AN;
213 }
214
215 #if USE_DT
216 unsigned
217 dt_get(int kdn,
218 int xn,
219 int n,
220 int ch)
221 {
222 int i;
223 #if DEBUG>=1
224 printf("*dt_get(%d)%c%c%c%c(%02X%02X:%02X%02X)->ch=%d",
225 n,
226 kd[kdn]->c[0], kd[kdn]->c[1], kd[kdn]->c[2], kd[kdn]->c[3],
227 kd[kdn]->c[0], kd[kdn]->c[1], kd[kdn]->c[2], kd[kdn]->c[3],
228 ch);
229 #endif
230 /* まずは数える */
231 for (i = 0;
232 kd[kdn][i].c[xn] == kd[kdn]->c[xn];
233 i++)
234 ;
235 assert(i > 0);
236 kd[kdn] += ch % i;
237 #if DEBUG>=1
238 printf("/%d\n dt_get: %c%c%c%c(%02X%02X:%02X%02X)->ch=%d\n",
239 i,
240 kd[kdn]->c[0], kd[kdn]->c[1], kd[kdn]->c[2], kd[kdn]->c[3],
241 kd[kdn]->c[0], kd[kdn]->c[1], kd[kdn]->c[2], kd[kdn]->c[3],
242 ch);
243 #endif
244 return kd[kdn]->c[n];
245 }
246 #endif
247
248 /* マップから文字を拾ってセット */
249 unsigned
250 key_set(int n, unsigned ch)
251 {
252 int cnt = 0, i;
253
254 #if USE_DT
255 if (3 <= n && n < 7 && kd[n - 3])
256 {
257 return dt_get(n - 3, 2, 3, ch);
258 return kd[n - 3]->c[3];
259 }
260 else if (2 <= n && n < 6 && kd[n - 2])
261 {
262 return dt_get(n - 2, 1, 2, ch);
263 return kd[n - 2]->c[2];
264 }
265 else if (1 <= n && n < 5 && kd[n - 1])
266 {
267 return dt_get(n - 1, 0, 1, ch);
268 return kd[n - 1]->c[1];
269 }
270 #endif
271
272 #if DEBUG>=3
273 if (cnt == 0)
274 {
275 printf("n=%d, ch=%d, (n-1)=%02X\n", n, ch, key[n - 1]);
276 int j;
277 for (i = 0; i < 16; i++)
278 {
279 printf("map[0x%02X] =", 16 * i);
280 for (j = 0; j < 16; j++)
281 printf(" %03X", kcls[n].map[16 * i + j]);
282 printf("\n");
283 }
284 }
285 #endif
286 for (i = 0; i < 256; i++)
287 {
288 if (kcls[n].map[i])
289 {
290 if (ch-- == 0)
291 return i;
292 cnt++;
293 }
294 if (n != 1 && n != 2 && i >= 127)
295 break;
296 }
297 /* 見つからなかったのでもいっぺん */
298 assert(cnt > 0);
299 ch %= cnt;
300 for (i = 0; i < 256; i++)
301 if (kcls[n].map[i])
302 {
303 if (ch-- == 0)
304 return i;
305 }
306 assert(!"not matched");
307 return 0;
308 }
309
310 /* bitwise key をセット */
311 static
312 void
313 key_set64(struct KEY *key64,
314 int n,
315 unsigned k,
316 unsigned vk,
317 unsigned sk)
318 {
319 int i, j;
320 if (!((vk | sk) & 0x7F))
321 return;
322
323 for (i = 0; i < 7; i++)
324 {
325 if (n == 7 && i < N_STRIDE) continue;
326 if (sk & (1 << i))
327 {
328 /* セット */
329 int o = tr_pc1[n][6 - i] - 1;
330 if (o < 28)
331 {
332 assert(o >= 0);
333 for (j = 0; j < N_ALU; j++)
334 key64->k[0][0][o].a[j]
335 = key64->k[0][1][o].a[j]
336 = -!!(k & (1 << i));
337 }
338 else
339 {
340 assert(o >= 28);
341 assert(o < 56);
342 for (j = 0; j < N_ALU; j++)
343 key64->k[1][0][o - 28].a[j]
344 = key64->k[1][1][o - 28].a[j]
345 = -!!(k & (1 << i));
346 }
347 }
348 else if (vk & (1 << i))
349 {
350 /* 反転 */
351 int o = tr_pc1[n][6 - i] - 1;
352 if (o < 28)
353 {
354 assert(o >= 0);
355 for (j = 0; j < N_ALU; j++)
356 key64->k[0][0][o].a[j]
357 = key64->k[0][1][o].a[j]
358 = ~key64->k[0][0][o].a[j];
359 }
360 else
361 {
362 assert(o >= 28);
363 assert(o < 56);
364 for (j = 0; j < N_ALU; j++)
365 key64->k[1][0][o - 28].a[j]
366 = key64->k[1][1][o - 28].a[j]
367 = ~key64->k[1][0][o - 28].a[j];
368 }
369 }
370 }
371 }
372
373 /* 指定されたクラスの開始値にリセット
374 直前の文字のクラスに縛られる */
375 int
376 key_reset(int n)
377 {
378 if (n >= 8)
379 return 1;
380 if (n == 7)
381 {
382 key[7] = 0;
383 return 1;
384 }
385
386 /* 0-2 文字目はランダムに決める
387 3 文字目以降は初期値に */
388 if (n >= 3)
389 key[n] = key_set(n, 0);
390 else
391 key[n] = key_set(n, rand());
392
393 #if DEBUG>=3
394 printf("key[%d]=%02X ncls=%04X\n", n, key[n], kcls[n].map[key[n]]);
395 #endif
396
397 /* セットされた文字を元に、次キャラの文字クラスを決める */
398 key_make_map(n);
399
400 return key_reset(n + 1);
401 }
402
403 /* 指定された鍵空間の中で、キーをひとつ進める
404 安全にインクリメントできた場合 true を返す */
405 static
406 int
407 key_inc(int n)
408 {
409 if (n >= 8)
410 return 0;
411 else if (n == 7)
412 {
413 /* 最後のバイト */
414 key[7] = (key[7] + (1 << N_STRIDE)) & 0x7F;
415 if (key[7]) return 1;
416 else return 0;
417 }
418 else if (key_inc(n + 1)
419 /*
420 && key_inc(n + 1)
421 && key_inc(n + 1)
422 && key_inc(n + 1)*/
423 )
424 return 1;
425
426 /* Salt はインクリメントしない約束にする */
427 if (n == 1 || n == 2)
428 return 1;
429
430 #if DEBUG>=3
431 printf("key_inc(n=%d,ck=%02X)\n", n, key[n]);
432 #endif
433
434 #if USE_DT
435 /* 辞書語はインクリメントしていい約束にする */
436 if (3 <= n && n < 7 && kd[n - 3])
437 {
438 if ((key[n - 3] & 0x7F) == ((kd[n - 3] + 1)->c[0] & 0x7F)
439 && (key[n - 2] & 0x7F) == ((kd[n - 3] + 1)->c[1] & 0x7F)
440 && (key[n - 1] & 0x7F) == ((kd[n - 3] + 1)->c[2] & 0x7F))
441 {
442 memcpy(&key[n - 3], &(++kd[n - 3])->c[0], 4);
443 #if DEBUG>=2
444 printf(">dt_get:%c%c%c%c(%02X%02X:%02X%02X)\n",
445 kd[n - 3]->c[0], kd[n - 3]->c[1], kd[n - 3]->c[2], kd[n - 3]->c[3],
446 kd[n - 3]->c[0], kd[n - 3]->c[1], kd[n - 3]->c[2], kd[n - 3]->c[3]);
447 #endif
448 return 1;
449 }
450 else
451 {
452 return 0;
453 }
454 }
455 else if (2 <= n && n < 6 && kd[n - 2])
456 {
457 if ((key[n - 2] & 0x7F) == ((kd[n - 2] + 1)->c[0] & 0x7F)
458 && (key[n - 1] & 0x7F) == ((kd[n - 2] + 1)->c[1] & 0x7F))
459 {
460 memcpy(&key[n - 2], &(++kd[n - 2])->c[0], 4);
461 #if DEBUG>=2
462 printf(">dt_get:%c%c%c%c(%02X%02X:%02X%02X)\n",
463 kd[n - 2]->c[0], kd[n - 2]->c[1], kd[n - 2]->c[2], kd[n - 2]->c[3],
464 kd[n - 2]->c[0], kd[n - 2]->c[1], kd[n - 2]->c[2], kd[n - 2]->c[3]);
465 #endif
466 return 1;
467 }
468 else
469 {
470 return 0;
471 }
472 if (kd[n - 2]->c[0] == key[n - 2])
473 return 1;
474 else
475 return 0;
476 }
477 else if (1 <= n && n < 5 && kd[n - 1])
478 {
479 unsigned c2 = kd[n - 1]->c[0];
480 if ((0x81 <= c2 && c2 <= 0x9F)
481 || (0xE0 <= c2 && c2 <= 0xFC))
482 {
483 kd[n - 1] = NULL;
484 #if 0
485 if (!(n == 1 && n == 2))
486 key[n] &= 0x7F;
487 if (!(n == 2 && n == 3))
488 key[n - 1] &= 0x7F;
489 #endif
490 key_make_map(n - 1);
491 }
492 else if ((key[n - 1] & 0x7F) == ((kd[n - 1] + 1)->c[0] & 0x7F))
493 {
494 memcpy(&key[n - 1], &(++kd[n - 1])->c[0], 4);
495 #if DEBUG>=2
496 printf(">dt_get:%c%c%c%c(%02X%02X:%02X%02X)\n",
497 kd[n - 1]->c[0], kd[n - 1]->c[1], kd[n - 1]->c[2], kd[n - 1]->c[3],
498 kd[n - 1]->c[0], kd[n - 1]->c[1], kd[n - 1]->c[2], kd[n - 1]->c[3]);
499 #endif
500 return 1;
501 }
502 else
503 {
504 return 0;
505 }
506 #if 0
507 if (kd[n - 1]->c[0] == key[n - 1])
508 return 1;
509 else
510 return 0;
511 #endif
512 }
513 else if (n < 4 && kd[n])
514 {
515 if (0 && kd[n]->c[0] == key[n])
516 return 1;
517 kd[n] = NULL;
518 #if 0
519 if (!(n == 1 || n == 2))
520 key[n] &= 0x7F;
521 #endif
522 }
523 #endif
524
525 /* 実際に増やしてみる */
526 assert(n >= 3);
527 for (;;)
528 {
529 if (n <= 3
530 && !(key[n] & 0x80)
531 && kcls[n].map[key[n] ^ 0x80] & (KCLS_DT0))
532 {
533 /* 半角カタカナの1バイト目 */
534 key[n] ^= 0x80;
535 }
536 else
537 {
538 key[n] = (key[n] & 0x7F) + 1;
539 if (key[n] >= 0x80)
540 return 0;
541 }
542
543 if (kcls[n].map[key[n]])
544 {
545 key_make_map(n);
546 key_reset(n + 1);
547 return 1;
548 }
549 }
550 while (++key[n] < 0x80)
551 {
552 if (kcls[n].map[key[n]])
553 {
554 key_make_map(n);
555 key_reset(n + 1);
556 return 1;
557 }
558 }
559 return 0;
560 }
561
562 /* 鍵を完全にリセットする
563 Saltもセットし直す */
564 static
565 void
566 key_init()
567 {
568 int i, j;
569
570 #if USE_DT
571 /* 辞書を、インデクスを作りながらマップにはめこんで逝く
572 辞書はコード順昇順に並んでるものとする */
573 for (i = 0; i < dtcnt; i++)
574 {
575 unsigned c = dt[i].c[0];
576
577 assert(dt[i].c[0]
578 && dt[i].c[1]
579 && dt[i].c[2]
580 && dt[i].c[3]);
581
582 /* BSD 鯖でしにそうな文字は残念ながら除外 */
583 assert((dt[i].c[0] & 0x7F)
584 && (dt[i].c[1] & 0x7F)
585 && (dt[i].c[2] & 0x7F)
586 && (dt[i].c[3] & 0x7F));
587
588 /* インデクス */
589 if (!dtidx[c])
590 dtidx[c] = &dt[i];
591
592 if ((0x81 <= c && c <= 0x9F)
593 || (0xE0 <= c && c <= 0xFC))
594 {
595 /* 全角なので、2バイトきまった時点で立てる */
596 cp932[256 * c + dt[i].c[1]] |= KCLS_DT1;
597 }
598 else if (0xA1 <= c && c <= 0xDF)
599 {
600 /* 半角カナ */
601 for (j = 0; j < 256; j++)
602 cp932[256 * c + j] |= KCLS_DT0;
603 }
604 }
605 /* ケツ、ちうか番人 */
606 dtidx[0x100] = &dt[i];
607 #endif
608
609 key[8] = 0;
610
611 /* 初期マップを組む */
612 for (i = 0; i < 256; i++)
613 {
614 unsigned bm = 0;
615 kcls[0].map[i] = 0;
616 for (j = 0; j < 256; j++)
617 bm |= cp932[256 * i + j];
618 kcls[0].map[i] = bm & (KCLS_AN | KCLS_KA | KCLS_K1
619 | KCLS_DT0
620 );
621 if (i >= 128)
622 kcls[0].map[i - 128] |= kcls[0].map[i];
623 }
624
625 key_reset(0);
626 }
627
628 /***************************************************************
629 *
630 * Salt のセット
631 * オペランドのオフセットを書き換えて回ってるので注意
632 *
633 */
634
635 #if N_STRIDE == 6
636 #define C(c,i,j,o) (*(int8_t *)((c) + (loo - crypt64_sta) + los[6 * (i) + (j) + (o)]))
637 #elif N_STRIDE == 7
638 #define C(c,i,j,o) (*(int32_t *)((c) + (loo - crypt64_sta) + los[6 * (i) + (j) + (o)]))
639 #endif
640
641 void
642 set_salt(signed char *code,
643 unsigned char const *k)
644 {
645 int i, j;
646
647 for (i = 0; i < 2; i++)
648 {
649 unsigned s = k[1 + i] & 255;
650 if (s > 'z')
651 s = 0;
652 else if (s >= 'a')
653 s = s - 'a' + 2 + 10 + 26;
654 else if (s >= 'A')
655 s = s - 'A' + 2 + 10;
656 else if (s >= '.')
657 s = s - '.';
658 else
659 s = 0;
660
661 #if DEBUG>=1
662 printf("Salt %d:%3o\n", i, s & 63);
663 #endif
664 for (j = 0; j < 6; j++)
665 {
666 #if DEBUG>=2
667 //printf("Salt %d:%d %+3d:%+3d",
668 printf("Salt %d:%d %08lX:%08lX",
669 i, j,
670 C(code, i, j, 0),
671 C(code, i, j, 24));
672 #endif
673 if (s & (1 << j))
674 {
675 C(code, i, j, 0) = sizeof(WS_T) * (((4 * i + j + 15) & 31) - 16);
676 C(code, i, j, 24) = sizeof(WS_T) * (((4 * i + j - 1) & 31) - 16);
677 }
678 else
679 {
680 C(code, i, j, 0) = sizeof(WS_T) * (((4 * i + j - 1) & 31) - 16);
681 C(code, i, j, 24) = sizeof(WS_T) * (((4 * i + j + 15) & 31) - 16);
682 }
683 C(code, i, j, 12) = sizeof(WS_T) * (((4 * i + j + 7) & 31) - 16);
684 C(code, i, j, 36) = sizeof(WS_T) * (((4 * i + j + 23) & 31) - 16);
685 #if DEBUG>=2
686 //printf(" => %+3d:%+3d\n",
687 printf(" => %08lX:%08lX\n",
688 C(code, i, j, 0),
689 C(code, i, j, 24));
690 #endif
691 }
692 }
693 }
694
695 static
696 int
697 log_printf(FILE *ofp, char const *fmt, ...)
698 {
699 int r;
700 va_list ap;
701 va_start(ap, fmt);
702 vfprintf(stdout, fmt, ap);
703 r = vfprintf(ofp, fmt, ap);
704 va_end(ap);
705 if (r > 0)
706 return r;
707 perror("log_printf");
708 exit(errno);
709 }
710
711 /***************************************************************
712 *
713 * メインループとか
714 *
715 */
716
717 /* 定数項 */
718 #if N_STRIDE == 7
719 static SLICE const sk6[N_STRIDE] =
720 {
721 {0xAAAAAAAAUL, 0xAAAAAAAAUL, 0xAAAAAAAAUL, 0xAAAAAAAAUL},
722 {0xCCCCCCCCUL, 0xCCCCCCCCUL, 0xCCCCCCCCUL, 0xCCCCCCCCUL},
723 {0xF0F0F0F0UL, 0xF0F0F0F0UL, 0xF0F0F0F0UL, 0xF0F0F0F0UL},
724 {0xFF00FF00UL, 0xFF00FF00UL, 0xFF00FF00UL, 0xFF00FF00UL},
725 {0xFFFF0000UL, 0xFFFF0000UL, 0xFFFF0000UL, 0xFFFF0000UL},
726 {0x00000000UL, 0xFFFFFFFFUL, 0x00000000UL, 0xFFFFFFFFUL},
727 {0x00000000UL, 0x00000000UL, 0xFFFFFFFFUL, 0xFFFFFFFFUL},
728 };
729 #elif N_STRIDE == 6
730 static SLICE const sk6[N_STRIDE] =
731 {
732 {0xAAAAAAAAUL, 0xAAAAAAAAUL},
733 {0xCCCCCCCCUL, 0xCCCCCCCCUL},
734 {0xF0F0F0F0UL, 0xF0F0F0F0UL},
735 {0xFF00FF00UL, 0xFF00FF00UL},
736 {0xFFFF0000UL, 0xFFFF0000UL},
737 {0x00000000UL, 0xFFFFFFFFUL},
738 };
739 #endif
740
741 ALIGN_PREFIX(16) struct KEY key64 ALIGN_SUFFIX(16);
742 ALIGN_PREFIX(16) struct PARAM param64 ALIGN_SUFFIX(16);
743
744 int
745 main(int argc, char *argv[])
746 {
747 int i;
748 int mincnt;
749 int nblk_hit, nblk_total;
750 int nap_hit, nap_total;
751 signed char *code = NULL;
752 FILE *ofp;
753 FILE *sfp; /* scoreboard */
754 struct ITREE *root_expr;
755 int cr;
756
757 int xhash_loaded;
758
759 #define LOOP_FACTOR 100000
760 #define UPDATE_INTERVAL 5 /* ステータス表示の最低間隔 */
761 struct status {
762 time_t startTime; /* 開始時刻 */
763 time_t lastTime; /* 最後に表示した時刻 */
764 unsigned long lastLoopCnt; /* 最後に表示した時の検索数 */
765 unsigned long loop; /* ループ回数 % LOOP_FACTOR */
766 unsigned long mloop; /* ループ回数 / LOOP_FACTOR */
767 } status;
768
769 #if 0
770 if (argc < 2)
771 {
772 fprintf(stderr, "式きぼんぬ\n");
773 return 1;
774 }
775 #endif
776
777 assert((1 << N_STRIDE) == N_ALU * ALU_BITS);
778
779 /* t[16] は、内部演算で使用する、all 1 が入っている */
780 for (i = 0; i < N_ALU; i++)
781 param64.t[T_INV].a[i] = -1;
782
783 /* 固定キーのコピー */
784 for (i = 0; i < N_STRIDE; i++)
785 {
786 int o = tr_pc1[7][6 - i] - 1;
787 if (o < 28)
788 {
789 key64.k[0][0][o] = key64.k[0][1][o] = sk6[i];
790 }
791 else
792 {
793 o -= 28;
794 key64.k[1][0][o] = key64.k[1][1][o] = sk6[i];
795 }
796 }
797
798 /* キースケジュールをここに押し込めておく
799 従来は crypt64.S 内で完結するように引いていた */
800 for (i = 0; i < 28; i++)
801 key64.ks[i].a[0] = sizeof(WS_T) * ks_ls[i];
802
803 /* タゲ読み込み */
804 root_expr = expr_parse("target.txt");
805
806 /* コードを生成・展開 */
807 sfp = scoreboard_open();
808 fwrite(crypt64_sta, 1, crypt64_end - crypt64_sta, sfp); /* prologue */
809 synth_synthesize(sfp, root_expr);
810 fwrite(crypt64_ep, 1, crypt64_ep_end - crypt64_ep, sfp); /* epilogue */
811
812 /* コードをメモリに貼り付ける */
813 code = scoreboard_map(sfp);
814
815 /* キーの初期化 */
816 srand(time(NULL));
817 key_init();
818 set_salt(code, key);
819 for (i = 0; i < 8; i++)
820 key_set64(&key64, i, okey[i] = key[i], 0, 0x7F);
821
822 #if DEBUG>=1
823 printf("key=%p param=%p\n", &key64, &param64);
824 #endif
825 assert(!((ptrdiff_t)&key64 & (sizeof(WS_T) - 1)));
826 assert(!((ptrdiff_t)&param64 & (sizeof(WS_T) - 1)));
827
828 if ((ofp = fopen("log.txt", "at")) == NULL)
829 {
830 perror("log.txt");
831 return errno;
832 }
833
834 setvbuf(ofp, NULL, _IONBF, BUFSIZ); /* XXX MSVCRT では _IOLBF が期待通りに動作しない */
835
836 mincnt = 0x7FFFFFFF;
837 nblk_hit = nblk_total = 0;
838 nap_hit = nap_total = 0;
839 cr = 0;
840 status.startTime = status.lastTime = time( NULL );
841 status.lastLoopCnt = 0;
842 status.loop = 0;
843 status.mloop = 0;
844 /* 探索ループだぞっと */
845 for (;;)
846 {
847 int32_t cnt;
848 int k, kk;
849
850 /* 鍵のセット */
851 for (i = 0; i < 8; i++)
852 {
853 key_set64(&key64, i, key[i], key[i] ^ okey[i], 0);
854 okey[i] = key[i];
855 }
856
857 /* 呼ぶ!
858 LR 初期化は、サブモジュール内で行うべし
859 FASTCALL に準じた呼び出しのため、
860 ホントはいろいろレジスタが破壊されるハズ…なんだが。 */
861 cnt = CALL_CRYPT64(code, key64.k, param64.lr);
862
863 #if DEBUG>=1
864 if (mincnt > cnt && cnt > 0)
865 {
866 mincnt = cnt;
867 if (cr)
868 fprintf(stderr, "\n");
869 cr = 0;
870 fprintf(stderr, "cycle=%d\n", (int)cnt);
871 }
872 #endif
873
874 /* ヒットしたときの処理 */
875 xhash_loaded = 0;
876 for (kk = 0; kk < N_ALU; kk++)
877 {
878 ALU_T t;
879 if (!(kk & (N_ALU / N_Q - 1)))
880 nblk_total++, xhash_loaded = 0;
881
882 t = param64.t[HIT_ANY].a[kk];
883 if (!t)
884 continue;
885
886 nap_total += ALU_BITS;
887
888 for (k = 0; k < ALU_BITS; k++)
889 {
890 static uint64_t xhash[64];
891 char hash[16];
892 unsigned char buf[32];
893
894 if (!(t & ((ALU_T)1 << k)))
895 continue;
896
897 nap_hit++;
898
899 /* 転置 */
900 if (!xhash_loaded)
901 {
902 nblk_hit++;
903 CALL_TR64(&param64.lr[0][0].q[kk / (N_ALU / N_Q)], xhash);
904 xhash_loaded = 1;
905 }
906
907 /* 辞書を調べる */
908 if (!((param64.t[HIT_BOOL].a[kk] & ((ALU_T)1 << k))
909 || wdict_ishit(param64.t,
910 kk, k,
911 xhash[(ALU_BITS * kk + k) & 0x3F])))
912 continue;
913
914 for (i = 1; i < 11; i++)
915 {
916 unsigned c = 0;
917 c = (xhash[(ALU_BITS * kk + k) & 63] >> (6 * (i - 1))) & 0x3F; /* XXX */
918 hash[i - 1] = C64[c];
919 }
920 hash[10] = 0;
921
922 memcpy(buf, key, 32);
923 buf[8] = buf[9] = 0;
924 buf[7] = (buf[7] & -(1 << N_STRIDE) & 0x7F) + ALU_BITS * kk + k;
925 if (translate(buf, 0, 1))
926 {
927 if (cr)
928 fprintf(stderr, "\n");
929 cr = 0;
930 log_printf(ofp,
931 "◆%s #%-10.10s(%02X %02X %02X %02X %02X %02X %02X %02X/%02X)\n",
932 hash,
933 buf,
934 buf[0], buf[1], buf[2], buf[3],
935 buf[4], buf[5], buf[6], buf[7],
936 buf[8]);
937 }
938 else
939 {
940 #if DEBUG>=1
941 if (cr)
942 fprintf(stderr, "\n");
943 cr = 0;
944 log_printf(ofp,
945 "◆%s (%02X %02X %02X %02X %02X %02X %02X %02X )\n",
946 hash,
947 buf[0], buf[1], buf[2], buf[3],
948 buf[4], buf[5], buf[6], buf[7]);
949 #endif
950 }
951 }
952 }
953
954 status.loop += N_ALU * ALU_BITS;
955 if ( status.loop >= LOOP_FACTOR )
956 {
957 time_t curTime;
958 status.mloop += ( status.loop / LOOP_FACTOR );
959 status.loop %= LOOP_FACTOR;
960 curTime = time( NULL );
961 if ( curTime >= (status.lastTime + (time_t)UPDATE_INTERVAL) )
962 {
963 int a, b;
964 unsigned long loopCnt, zoneCnt;
965 loopCnt = status.mloop * LOOP_FACTOR + status.loop;
966 zoneCnt = loopCnt - status.lastLoopCnt;
967 a = loopCnt / (curTime - status.startTime);
968 b = zoneCnt / (curTime - status.lastTime);
969 #if DEBUG>=1
970 fprintf(stderr,
971 "%5d/%5d(%3d%%)",
972 nblk_hit, nblk_total, 100 * nblk_hit / nblk_total);
973 nblk_hit = nblk_total = 0;
974 if (nap_total)
975 fprintf(stderr,
976 " %5d/%5d(%3d%%)",
977 nap_hit, nap_total, 100 * nap_hit / nap_total);
978 else
979 fprintf(stderr,
980 " -----/-----(---%%)");
981 nap_hit = nap_total = 0;
982 #endif
983 fprintf( stderr,
984 "%8d.%03dktrips/s [%8d.%03dktrips/s]\r",
985 a / 1000, a % 1000,
986 b / 1000, b % 1000 );
987 status.lastTime = curTime;
988 status.lastLoopCnt = loopCnt;
989 cr++;
990 }
991 }
992 #if 1
993 if (!key_inc(3))
994 {
995 #if DEBUG>=2
996 printf("********************************\n");
997 #endif
998 key_reset(0);
999 set_salt(code, key);
1000 for (i = 0; i < 8; i++)
1001 key_set64(&key64, i, okey[i] = key[i], 0, 0x7F);
1002 }
1003 #endif
1004 }
1005
1006 return 0;
1007 }
1008
1009 /*
1010 * Local Variables:
1011 * tab-width: 4
1012 * End:
1013 *
1014 * EOF */

Properties

Name Value
svn:eol-style native
svn:keywords Author Date Id Rev URL

Back to OSDN">Back to OSDN
ViewVC Help
Powered by ViewVC 1.1.26