#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ printf(char *, int);
int triplet2index(int, int, int, int);
int motif_triplet_distribution(int **matrix, int n_rows, int max_n_motifs,
                               int *motifs_triplets_list, int n_motifs) {
  int index;
  int i, j, k, ll;
  int i1, i2, i3;
  int n = 0;
  int max_index;
  if (max_n_motifs < 3) {
    printf("sorry. i cannot compute a distribution of triplets with "
           "max_n_motifs=%d\n",
           max_n_motifs);
    return n;
  }
  max_index = n_motifs * n_motifs * n_motifs;
  for (i = 1; i <= max_index; i++)
    motifs_triplets_list[i] = 0;
  for (i = 1; i <= n_rows; i++) {
    for (j = 1; j < max_n_motifs; j++) {
      i1 = matrix[i][j];
      if (i1 > 0) {
        for (k = (j + 1); k <= max_n_motifs; k++) {
          i2 = matrix[i][k];
          if (i2 > 0) {
            for (ll = (k + 1); ll <= max_n_motifs; ll++) {
              i3 = matrix[i][ll];
              if (i3 > 0) {
                index = triplet2index(i1, i2, i3, n_motifs);
                motifs_triplets_list[index] = motifs_triplets_list[index] + 1;
                n++;
              }
            }
          }
        }
      }
    }
  }
  return n;
}


x86 assembly (clang, -O0):
<motif_triplet_distribution>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %edx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
movl   $0x0,-0x48(%rbp)
cmpl   $0x3,-0x18(%rbp)
jge    1165 <motif_triplet_distribution+0x45>
mov    -0x18(%rbp),%esi
lea    0xeab(%rip),%rdi
call   1030 <printf@plt>
mov    -0x48(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    12e5 <motif_triplet_distribution+0x1c5>
mov    -0x24(%rbp),%eax
imul   -0x24(%rbp),%eax
imul   -0x24(%rbp),%eax
mov    %eax,-0x4c(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4c(%rbp),%eax
jg     11a3 <motif_triplet_distribution+0x83>
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movl   $0x0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    117a <motif_triplet_distribution+0x5a>
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     12df <motif_triplet_distribution+0x1bf>
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    12cc <motif_triplet_distribution+0x1ac>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x3c(%rbp)
cmpl   $0x0,-0x3c(%rbp)
jle    12b9 <motif_triplet_distribution+0x199>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     12b4 <motif_triplet_distribution+0x194>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x40(%rbp)
cmpl   $0x0,-0x40(%rbp)
jle    12a1 <motif_triplet_distribution+0x181>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     129c <motif_triplet_distribution+0x17c>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x44(%rbp)
cmpl   $0x0,-0x44(%rbp)
jle    1289 <motif_triplet_distribution+0x169>
mov    -0x3c(%rbp),%edi
mov    -0x40(%rbp),%esi
mov    -0x44(%rbp),%edx
mov    -0x24(%rbp),%ecx
call   1040 <triplet2index@plt>
mov    %eax,-0x28(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
add    $0x1,%edx
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    128e <motif_triplet_distribution+0x16e>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1227 <motif_triplet_distribution+0x107>
jmp    12a1 <motif_triplet_distribution+0x181>
jmp    12a6 <motif_triplet_distribution+0x186>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11f2 <motif_triplet_distribution+0xd2>
jmp    12b9 <motif_triplet_distribution+0x199>
jmp    12be <motif_triplet_distribution+0x19e>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11bd <motif_triplet_distribution+0x9d>
jmp    12d1 <motif_triplet_distribution+0x1b1>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11aa <motif_triplet_distribution+0x8a>
mov    -0x48(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<motif_triplet_distribution>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %rcx,0x20(%rsp)
mov    %edx,%ebp
mov    %rdi,0x8(%rsp)
cmp    $0x2,%edx
jg     1162 <motif_triplet_distribution+0x32>
lea    0xeaa(%rip),%rdi
mov    %ebp,%esi
call   1030 <printf@plt>
jmp    1299 <motif_triplet_distribution+0x169>
mov    %r8d,%r14d
mov    %esi,%ebx
mov    %r8d,%eax
imul   %r14d,%eax
imul   %r8d,%eax
test   %eax,%eax
jle    118c <motif_triplet_distribution+0x5c>
mov    0x20(%rsp),%rcx
lea    0x4(%rcx),%rdi
mov    %eax,%edx
shl    $0x2,%rdx
xor    %esi,%esi
call   1040 <memset@plt>
test   %ebx,%ebx
jle    1299 <motif_triplet_distribution+0x169>
mov    %r14d,0x14(%rsp)
add    $0x1,%ebx
mov    %rbx,0x28(%rsp)
mov    %ebp,%edx
xor    %ebx,%ebx
mov    $0x1,%r12d
mov    %rdx,0x38(%rsp)
jmp    11cf <motif_triplet_distribution+0x9f>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%r12
cmp    0x28(%rsp),%r12
je     129b <motif_triplet_distribution+0x16b>
mov    $0x3,%eax
mov    %rax,0x18(%rsp)
mov    $0x1,%ecx
jmp    11f4 <motif_triplet_distribution+0xc4>
mov    0x30(%rsp),%rcx
add    $0x1,%rcx
addq   $0x1,0x18(%rsp)
cmp    %rdx,%rcx
je     11c0 <motif_triplet_distribution+0x90>
mov    0x8(%rsp),%rax
mov    (%rax,%r12,8),%rax
mov    (%rax,%rcx,4),%r15d
test   %r15d,%r15d
jle    11e5 <motif_triplet_distribution+0xb5>
mov    %rcx,0x30(%rsp)
mov    0x18(%rsp),%r13
lea    0x1(%rcx),%rsi
mov    %rsi,(%rsp)
mov    0x4(%rax,%rcx,4),%r14d
test   %r14d,%r14d
jle    127a <motif_triplet_distribution+0x14a>
add    $0x2,%ecx
cmp    %ebp,%ecx
jg     127a <motif_triplet_distribution+0x14a>
mov    %r13,0x40(%rsp)
xchg   %ax,%ax
mov    (%rax,%r13,4),%edx
test   %edx,%edx
jle    1255 <motif_triplet_distribution+0x125>
mov    %r15d,%edi
mov    %r14d,%esi
mov    0x14(%rsp),%ecx
call   1050 <triplet2index@plt>
cltq
mov    0x20(%rsp),%rcx
addl   $0x1,(%rcx,%rax,4)
add    $0x1,%ebx
cmp    %ebp,%r13d
jge    1270 <motif_triplet_distribution+0x140>
add    $0x1,%r13
mov    0x8(%rsp),%rax
mov    (%rax,%r12,8),%rax
jmp    1230 <motif_triplet_distribution+0x100>
nopl   0x0(%rax)
mov    0x38(%rsp),%rdx
mov    0x40(%rsp),%r13
mov    (%rsp),%rcx
cmp    %rdx,%rcx
je     11e0 <motif_triplet_distribution+0xb0>
add    $0x1,%r13
mov    0x8(%rsp),%rax
mov    (%rax,%r12,8),%rax
jmp    1210 <motif_triplet_distribution+0xe0>
xor    %ebx,%ebx
mov    %ebx,%eax
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [73 x i8] c"sorry. i cannot compute a distribution of triplets with max_n_motifs=%d\0A\00"
define i32 @motif_triplet_distribution(i32** %0, i32 %1, i32 %2, i32* %3, i32 %4) {
%6 = alloca i32
%7 = alloca i32**
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32*
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
store i32** %0, i32*** %7
store i32 %1, i32* %8
store i32 %2, i32* %9
store i32* %3, i32** %10
store i32 %4, i32* %11
store i32 0, i32* %20
%22 = load i32, i32* %9
%23 = icmp slt i32 %22, 3
br i1 %23, label %24, label %28
24: ; preds = %5
%25 = load i32, i32* %9
%26 = call i32 @printf(i8* getelementptr inbounds ([73 x i8], [73 x i8]* @.str, i64 0, i64 0), i32 %25)
%27 = load i32, i32* %20
store i32 %27, i32* %6
br label %144
28: ; preds = %5
%29 = load i32, i32* %11
%30 = load i32, i32* %11
%31 = mul nsw i32 %29, %30
%32 = load i32, i32* %11
%33 = mul nsw i32 %31, %32
store i32 %33, i32* %21
store i32 1, i32* %13
br label %34
34: ; preds = %43, %28
%35 = load i32, i32* %13
%36 = load i32, i32* %21
%37 = icmp sle i32 %35, %36
br i1 %37, label %38, label %46
38: ; preds = %34
%39 = load i32*, i32** %10
%40 = load i32, i32* %13
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds i32, i32* %39, i64 %41
store i32 0, i32* %42
br label %43
43: ; preds = %38
%44 = load i32, i32* %13
%45 = add nsw i32 %44, 1
store i32 %45, i32* %13
br label %34
46: ; preds = %34
store i32 1, i32* %13
br label %47
47: ; preds = %139, %46
%48 = load i32, i32* %13
%49 = load i32, i32* %8
%50 = icmp sle i32 %48, %49
br i1 %50, label %51, label %142
51: ; preds = %47
store i32 1, i32* %14
br label %52
52: ; preds = %135, %51
%53 = load i32, i32* %14
%54 = load i32, i32* %9
%55 = icmp slt i32 %53, %54
br i1 %55, label %56, label %138
56: ; preds = %52
%57 = load i32**, i32*** %7
%58 = load i32, i32* %13
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds i32*, i32** %57, i64 %59
%61 = load i32*, i32** %60
%62 = load i32, i32* %14
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds i32, i32* %61, i64 %63
%65 = load i32, i32* %64
store i32 %65, i32* %17
%66 = load i32, i32* %17
%67 = icmp sgt i32 %66, 0
br i1 %67, label %68, label %134
68: ; preds = %56
%69 = load i32, i32* %14
%70 = add nsw i32 %69, 1
store i32 %70, i32* %15
br label %71
71: ; preds = %130, %68
%72 = load i32, i32* %15
%73 = load i32, i32* %9
%74 = icmp sle i32 %72, %73
br i1 %74, label %75, label %133
75: ; preds = %71
%76 = load i32**, i32*** %7
%77 = load i32, i32* %13
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds i32*, i32** %76, i64 %78
%80 = load i32*, i32** %79
%81 = load i32, i32* %15
%82 = sext i32 %81 to i64
%83 = getelementptr inbounds i32, i32* %80, i64 %82
%84 = load i32, i32* %83
store i32 %84, i32* %18
%85 = load i32, i32* %18
%86 = icmp sgt i32 %85, 0
br i1 %86, label %87, label %129
87: ; preds = %75
%88 = load i32, i32* %15
%89 = add nsw i32 %88, 1
store i32 %89, i32* %16
br label %90
90: ; preds = %125, %87
%91 = load i32, i32* %16
%92 = load i32, i32* %9
%93 = icmp sle i32 %91, %92
br i1 %93, label %94, label %128
94: ; preds = %90
%95 = load i32**, i32*** %7
%96 = load i32, i32* %13
%97 = sext i32 %96 to i64
%98 = getelementptr inbounds i32*, i32** %95, i64 %97
%99 = load i32*, i32** %98
%100 = load i32, i32* %16
%101 = sext i32 %100 to i64
%102 = getelementptr inbounds i32, i32* %99, i64 %101
%103 = load i32, i32* %102
store i32 %103, i32* %19
%104 = load i32, i32* %19
%105 = icmp sgt i32 %104, 0
br i1 %105, label %106, label %124
106: ; preds = %94
%107 = load i32, i32* %17
%108 = load i32, i32* %18
%109 = load i32, i32* %19
%110 = load i32, i32* %11
%111 = call i32 @triplet2index(i32 %107, i32 %108, i32 %109, i32 %110)
store i32 %111, i32* %12
%112 = load i32*, i32** %10
%113 = load i32, i32* %12
%114 = sext i32 %113 to i64
%115 = getelementptr inbounds i32, i32* %112, i64 %114
%116 = load i32, i32* %115
%117 = add nsw i32 %116, 1
%118 = load i32*, i32** %10
%119 = load i32, i32* %12
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds i32, i32* %118, i64 %120
store i32 %117, i32* %121
%122 = load i32, i32* %20
%123 = add nsw i32 %122, 1
store i32 %123, i32* %20
br label %124
124: ; preds = %106, %94
br label %125
125: ; preds = %124
%126 = load i32, i32* %16
%127 = add nsw i32 %126, 1
store i32 %127, i32* %16
br label %90
128: ; preds = %90
br label %129
129: ; preds = %128, %75
br label %130
130: ; preds = %129
%131 = load i32, i32* %15
%132 = add nsw i32 %131, 1
store i32 %132, i32* %15
br label %71
133: ; preds = %71
br label %134
134: ; preds = %133, %56
br label %135
135: ; preds = %134
%136 = load i32, i32* %14
%137 = add nsw i32 %136, 1
store i32 %137, i32* %14
br label %52
138: ; preds = %52
br label %139
139: ; preds = %138
%140 = load i32, i32* %13
%141 = add nsw i32 %140, 1
store i32 %141, i32* %13
br label %47
142: ; preds = %47
%143 = load i32, i32* %20
store i32 %143, i32* %6
br label %144
144: ; preds = %142, %24
%145 = load i32, i32* %6
ret i32 %145
}

LLVM IR (-O2):
@.str = [73 x i8] c"sorry. i cannot compute a distribution of triplets with max_n_motifs=%d\0A\00"
define i32 @motif_triplet_distribution(i32** nocapture %0, i32 %1, i32 %2, i32* nocapture %3, i32 %4) {
%6 = icmp slt i32 %2, 3
br i1 %6, label %7, label %9
7: ; preds = %5
%8 = tail call i32 @printf(i8* getelementptr inbounds ([73 x i8], [73 x i8]* @.str, i64 0, i64 0), i32 %2)
br label %84
9: ; preds = %5
%10 = mul nsw i32 %4, %4
%11 = mul i32 %10, %4
%12 = icmp slt i32 %11, 1
br i1 %12, label %18, label %13
13: ; preds = %9
%14 = getelementptr i32, i32* %3, i64 1
%15 = bitcast i32* %14 to i8*
%16 = zext i32 %11 to i64
%17 = shl nuw nsw i64 %16, 2
br label %18
18: ; preds = %13, %9
%19 = icmp slt i32 %1, 1
br i1 %19, label %84, label %20
20: ; preds = %18
%21 = add nuw i32 %1, 1
%22 = zext i32 %21 to i64
%23 = zext i32 %2 to i64
br label %24
24: ; preds = %20, %81
%25 = phi i64 [ 1, %20 ], [ %82, %81 ]
%26 = phi i32 [ 0, %20 ], [ %77, %81 ]
%27 = getelementptr inbounds i32*, i32** %0, i64 %25
br label %28
28: ; preds = %24, %76
%29 = phi i64 [ 1, %24 ], [ %78, %76 ]
%30 = phi i64 [ 3, %24 ], [ %79, %76 ]
%31 = phi i32 [ %26, %24 ], [ %77, %76 ]
%32 = load i32*, i32** %27
%33 = getelementptr inbounds i32, i32* %32, i64 %29
%34 = load i32, i32* %33
%35 = icmp sgt i32 %34, 0
br i1 %35, label %36, label %76
36: ; preds = %28, %73
%37 = phi i32* [ %75, %73 ], [ %32, %28 ]
%38 = phi i64 [ %41, %73 ], [ %29, %28 ]
%39 = phi i64 [ %74, %73 ], [ %30, %28 ]
%40 = phi i32 [ %71, %73 ], [ %31, %28 ]
%41 = add nuw nsw i64 %38, 1
%42 = getelementptr inbounds i32, i32* %37, i64 %41
%43 = load i32, i32* %42
%44 = icmp slt i32 %43, 1
%45 = trunc i64 %38 to i32
%46 = add i32 %45, 2
%47 = icmp sgt i32 %46, %2
%48 = select i1 %44, i1 true, i1 %47
br i1 %48, label %70, label %49
49: ; preds = %36, %67
%50 = phi i32* [ %69, %67 ], [ %37, %36 ]
%51 = phi i64 [ %68, %67 ], [ %39, %36 ]
%52 = phi i32 [ %64, %67 ], [ %40, %36 ]
%53 = getelementptr inbounds i32, i32* %50, i64 %51
%54 = load i32, i32* %53
%55 = icmp sgt i32 %54, 0
br i1 %55, label %56, label %63
56: ; preds = %49
%57 = tail call i32 @triplet2index(i32 %34, i32 %43, i32 %54, i32 %4)
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds i32, i32* %3, i64 %58
%60 = load i32, i32* %59
%61 = add nsw i32 %60, 1
store i32 %61, i32* %59
%62 = add nsw i32 %52, 1
br label %63
63: ; preds = %49, %56
%64 = phi i32 [ %62, %56 ], [ %52, %49 ]
%65 = trunc i64 %51 to i32
%66 = icmp slt i32 %65, %2
br i1 %66, label %67, label %70
67: ; preds = %63
%68 = add nuw nsw i64 %51, 1
%69 = load i32*, i32** %27
br label %49
70: ; preds = %63, %36
%71 = phi i32 [ %40, %36 ], [ %64, %63 ]
%72 = icmp eq i64 %41, %23
br i1 %72, label %76, label %73
73: ; preds = %70
%74 = add nuw nsw i64 %39, 1
%75 = load i32*, i32** %27
br label %36
76: ; preds = %70, %28
%77 = phi i32 [ %31, %28 ], [ %71, %70 ]
%78 = add nuw nsw i64 %29, 1
%79 = add nuw nsw i64 %30, 1
%80 = icmp eq i64 %78, %23
br i1 %80, label %81, label %28
81: ; preds = %76
%82 = add nuw nsw i64 %25, 1
%83 = icmp eq i64 %82, %22
br i1 %83, label %84, label %24
84: ; preds = %81, %18, %7
%85 = phi i32 [ 0, %7 ], [ 0, %18 ], [ %77, %81 ]
ret i32 %85
}