diff options
-rw-r--r-- | gp-scripts/m+1stPrice | 35 |
1 files changed, 17 insertions, 18 deletions
diff --git a/gp-scripts/m+1stPrice b/gp-scripts/m+1stPrice index 8407e96..a6ac17c 100644 --- a/gp-scripts/m+1stPrice +++ b/gp-scripts/m+1stPrice | |||
@@ -5,6 +5,8 @@ | |||
5 | \\ Adapt the following values to your needs | 5 | \\ Adapt the following values to your needs |
6 | \\\\\\\\\\\\ | 6 | \\\\\\\\\\\\ |
7 | 7 | ||
8 | \\ auction parameter | ||
9 | M = 1 | ||
8 | \\ amount of bidders | 10 | \\ amount of bidders |
9 | n = 2^3 | 11 | n = 2^3 |
10 | \\ amount of possible prices | 12 | \\ amount of possible prices |
@@ -16,20 +18,18 @@ bid = vector(n,i,random(k)+1) | |||
16 | \\bid = vector(n,i,(i+1)%2) \\ second bidder wins (with ties) | 18 | \\bid = vector(n,i,(i+1)%2) \\ second bidder wins (with ties) |
17 | 19 | ||
18 | \\ prime finite field setup (result may be ambiguous if your prime is too small, 4*n*k seems to work fine) | 20 | \\ prime finite field setup (result may be ambiguous if your prime is too small, 4*n*k seems to work fine) |
19 | q = prime(4*n*k) | 21 | \\p = 263 |
20 | 22 | \\q = (p-1)/2 | |
21 | \\\\\\\\\\\\ | 23 | \\ use save prime p |
22 | \\ SETUP | 24 | q = prime(2^12) |
23 | \\\\\\\\\\\\ | 25 | p = 2*q + 1 |
24 | 26 | ||
25 | \\ p not needed? wat? | 27 | \\ get generator / primitive element for G_q |
26 | \\p = 47 | 28 | \\var = 'x \\ copy pasta from internet |
27 | 29 | \\pe=ffgen(minpoly(ffprimroot(ffgen(ffinit(p,1))),var),var) \\ get primitive element | |
28 | \\ get generator / primitive element for Z_q | 30 | \\1/(fforder(pe) == p-1) \\ error out, if ord(pe) is wrong |
29 | var = 'x \\ copy pasta from internet | 31 | \\g = Mod(eval(Str(pe))^2, p) \\ dirty hack to convert t_FFELEM to t_INT |
30 | pe=ffgen(minpoly(ffprimroot(ffgen(ffinit(q,1))),var),var) \\ get primitive element | 32 | g = Mod(4, p) |
31 | 1/(fforder(pe) == q-1) \\ error out, if ord(pe) is wrong | ||
32 | g = Mod(eval(Str(pe)), q) \\ dirty hack to convert t_FFELEM to t_INT | ||
33 | 33 | ||
34 | \\\\\\\\\\\\ | 34 | \\\\\\\\\\\\ |
35 | \\ PROLOG | 35 | \\ PROLOG |
@@ -66,8 +66,8 @@ beta = matrix(n,k,i,j, g^r[i,j]) | |||
66 | \\ multiplicative shares | 66 | \\ multiplicative shares |
67 | \\ first index level = owning agent id (multiplicative share) | 67 | \\ first index level = owning agent id (multiplicative share) |
68 | \\ second index level = agent id, price id | 68 | \\ second index level = agent id, price id |
69 | Gamma = vector(n,a,matrix(n,k,i,j, ( prod(h=1,n,prod(d=j+1,k,alpha[h,d])) * prod(d=1,j-1,alpha[i,d]) * prod(h=1,i-1,alpha[h,j]) )^m[a][i,j] )) | 69 | Gamma = vector(n,a,matrix(n,k,i,j, ( prod(h=1,n,prod(d=j,k,alpha[h,d])*prod(d=j+1,k,alpha[h,d])) * prod(d=1,j,alpha[i,d])^(2*M+2) / g^(2*M+1) )^(m[a][i,j]) )) |
70 | Delta = vector(n,a,matrix(n,k,i,j, ( prod(h=1,n,prod(d=j+1,k, beta[h,d])) * prod(d=1,j-1, beta[i,d]) * prod(h=1,i-1, beta[h,j]) )^m[a][i,j] )) | 70 | Delta = vector(n,a,matrix(n,k,i,j, ( prod(h=1,n,prod(d=j,k, beta[h,d])*prod(d=j+1,k, beta[h,d])) * prod(d=1,j, beta[i,d])^(2*M+2) )^(m[a][i,j]) )) |
71 | 71 | ||
72 | \\\\\\\\\\\\ | 72 | \\\\\\\\\\\\ |
73 | \\ ROUND3 | 73 | \\ ROUND3 |
@@ -84,8 +84,7 @@ Phi = vector(n,a,matrix(n,k,i,j, prod(h=1,n,Delta[h][i,j])^x[a] )) | |||
84 | 84 | ||
85 | \\ winner matrix | 85 | \\ winner matrix |
86 | v = matrix(n,k,a,j, prod(i=1,n,Gamma[i][a,j]) / prod(i=1,n,Phi[i][a,j]) ) | 86 | v = matrix(n,k,a,j, prod(i=1,n,Gamma[i][a,j]) / prod(i=1,n,Phi[i][a,j]) ) |
87 | |||
88 | vi = lift(v) | 87 | vi = lift(v) |
89 | 88 | ||
90 | print("bids are: ", bid) | 89 | print("bids are: ", bid) |
91 | for(X=1,n, if(vecmin(vi[X,])==1, print("And the winner is ", X) )) | 90 | for(X=1,n, if(vecmin(vi[X,],&i)==1, print("And the winner is ", X, " with the price ", i) )) |