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
|
/* General utility functions */
#ifndef UTILS_H
#define UTILS_H
#include <stdbool.h>
#define min(a,b) (((a) < (b)) ? (a) : (b))
#define max(a,b) (((a) > (b)) ? (a) : (b))
/* Some useful constants */
#define pow2to11 2048
#define pow2to12 4096
#define pow3to7 2187
#define pow3to8 6561
#define pow12to4 20736
#define factorial4 24
#define factorial6 720
#define factorial8 40320
#define factorial12 479001600
#define binom12on4 495
#define binom8on4 70
/* Generic utility functions */
void swap(int *a, int *b);
void intarrcopy(int *src, int *dst, int n);
int sum(int *a, int n);
bool is_perm(int *a, int n);
bool is_perm(int *a, int n);
/* Standard mathematical functions */
int powint(int a, int b);
int factorial(int n);
int binomial(int n, int k);
/* Converts the integer a to its representation in base b (first n digits
* only) and saves the result in r. */
void int_to_digit_array(int a, int b, int n, int *r);
int digit_array_to_int(int *a, int n, int b);
/* Converts the first n-1 digits of a number to an array a of digits in base b;
* then adds one element to the array, so that the sum of the elements of a is
* zero modulo b.
* This is used for determing the edge orientation from an 11-bits integer or
* the corner orientation from a 7-trits integer. */
void int_to_sum_zero_array(int x, int b, int n, int *a);
/* Converts a permutation on [0..(n-1)] into the integer i which is the index
* of the permutation in the sorted list of all n! such permutations. */
int perm_to_index(int *a, int n);
void index_to_perm(int p, int n, int *r);
/* Determine the sign of a permutation */
int perm_sign(int a[], int n);
/* Converts a k-element subset of a set from an array of n elements, of which k
* are 1 and n-k are 0, to its index in the sorted list of all such subsets. */
int subset_to_index(int *a, int n, int k);
void index_to_subset(int s, int n, int k, int *r);
int ordered_subset_to_index(int *a, int n, int k);
void index_to_ordered_subset(int s, int n, int k, int *r);
void apply_permutation(int *perm, int *set, int n);
/* b[i] = (a[i]+b[i])%m for i=1,...,n */
void sum_arrays_mod(int *a, int *b, int n, int m);
#endif
|
