#include "unif01.h" #include "util.h" #include "addstr.h" #include "gdef.h" #include "num.h" #include "bbattery.h" #include <stdlib.h> #include <stdio.h> #include <string.h> #include <math.h> typedef unsigned long long b64; typedef struct { int p, k, m1, m2, d; b64 *seed; double minv; }md_para; typedef struct { int mdi,mda,mdb,mdc,mdj; }md_index; (中刪) static unsigned long mdx_Bits(void *vpar,void *vind) { return (unsigned long)(4294967296*mdx_pro(vpar,vind)); } (中刪) void Chris_DeleteMdx(unif01_Gen *gen) { md_para *index; if(NULL==gen)return; index = gen->state; util_Free(para->seed); gen->state = util_Free(gen->state); gen->param = util_Free(gen->param); gen->name = util_Free(gen->name); util_Free(gen); } ============================================== cygwin環境 compile之後，出現的錯誤訊息 warning: integer constant is too large for "long" type error 'para' undeclcared (first use in this function) (Each undeclared identifier is reported only once for each function it appears in.) 請問一下，這兩段的問題出在哪?? bug一直de不掉，仿另一篇，但原始的不會有第2項錯誤.... 感恩 --

※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 123.195.137.179 ※ 編輯: chrisjon 來自: 123.195.137.179 (06/11 00:38) 這我之前有問過，有兩種宣告方式 long long和int64_t 不過聽說int64_t要新版的才有支援? 怪了...通常計算後的最大值是2^31-1 而mdx_pro產出來的值是X/p，p為2^31-1，而X範圍在0~p-1 (也就是X/p是介於0<=x<1之間的小數，也就是"不可能出現1") 照理用unsigned int 應該可用數是0~2^32-1...也就是上面說的4294967296-1 (2^32)*x，不可能會出現2^32(因為x一定是小於1的小數) 我把mdx_pro的程式碼貼回來好了... (因為指標太多，看起來很花，問完就刪掉了^^"把括號加上顏色應該比較好看範圍) 以下接在上程式碼中，紅色"中刪"的位置 #define DMOD(n,p)((n)%(p)) static double mdx_pro(void *vpar,void *vind) { md_para *para = vpar; md_index *index = vind; if(++index->mdi==para->k){index->mdi=0;} //計算第i-1項指標 index->mdj=(index->mdi+1)%para->k; //計算計算第i-k-1項指標 index->mda=(index->mdi+2)%para->k; //計算計算第i-k項指標 if(index->mdi==0) //計算第i-2項指標 {index->mdb=para->k-1;} else{index->mdb=index->mdi-1;} para->seed[index->mdj] = DMOD(DMOD(para->seed[index->mdi] + para->m1*(para->seed[index->mdi] + para->p - para->seed[index->mdj]),para->p)+ DMOD(para->d*(para->seed[index->mda]+para->p-para->seed[index->mdb]),para->p) ,para->p); return ((double)para->seed[index->mdj]+0.5)*para->minv; } 麻煩各位先進再幫忙看了，感謝 ※ 編輯: chrisjon 來自: 123.195.137.179 (06/11 12:11) 不好意思....趕著上課，沒寫清楚.. 2^32 = 42949672956 ※ 編輯: chrisjon 來自: 120.125.72.28 (06/11 16:07)