編集履歴

質問編集履歴

p,q,Nが検証フェーズ前では正しい値が与えられたことを追加しました．

2021/01/29 06:38

投稿

Q3fdxrGzWzu0u5n

スコア8

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -6,6 +6,8 @@
 コードの下に実行結果を載せました．
+検証フェーズの直前で，p,q,Nの値を確認したところ正しい値が出力されました．
 ```C

出力結果を更新しました．

2021/01/29 06:38

投稿

Q3fdxrGzWzu0u5n

スコア8

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -568,6 +568,58 @@
 ```出力結果
+initialize.c:59:15: warning: implicit declaration of function 'pow1' is invalid in C99 [-Wimplicit-function-declaration]
+      R[i][j]=pow1(g,s1[i][j],N);
+              ^
+initialize.c:89:29: warning: implicit declaration of function 'phi' is invalid in C99 [-Wimplicit-function-declaration]
+  printf("Φ(N)=%d, f=%d\n",phi(N),f);
+                           ^
+initialize.c:91:10: warning: implicit declaration of function 'pow1' is invalid in C99 [-Wimplicit-function-declaration]
+  int R0=pow1(g,s0,N);
+         ^
+initialize.c:112:10: warning: implicit declaration of function 'jacobi' is invalid in C99 [-Wimplicit-function-declaration]
+      if(jacobi(K[i],P)==1){
+         ^
+initialize.c:130:24: warning: implicit declaration of function 'h' is invalid in C99 [-Wimplicit-function-declaration]
+          r[i][j]=fmod(h(I[i][j],R1[i],t[i],P),P);
+                       ^
+initialize.c:142:20: warning: implicit declaration of function 'h' is invalid in C99 [-Wimplicit-function-declaration]
+      y[i][j]=fmod(h(I[i][j],S[i],t[i],P),P);
+                   ^
+initialize.c:161:45: warning: more '%' conversions than data arguments [-Wformat]
+      printf("R0=%d, s[%d][%d]=%d, N=%d, p=%d, q=%d\n",R0,s1[i][j],N,p,q);
+                                           ~^
+initialize.c:163:46: warning: format specifies type 'int' but the argument has type 'double' [-Wformat]
+      printf("I'[%d][%d]=%d,f=%d,N=%d\n",i,j,I1[i][j],f,N);
+                         ~~                  ^~~~~~~~
+                         %f
+8 warnings generated.
 m=3
 t[0]=2, n[0]=2

どこで宣言して，どこがおかしくなったのか具体的な行数を追加し，出力結果を載せました．

2021/01/29 06:35

投稿

Q3fdxrGzWzu0u5n

スコア8

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -2,6 +2,12 @@
 なお，以前同じ内容で質問してベストアンサーを決めて解決済みになってしまったため，もう一度質問し直しました．また，以前ご指摘のあった，領域の問題は解決したと思います．
+40,41行目でp,q,Nを宣言して，151行目以降の検証フェーズにて，値がおかしくなっています．
+コードの下に実行結果を載せました．
 ```C
 //逐次秘密分散法の実装
@@ -318,7 +324,13 @@
       //I'ijを計算する
+//N,p,qの値がおかしい
-      printf("R0=%d, s[%d][%d]=%d, N=%d\n",R0,s1[i][j],N);
+      printf("R0=%d, s[%d][%d]=%d, N=%d, p=%d, q=%d\n",R0,s1[i][j],N,p,q);
       I1[i][j]=pow1(R0,s1[i][j],N);
@@ -326,12 +338,6 @@
       //Rij==I'ij^f mod Nを検証する
-ここでN,p,qの値がおかしくなっている．
       int h=pow1(I1[i][j],f,N);
       printf("I'[%d][%d]^f mod N=%d\n",i,j,h);
@@ -352,9 +358,63 @@
   if(countv==0)printf("全て正しいシェアでした\n");
+  /*
+      //復元フェーズ
+    //プレイヤー
+  double S1[m],r[m]; //求めるSi
+  double y1[m][P];
+  for(i=1;i<=m;i++){
+    S1[i]=0;
+    for(j=1;j<=t[i];j++){
+      //ラグランジュ補間を使用して，Siを求める
+      y1[i][j]=fmod(h(I[i][j],S[i],t[i],P),P);
+      printf("y[%d][%d]=%lf, ",i,j,y1[i][j]);
+      r[i]=y1[i][j];
+      for(k=1;k<=t[i];k++){
+        if(j!=k){
+          r[i]*=(0-I1[i][k])/(I1[i][j]-I1[i][k]);
+        }
+      }
+      S1[i]+=r[i];
+      printf("S[%d]=%lf ",i,S1[i]);
+      S1[i]=fmod(S1[i],P);
+    }
+    S1[i]=round(S1[i]);
+    printf("S[%d]=%lf\n",i,S1[i]);
+  }
+*/
   return 0;
@@ -502,8 +562,168 @@
 }
 ```
+```出力結果
+m=3
+t[0]=2, n[0]=2
+t[1]=2, n[1]=3
+t[2]=3, n[2]=4
+P=23
+s[0]=3, s[1]=5, s[2]=5,
+p=7, q=3, N=21
+g=5
+s[0][0]=1, s[0][1]=2,
+s[1][0]=1, s[1][1]=2, s[1][2]=3,
+s[2][0]=1, s[2][1]=2, s[2][2]=3, s[2][3]=4,
+R[0][0]=5, R[0][1]=4,
+R[1][0]=5, R[1][1]=4, R[1][2]=20,
+R[2][0]=5, R[2][1]=4, R[2][2]=20, R[2][3]=16,
+レベル1,2,3,で一意である．
+s0=5
+Φ(N)=12, f=17
+R0=17
+I[0][0]=17, I[0][1]=16,
+y[0][0]=20, y[0][1]=19,
+S[0]=3, s'[0]=10
+I[1][0]=17, I[1][1]=16, I[1][2]=20,
+K[1]=4
+K[1]は平方剰余である
+S[1]=2
+y[1][0]=19, y[1][1]=18, y[1][2]=22,
+S[1]=2, s'[1]=2
+I[2][0]=17, I[2][1]=16, I[2][2]=20, I[2][3]=4,
+K[2]=10
+K[2]は平方剰余でない
+R[2]=5
+r[2][0]=2, r[2][1]=4, r[2][2]=20, r[2][3]=18,
+y[2][0]=22, y[2][1]=1, y[2][2]=17, y[2][3]=15,
+S[2]=2
+R0=17, s[1][21]=7, N=3, p=23, q=0
+I'[0][0]=17,f=21,N=0
+I'[0][0]^f mod N=-2147483648
+ R[0][0]=5
+レベル0の0番目のプレイヤーは，誤ったシェアを提出しました
+R0=17, s[2][21]=7, N=3, p=0, q=0
+I'[0][1]=17,f=21,N=0
+I'[0][1]^f mod N=-2147483648
+ R[0][1]=4
+レベル0の1番目のプレイヤーは，誤ったシェアを提出しました
+R0=17, s[1][21]=7, N=3, p=0, q=0
+I'[1][0]=17,f=21,N=0
+I'[1][0]^f mod N=-2147483648
+ R[1][0]=5
+レベル1の0番目のプレイヤーは，誤ったシェアを提出しました
+R0=17, s[2][21]=7, N=3, p=0, q=0
+I'[1][1]=17,f=21,N=0
+I'[1][1]^f mod N=-2147483648
+ R[1][1]=4
+レベル1の1番目のプレイヤーは，誤ったシェアを提出しました
+R0=17, s[3][21]=7, N=3, p=0, q=0
+I'[1][2]=17,f=21,N=0
+I'[1][2]^f mod N=-2147483648
+ R[1][2]=20
+レベル1の2番目のプレイヤーは，誤ったシェアを提出しました
+R0=17, s[1][21]=7, N=3, p=0, q=0
+I'[2][0]=17,f=21,N=0
+I'[2][0]^f mod N=-2147483648
+ R[2][0]=5
+レベル2の0番目のプレイヤーは，誤ったシェアを提出しました
+R0=17, s[2][21]=7, N=3, p=0, q=0
+I'[2][1]=17,f=21,N=0
+I'[2][1]^f mod N=-2147483648
+ R[2][1]=4
+レベル2の1番目のプレイヤーは，誤ったシェアを提出しました
+R0=17, s[3][21]=7, N=3, p=0, q=0
+I'[2][2]=17,f=21,N=0
+I'[2][2]^f mod N=-2147483648
+ R[2][2]=20
+レベル2の2番目のプレイヤーは，誤ったシェアを提出しました
+R0=17, s[4][21]=7, N=3, p=0, q=0
+I'[2][3]=17,f=21,N=0
+I'[2][3]^f mod N=-2147483648
+ R[2][3]=16
+レベル2の3番目のプレイヤーは，誤ったシェアを提出しました
+```