int sieve(int *x, int *y, int b, int n, int nx) {
  int i;
  int ny = 1;
  for (i = 1; i < nx; i++) {
    if (i != b && i % n == b % n)
      continue;
    y[ny++] = x[i];
  }
  for (i = 1; i < ny; i++) {
    x[i] = y[i];
  }
  return ny;
}


x86 assembly (clang, -O0):
<sieve>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
movl   $0x1,-0x24(%rbp)
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    118d <sieve+0x8d>
mov    -0x20(%rbp),%eax
cmp    -0x14(%rbp),%eax
je     115f <sieve+0x5f>
mov    -0x20(%rbp),%eax
cltd
idivl  -0x18(%rbp)
mov    %edx,-0x28(%rbp)
mov    -0x14(%rbp),%eax
cltd
idivl  -0x18(%rbp)
mov    %edx,%eax
mov    -0x28(%rbp),%edx
cmp    %eax,%edx
jne    115f <sieve+0x5f>
jmp    117f <sieve+0x7f>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    %ecx,%esi
add    $0x1,%esi
mov    %esi,-0x24(%rbp)
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1124 <sieve+0x24>
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    11c4 <sieve+0xc4>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1194 <sieve+0x94>
mov    -0x24(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<sieve>:
push   %rbx
mov    $0x1,%r11d
cmp    $0x2,%r8d
jl     12c7 <sieve+0x1c7>
mov    %edx,%r9d
mov    %r8d,%ebx
mov    $0x1,%r11d
mov    $0x1,%r10d
jmp    1147 <sieve+0x47>
cs nopw 0x0(%rax,%rax,1)
nop
mov    (%rdi,%r10,4),%eax
movslq %r11d,%rdx
add    $0x1,%r11d
mov    %eax,(%rsi,%rdx,4)
add    $0x1,%r10
cmp    %r10,%rbx
je     1162 <sieve+0x62>
cmp    %r10,%r9
je     1130 <sieve+0x30>
mov    %r10d,%eax
cltd
idiv   %ecx
mov    %edx,%r8d
mov    %r9d,%eax
cltd
idiv   %ecx
cmp    %edx,%r8d
jne    1130 <sieve+0x30>
jmp    113e <sieve+0x3e>
cmp    $0x2,%r11d
jl     12c7 <sieve+0x1c7>
mov    %r11d,%r9d
lea    -0x1(%r9),%r8
mov    $0x1,%ecx
cmp    $0x8,%r8
jb     1261 <sieve+0x161>
lea    0x4(%rdi),%rax
lea    (%rsi,%r9,4),%rdx
cmp    %rdx,%rax
jae    11a0 <sieve+0xa0>
lea    (%rdi,%r9,4),%rax
lea    0x4(%rsi),%rdx
cmp    %rax,%rdx
jb     1261 <sieve+0x161>
mov    %r8,%rcx
and    $0xfffffffffffffff8,%rcx
lea    -0x8(%rcx),%rax
mov    %rax,%rdx
shr    $0x3,%rdx
add    $0x1,%rdx
mov    %edx,%ebx
and    $0x3,%ebx
cmp    $0x18,%rax
jae    11c5 <sieve+0xc5>
xor    %eax,%eax
jmp    122a <sieve+0x12a>
and    $0xfffffffffffffffc,%rdx
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
movups 0x4(%rsi,%rax,4),%xmm0
movups 0x14(%rsi,%rax,4),%xmm1
movups %xmm0,0x4(%rdi,%rax,4)
movups %xmm1,0x14(%rdi,%rax,4)
movups 0x24(%rsi,%rax,4),%xmm0
movups 0x34(%rsi,%rax,4),%xmm1
movups %xmm0,0x24(%rdi,%rax,4)
movups %xmm1,0x34(%rdi,%rax,4)
movups 0x44(%rsi,%rax,4),%xmm0
movups 0x54(%rsi,%rax,4),%xmm1
movups %xmm0,0x44(%rdi,%rax,4)
movups %xmm1,0x54(%rdi,%rax,4)
movups 0x64(%rsi,%rax,4),%xmm0
movups 0x74(%rsi,%rax,4),%xmm1
movups %xmm0,0x64(%rdi,%rax,4)
movups %xmm1,0x74(%rdi,%rax,4)
add    $0x20,%rax
add    $0xfffffffffffffffc,%rdx
jne    11d0 <sieve+0xd0>
test   %rbx,%rbx
je     1258 <sieve+0x158>
nop
lea    0x0(,%rax,4),%rdx
or     $0x4,%rdx
movups (%rsi,%rdx,1),%xmm0
movups 0x10(%rsi,%rdx,1),%xmm1
movups %xmm0,(%rdi,%rdx,1)
movups %xmm1,0x10(%rdi,%rdx,1)
add    $0x8,%rax
add    $0xffffffffffffffff,%rbx
jne    1230 <sieve+0x130>
cmp    %rcx,%r8
je     12c7 <sieve+0x1c7>
or     $0x1,%rcx
mov    %r9d,%edx
sub    %ecx,%edx
mov    %rcx,%rax
not    %rax
add    %r9,%rax
and    $0x3,%rdx
je     1290 <sieve+0x190>
cs nopw 0x0(%rax,%rax,1)
nop
mov    (%rsi,%rcx,4),%ebx
mov    %ebx,(%rdi,%rcx,4)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rdx
jne    1280 <sieve+0x180>
cmp    $0x3,%rax
jb     12c7 <sieve+0x1c7>
cs nopw 0x0(%rax,%rax,1)
mov    (%rsi,%rcx,4),%eax
mov    %eax,(%rdi,%rcx,4)
mov    0x4(%rsi,%rcx,4),%eax
mov    %eax,0x4(%rdi,%rcx,4)
mov    0x8(%rsi,%rcx,4),%eax
mov    %eax,0x8(%rdi,%rcx,4)
mov    0xc(%rsi,%rcx,4),%eax
mov    %eax,0xc(%rdi,%rcx,4)
add    $0x4,%rcx
cmp    %rcx,%r9
jne    12a0 <sieve+0x1a0>
mov    %r11d,%eax
pop    %rbx
ret

LLVM IR (-O0):
define i32 @sieve(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
store i32* %0, i32** %6
store i32* %1, i32** %7
store i32 %2, i32* %8
store i32 %3, i32* %9
store i32 %4, i32* %10
store i32 1, i32* %12
store i32 1, i32* %11
br label %13
13: ; preds = %41, %5
%14 = load i32, i32* %11
%15 = load i32, i32* %10
%16 = icmp slt i32 %14, %15
br i1 %16, label %17, label %44
17: ; preds = %13
%18 = load i32, i32* %11
%19 = load i32, i32* %8
%20 = icmp ne i32 %18, %19
br i1 %20, label %21, label %30
21: ; preds = %17
%22 = load i32, i32* %11
%23 = load i32, i32* %9
%24 = srem i32 %22, %23
%25 = load i32, i32* %8
%26 = load i32, i32* %9
%27 = srem i32 %25, %26
%28 = icmp eq i32 %24, %27
br i1 %28, label %29, label %30
29: ; preds = %21
br label %41
30: ; preds = %21, %17
%31 = load i32*, i32** %6
%32 = load i32, i32* %11
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds i32, i32* %31, i64 %33
%35 = load i32, i32* %34
%36 = load i32*, i32** %7
%37 = load i32, i32* %12
%38 = add nsw i32 %37, 1
store i32 %38, i32* %12
%39 = sext i32 %37 to i64
%40 = getelementptr inbounds i32, i32* %36, i64 %39
store i32 %35, i32* %40
br label %41
41: ; preds = %30, %29
%42 = load i32, i32* %11
%43 = add nsw i32 %42, 1
store i32 %43, i32* %11
br label %13
44: ; preds = %13
store i32 1, i32* %11
br label %45
45: ; preds = %59, %44
%46 = load i32, i32* %11
%47 = load i32, i32* %12
%48 = icmp slt i32 %46, %47
br i1 %48, label %49, label %62
49: ; preds = %45
%50 = load i32*, i32** %7
%51 = load i32, i32* %11
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds i32, i32* %50, i64 %52
%54 = load i32, i32* %53
%55 = load i32*, i32** %6
%56 = load i32, i32* %11
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds i32, i32* %55, i64 %57
store i32 %54, i32* %58
br label %59
59: ; preds = %49
%60 = load i32, i32* %11
%61 = add nsw i32 %60, 1
store i32 %61, i32* %11
br label %45
62: ; preds = %45
%63 = load i32, i32* %12
ret i32 %63
}

LLVM IR (-O2):
define i32 @sieve(i32* nocapture %0, i32* nocapture %1, i32 %2, i32 %3, i32 %4) {
%6 = icmp sgt i32 %4, 1
br i1 %6, label %7, label %163
7: ; preds = %5
%8 = zext i32 %2 to i64
%9 = zext i32 %4 to i64
br label %125
10: ; preds = %140
%11 = icmp sgt i32 %141, 1
br i1 %11, label %12, label %163
12: ; preds = %10
%13 = zext i32 %141 to i64
%14 = add nsw i64 %13, -1
%15 = icmp ult i64 %14, 8
br i1 %15, label %106, label %16
16: ; preds = %12
%17 = getelementptr i32, i32* %0, i64 1
%18 = getelementptr i32, i32* %0, i64 %13
%19 = getelementptr i32, i32* %1, i64 1
%20 = getelementptr i32, i32* %1, i64 %13
%21 = icmp ult i32* %17, %20
%22 = icmp ult i32* %19, %18
%23 = and i1 %21, %22
br i1 %23, label %106, label %24
24: ; preds = %16
%25 = and i64 %14, -8
%26 = or i64 %25, 1
%27 = add nsw i64 %25, -8
%28 = lshr exact i64 %27, 3
%29 = add nuw nsw i64 %28, 1
%30 = and i64 %29, 3
%31 = icmp ult i64 %27, 24
br i1 %31, label %84, label %32
32: ; preds = %24
%33 = and i64 %29, 4611686018427387900
br label %34
34: ; preds = %34, %32
%35 = phi i64 [ 0, %32 ], [ %81, %34 ]
%36 = phi i64 [ 0, %32 ], [ %82, %34 ]
%37 = or i64 %35, 1
%38 = getelementptr inbounds i32, i32* %1, i64 %37
%39 = bitcast i32* %38 to <4 x i32>*
%40 = load <4 x i32>, <4 x i32>* %39
%41 = getelementptr inbounds i32, i32* %38, i64 4
%42 = bitcast i32* %41 to <4 x i32>*
%43 = load <4 x i32>, <4 x i32>* %42
%44 = getelementptr inbounds i32, i32* %0, i64 %37
%45 = bitcast i32* %44 to <4 x i32>*
store <4 x i32> %40, <4 x i32>* %45
%46 = getelementptr inbounds i32, i32* %44, i64 4
%47 = bitcast i32* %46 to <4 x i32>*
store <4 x i32> %43, <4 x i32>* %47
%48 = or i64 %35, 9
%49 = getelementptr inbounds i32, i32* %1, i64 %48
%50 = bitcast i32* %49 to <4 x i32>*
%51 = load <4 x i32>, <4 x i32>* %50
%52 = getelementptr inbounds i32, i32* %49, i64 4
%53 = bitcast i32* %52 to <4 x i32>*
%54 = load <4 x i32>, <4 x i32>* %53
%55 = getelementptr inbounds i32, i32* %0, i64 %48
%56 = bitcast i32* %55 to <4 x i32>*
store <4 x i32> %51, <4 x i32>* %56
%57 = getelementptr inbounds i32, i32* %55, i64 4
%58 = bitcast i32* %57 to <4 x i32>*
store <4 x i32> %54, <4 x i32>* %58
%59 = or i64 %35, 17
%60 = getelementptr inbounds i32, i32* %1, i64 %59
%61 = bitcast i32* %60 to <4 x i32>*
%62 = load <4 x i32>, <4 x i32>* %61
%63 = getelementptr inbounds i32, i32* %60, i64 4
%64 = bitcast i32* %63 to <4 x i32>*
%65 = load <4 x i32>, <4 x i32>* %64
%66 = getelementptr inbounds i32, i32* %0, i64 %59
%67 = bitcast i32* %66 to <4 x i32>*
store <4 x i32> %62, <4 x i32>* %67
%68 = getelementptr inbounds i32, i32* %66, i64 4
%69 = bitcast i32* %68 to <4 x i32>*
store <4 x i32> %65, <4 x i32>* %69
%70 = or i64 %35, 25
%71 = getelementptr inbounds i32, i32* %1, i64 %70
%72 = bitcast i32* %71 to <4 x i32>*
%73 = load <4 x i32>, <4 x i32>* %72
%74 = getelementptr inbounds i32, i32* %71, i64 4
%75 = bitcast i32* %74 to <4 x i32>*
%76 = load <4 x i32>, <4 x i32>* %75
%77 = getelementptr inbounds i32, i32* %0, i64 %70
%78 = bitcast i32* %77 to <4 x i32>*
store <4 x i32> %73, <4 x i32>* %78
%79 = getelementptr inbounds i32, i32* %77, i64 4
%80 = bitcast i32* %79 to <4 x i32>*
store <4 x i32> %76, <4 x i32>* %80
%81 = add nuw i64 %35, 32
%82 = add i64 %36, 4
%83 = icmp eq i64 %82, %33
br i1 %83, label %84, label %34
84: ; preds = %34, %24
%85 = phi i64 [ 0, %24 ], [ %81, %34 ]
%86 = icmp eq i64 %30, 0
br i1 %86, label %104, label %87
87: ; preds = %84, %87
%88 = phi i64 [ %101, %87 ], [ %85, %84 ]
%89 = phi i64 [ %102, %87 ], [ 0, %84 ]
%90 = or i64 %88, 1
%91 = getelementptr inbounds i32, i32* %1, i64 %90
%92 = bitcast i32* %91 to <4 x i32>*
%93 = load <4 x i32>, <4 x i32>* %92
%94 = getelementptr inbounds i32, i32* %91, i64 4
%95 = bitcast i32* %94 to <4 x i32>*
%96 = load <4 x i32>, <4 x i32>* %95
%97 = getelementptr inbounds i32, i32* %0, i64 %90
%98 = bitcast i32* %97 to <4 x i32>*
store <4 x i32> %93, <4 x i32>* %98
%99 = getelementptr inbounds i32, i32* %97, i64 4
%100 = bitcast i32* %99 to <4 x i32>*
store <4 x i32> %96, <4 x i32>* %100
%101 = add nuw i64 %88, 8
%102 = add i64 %89, 1
%103 = icmp eq i64 %102, %30
br i1 %103, label %104, label %87
104: ; preds = %87, %84
%105 = icmp eq i64 %14, %25
br i1 %105, label %163, label %106
106: ; preds = %16, %12, %104
%107 = phi i64 [ 1, %16 ], [ 1, %12 ], [ %26, %104 ]
%108 = sub nsw i64 %13, %107
%109 = xor i64 %107, -1
%110 = add nsw i64 %109, %13
%111 = and i64 %108, 3
%112 = icmp eq i64 %111, 0
br i1 %112, label %122, label %113
113: ; preds = %106, %113
%114 = phi i64 [ %119, %113 ], [ %107, %106 ]
%115 = phi i64 [ %120, %113 ], [ 0, %106 ]
%116 = getelementptr inbounds i32, i32* %1, i64 %114
%117 = load i32, i32* %116
%118 = getelementptr inbounds i32, i32* %0, i64 %114
store i32 %117, i32* %118
%119 = add nuw nsw i64 %114, 1
%120 = add i64 %115, 1
%121 = icmp eq i64 %120, %111
br i1 %121, label %122, label %113
122: ; preds = %113, %106
%123 = phi i64 [ %107, %106 ], [ %119, %113 ]
%124 = icmp ult i64 %110, 3
br i1 %124, label %163, label %144
125: ; preds = %7, %140
%126 = phi i64 [ 1, %7 ], [ %142, %140 ]
%127 = phi i32 [ 1, %7 ], [ %141, %140 ]
%128 = icmp eq i64 %126, %8
br i1 %128, label %134, label %129
129: ; preds = %125
%130 = trunc i64 %126 to i32
%131 = srem i32 %130, %3
%132 = srem i32 %2, %3
%133 = icmp eq i32 %131, %132
br i1 %133, label %140, label %134
134: ; preds = %129, %125
%135 = getelementptr inbounds i32, i32* %0, i64 %126
%136 = load i32, i32* %135
%137 = add nsw i32 %127, 1
%138 = sext i32 %127 to i64
%139 = getelementptr inbounds i32, i32* %1, i64 %138
store i32 %136, i32* %139
br label %140
140: ; preds = %129, %134
%141 = phi i32 [ %127, %129 ], [ %137, %134 ]
%142 = add nuw nsw i64 %126, 1
%143 = icmp eq i64 %142, %9
br i1 %143, label %10, label %125
144: ; preds = %122, %144
%145 = phi i64 [ %161, %144 ], [ %123, %122 ]
%146 = getelementptr inbounds i32, i32* %1, i64 %145
%147 = load i32, i32* %146
%148 = getelementptr inbounds i32, i32* %0, i64 %145
store i32 %147, i32* %148
%149 = add nuw nsw i64 %145, 1
%150 = getelementptr inbounds i32, i32* %1, i64 %149
%151 = load i32, i32* %150
%152 = getelementptr inbounds i32, i32* %0, i64 %149
store i32 %151, i32* %152
%153 = add nuw nsw i64 %145, 2
%154 = getelementptr inbounds i32, i32* %1, i64 %153
%155 = load i32, i32* %154
%156 = getelementptr inbounds i32, i32* %0, i64 %153
store i32 %155, i32* %156
%157 = add nuw nsw i64 %145, 3
%158 = getelementptr inbounds i32, i32* %1, i64 %157
%159 = load i32, i32* %158
%160 = getelementptr inbounds i32, i32* %0, i64 %157
store i32 %159, i32* %160
%161 = add nuw nsw i64 %145, 4
%162 = icmp eq i64 %161, %13
br i1 %162, label %163, label %144
163: ; preds = %122, %144, %104, %5, %10
%164 = phi i32 [ %141, %10 ], [ 1, %5 ], [ %141, %104 ], [ %141, %144 ], [ %141, %122 ]
ret i32 %164
}