#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 *, long long);
long long sqrtl(long double);
int main() {
  long long b, bmax = sqrtl(2 * 10000000000ll - 1), s = 0;
#pragma omp parallel for reduction(+ : s) schedule(dynamic)
  for (b = 2; b <= bmax; b += 2) {
    long long c, cmax, t, diff;
    cmax = sqrtl((((long double)10000000000ll) * 4 * 10000000000ll + 1) /
                     (1 + b * b) -
                 1);
    t = 1 + b * b;
    diff = -1;
    for (c = b; c <= cmax; diff += (4 * c + 4) * (1 + b * b), c += 2) {
      while (diff > 0) {
        diff -= 2 * t + 1, t++;
      }
      if (diff == 0 && t % 2 == 0)
        s += t / 2;
      diff -= 4 * b * t + 4 * b * b, t += 2 * b;
    }
  }
  printf("%lld\n", s);
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
movl   $0x0,-0x4(%rbp)
mov    %rsp,%rax
fldl   0xec8(%rip)
fstpt  (%rax)
call   1040 <sqrtl@plt>
mov    %rax,-0x18(%rbp)
movq   $0x0,-0x20(%rbp)
movq   $0x2,-0x10(%rbp)
mov    -0x10(%rbp),%rax
cmp    -0x18(%rbp),%rax
jg     12b4 <main+0x194>
mov    -0x10(%rbp),%rax
imul   %rax,%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
fildll -0x48(%rbp)
fldl   0xe8e(%rip)
fdivp  %st,%st(1)
fld1
fchs
faddp  %st,%st(1)
mov    %rsp,%rax
fstpt  (%rax)
call   1040 <sqrtl@plt>
mov    %rax,-0x30(%rbp)
mov    -0x10(%rbp),%rax
imul   -0x10(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
movq   $0xffffffffffffffff,-0x40(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x30(%rbp),%rax
jg     129e <main+0x17e>
jmp    11c4 <main+0xa4>
cmpq   $0x0,-0x40(%rbp)
jle    11f7 <main+0xd7>
mov    -0x38(%rbp),%rcx
shl    $0x1,%rcx
add    $0x1,%rcx
mov    -0x40(%rbp),%rax
sub    %rcx,%rax
mov    %rax,-0x40(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    11c4 <main+0xa4>
cmpq   $0x0,-0x40(%rbp)
jne    1230 <main+0x110>
mov    -0x38(%rbp),%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
cmp    $0x0,%rdx
jne    1230 <main+0x110>
mov    -0x38(%rbp),%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
add    -0x20(%rbp),%rax
mov    %rax,-0x20(%rbp)
mov    -0x10(%rbp),%rcx
shl    $0x2,%rcx
imul   -0x38(%rbp),%rcx
mov    -0x10(%rbp),%rax
shl    $0x2,%rax
imul   -0x10(%rbp),%rax
add    %rax,%rcx
mov    -0x40(%rbp),%rax
sub    %rcx,%rax
mov    %rax,-0x40(%rbp)
mov    -0x10(%rbp),%rax
shl    $0x1,%rax
add    -0x38(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x28(%rbp),%rax
shl    $0x2,%rax
add    $0x4,%rax
mov    -0x10(%rbp),%rcx
imul   -0x10(%rbp),%rcx
add    $0x1,%rcx
imul   %rcx,%rax
add    -0x40(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x28(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x28(%rbp)
jmp    11b1 <main+0x91>
jmp    12a3 <main+0x183>
mov    -0x10(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x10(%rbp)
jmp    1153 <main+0x33>
mov    -0x20(%rbp),%rsi
lea    0xd51(%rip),%rdi
call   1030 <printf@plt>
mov    -0x4(%rbp),%eax
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x30,%rsp
fldl   0xecd(%rip)
fstpt  (%rsp)
call   1040 <sqrtl@plt>
cmp    $0x2,%rax
jge    1163 <main+0x43>
xor    %r15d,%r15d
lea    0xec5(%rip),%rdi
mov    %r15,%rsi
call   1030 <printf@plt>
xor    %eax,%eax
add    $0x30,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret
mov    %rax,%r14
mov    $0x2,%r13d
xor    %r15d,%r15d
fld1
fchs
fstpt  0x24(%rsp)
jmp    1189 <main+0x69>
nopl   0x0(%rax)
add    $0x2,%r13
cmp    %r14,%r13
jg     1144 <main+0x24>
mov    %r13,%r12
imul   %r13,%r12
or     $0x1,%r12
mov    %r12,0x18(%rsp)
fildll 0x18(%rsp)
fdivrl 0xe65(%rip)
fldt   0x24(%rsp)
faddp  %st,%st(1)
fstpt  (%rsp)
call   1040 <sqrtl@plt>
cmp    %rax,%r13
jg     1180 <main+0x60>
lea    0x0(,%r13,4),%r9
neg    %r9
mov    %r13,%r8
add    %r13,%r8
mov    $0xffffffffffffffff,%rsi
mov    %r12,%rdi
mov    %r13,%rbx
jmp    120e <main+0xee>
cs nopw 0x0(%rax,%rax,1)
lea    (%rdi,%r13,1),%rdx
imul   %r9,%rdx
add    %r8,%rdi
lea    0x0(,%rbx,4),%rcx
or     $0x4,%rcx
imul   %r12,%rcx
add    %rcx,%rsi
add    %rdx,%rsi
add    $0x2,%rbx
cmp    %rax,%rbx
jg     1180 <main+0x60>
test   %rsi,%rsi
jle    1230 <main+0x110>
lea    (%rdi,%rdi,1),%rdx
not    %rdx
nopw   0x0(%rax,%rax,1)
add    %rdx,%rsi
add    $0x1,%rdi
add    $0xfffffffffffffffe,%rdx
test   %rsi,%rsi
jg     1220 <main+0x100>
mov    %edi,%edx
and    $0x1,%edx
or     %rsi,%rdx
jne    11e0 <main+0xc0>
mov    %rdi,%rdx
shr    $0x3f,%rdx
add    %rdi,%rdx
sar    %rdx
add    %rdx,%r15
jmp    11e0 <main+0xc0>

LLVM IR (-O0):
@.str = [6 x i8] c"%lld\0A\00"
define i32 @main() {
%1 = alloca i32
%2 = alloca i64
%3 = alloca i64
%4 = alloca i64
%5 = alloca i64
%6 = alloca i64
%7 = alloca i64
%8 = alloca i64
store i32 0, i32* %1
%9 = call i64 @sqrtl(x86_fp80 0xK40219502F8FFE0000000)
store i64 %9, i64* %3
store i64 0, i64* %4
store i64 2, i64* %2
br label %10
10: ; preds = %86, %0
%11 = load i64, i64* %2
%12 = load i64, i64* %3
%13 = icmp sle i64 %11, %12
br i1 %13, label %14, label %89
14: ; preds = %10
%15 = load i64, i64* %2
%16 = load i64, i64* %2
%17 = mul nsw i64 %15, %16
%18 = add nsw i64 1, %17
%19 = sitofp i64 %18 to x86_fp80
%20 = fdiv x86_fp80 0xK4043AD78EBC5AC620000, %19
%21 = fsub x86_fp80 %20, 0xK3FFF8000000000000000
%22 = call i64 @sqrtl(x86_fp80 %21)
store i64 %22, i64* %6
%23 = load i64, i64* %2
%24 = load i64, i64* %2
%25 = mul nsw i64 %23, %24
%26 = add nsw i64 1, %25
store i64 %26, i64* %7
store i64 -1, i64* %8
%27 = load i64, i64* %2
store i64 %27, i64* %5
br label %28
28: ; preds = %72, %14
%29 = load i64, i64* %5
%30 = load i64, i64* %6
%31 = icmp sle i64 %29, %30
br i1 %31, label %32, label %85
32: ; preds = %28
br label %33
33: ; preds = %36, %32
%34 = load i64, i64* %8
%35 = icmp sgt i64 %34, 0
br i1 %35, label %36, label %44
36: ; preds = %33
%37 = load i64, i64* %7
%38 = mul nsw i64 2, %37
%39 = add nsw i64 %38, 1
%40 = load i64, i64* %8
%41 = sub nsw i64 %40, %39
store i64 %41, i64* %8
%42 = load i64, i64* %7
%43 = add nsw i64 %42, 1
store i64 %43, i64* %7
br label %33
44: ; preds = %33
%45 = load i64, i64* %8
%46 = icmp eq i64 %45, 0
br i1 %46, label %47, label %56
47: ; preds = %44
%48 = load i64, i64* %7
%49 = srem i64 %48, 2
%50 = icmp eq i64 %49, 0
br i1 %50, label %51, label %56
51: ; preds = %47
%52 = load i64, i64* %7
%53 = sdiv i64 %52, 2
%54 = load i64, i64* %4
%55 = add nsw i64 %54, %53
store i64 %55, i64* %4
br label %56
56: ; preds = %51, %47, %44
%57 = load i64, i64* %2
%58 = mul nsw i64 4, %57
%59 = load i64, i64* %7
%60 = mul nsw i64 %58, %59
%61 = load i64, i64* %2
%62 = mul nsw i64 4, %61
%63 = load i64, i64* %2
%64 = mul nsw i64 %62, %63
%65 = add nsw i64 %60, %64
%66 = load i64, i64* %8
%67 = sub nsw i64 %66, %65
store i64 %67, i64* %8
%68 = load i64, i64* %2
%69 = mul nsw i64 2, %68
%70 = load i64, i64* %7
%71 = add nsw i64 %70, %69
store i64 %71, i64* %7
br label %72
72: ; preds = %56
%73 = load i64, i64* %5
%74 = mul nsw i64 4, %73
%75 = add nsw i64 %74, 4
%76 = load i64, i64* %2
%77 = load i64, i64* %2
%78 = mul nsw i64 %76, %77
%79 = add nsw i64 1, %78
%80 = mul nsw i64 %75, %79
%81 = load i64, i64* %8
%82 = add nsw i64 %81, %80
store i64 %82, i64* %8
%83 = load i64, i64* %5
%84 = add nsw i64 %83, 2
store i64 %84, i64* %5
br label %28
85: ; preds = %28
br label %86
86: ; preds = %85
%87 = load i64, i64* %2
%88 = add nsw i64 %87, 2
store i64 %88, i64* %2
br label %10
89: ; preds = %10
%90 = load i64, i64* %4
%91 = call i32 @printf(i8* getelementptr inbounds ([6 x i8], [6 x i8]* @.str, i64 0, i64 0), i64 %90)
%92 = load i32, i32* %1
ret i32 %92
}

LLVM IR (-O2):
@.str = [6 x i8] c"%lld\0A\00"
define i32 @main() {
%1 = tail call i64 @sqrtl(x86_fp80 0xK40219502F8FFE0000000)
%2 = icmp slt i64 %1, 2
br i1 %2, label %56, label %3
3: ; preds = %0, %52
%4 = phi i64 [ %54, %52 ], [ 2, %0 ]
%5 = phi i64 [ %53, %52 ], [ 0, %0 ]
%6 = mul nsw i64 %4, %4
%7 = or i64 %6, 1
%8 = sitofp i64 %7 to x86_fp80
%9 = fdiv x86_fp80 0xK4043AD78EBC5AC620000, %8
%10 = fadd x86_fp80 %9, 0xKBFFF8000000000000000
%11 = tail call i64 @sqrtl(x86_fp80 %10)
%12 = icmp sgt i64 %4, %11
br i1 %12, label %52, label %13
13: ; preds = %3
%14 = mul i64 %4, -4
%15 = shl nuw nsw i64 %4, 1
br label %16
16: ; preds = %13, %40
%17 = phi i64 [ -1, %13 ], [ %49, %40 ]
%18 = phi i64 [ %7, %13 ], [ %44, %40 ]
%19 = phi i64 [ %4, %13 ], [ %50, %40 ]
%20 = phi i64 [ %5, %13 ], [ %41, %40 ]
%21 = icmp sgt i64 %17, 0
br i1 %21, label %22, label %30
22: ; preds = %16, %22
%23 = phi i64 [ %27, %22 ], [ %17, %16 ]
%24 = phi i64 [ %28, %22 ], [ %18, %16 ]
%25 = shl nsw i64 %24, 1
%26 = xor i64 %25, -1
%27 = add i64 %23, %26
%28 = add nsw i64 %24, 1
%29 = icmp sgt i64 %27, 0
br i1 %29, label %22, label %30
30: ; preds = %22, %16
%31 = phi i64 [ %18, %16 ], [ %28, %22 ]
%32 = phi i64 [ %17, %16 ], [ %27, %22 ]
%33 = icmp eq i64 %32, 0
%34 = and i64 %31, 1
%35 = icmp eq i64 %34, 0
%36 = select i1 %33, i1 %35, i1 false
br i1 %36, label %37, label %40
37: ; preds = %30
%38 = sdiv i64 %31, 2
%39 = add nsw i64 %38, %20
br label %40
40: ; preds = %37, %30
%41 = phi i64 [ %39, %37 ], [ %20, %30 ]
%42 = add i64 %31, %4
%43 = mul i64 %14, %42
%44 = add nsw i64 %31, %15
%45 = shl nsw i64 %19, 2
%46 = or i64 %45, 4
%47 = mul nsw i64 %46, %7
%48 = add i64 %32, %47
%49 = add i64 %48, %43
%50 = add nuw nsw i64 %19, 2
%51 = icmp sgt i64 %50, %11
br i1 %51, label %52, label %16
52: ; preds = %40, %3
%53 = phi i64 [ %5, %3 ], [ %41, %40 ]
%54 = add nuw nsw i64 %4, 2
%55 = icmp sgt i64 %54, %1
br i1 %55, label %56, label %3
56: ; preds = %52, %0
%57 = phi i64 [ 0, %0 ], [ %53, %52 ]
%58 = tail call i32 @printf(i8* getelementptr inbounds ([6 x i8], [6 x i8]* @.str, i64 0, i64 0), i64 %57)
ret i32 0
}