! permとfactの前計算は ~/repos/contest/atcoder_ojt/abc156_e/fortran.f08 参照
program main
  implicit none
  integer(8) :: n, k, i, res
  integer(8), parameter :: m = 10 ** 9 + 7
  integer(8), allocatable :: a(:)
  integer(8), allocatable :: fact(:), ifact(:)
  read *, n, k
  allocate (a(n), fact(0:n), ifact(0:n))
  read *, a
  call quicksort(a)
  call calc_fact

  res = 0

  ! Add sum(max(s)) to res
  do i = 1, n
    res = res + modulo(a(i) * combi(i - 1, k - 1), m)
    res = modulo(res, m)
  end do

  ! Subtract sum(min(s)) from res
  do i = 1, n
    res = res - modulo(a(i) * combi(n - i, k - 1), m)
    res = modulo(res, m)
  end do

  print "(i0)", res
  contains

  recursive subroutine quicksort(a)
    integer(8), intent(inout) :: a(:)
    integer(8) :: pivot, i, j, k, tmp, n
    real(8) :: r
    n = size(a)
    if (n <= 1) return
    i = 1
    j = n
    call random_number(r)
    k = 1 + floor(r * n)
    k = max(1_8, min(k, n))
    pivot = a(k)
    do
      do while (.not. pivot <= a(i))
        i = i + 1
      end do
      do while (.not. a(j) <= pivot)
        j = j - 1
      end do
      if (i >= j) exit
      tmp = a(i)
      a(i) = a(j)
      a(j) = tmp
      i = i + 1
      j = j - 1
    end do
    call quicksort(a(1:i - 1))
    call quicksort(a(j + 1:n))
  end subroutine quicksort

  integer(8) function combi(n, k)
    integer(8), intent(in) :: n, k
    if (n < k) then
      combi = 0
      return
    end if
    combi = fact(n)
    combi = combi * ifact(k)
    combi = modulo(combi, m)
    combi = combi * ifact(n - k)
    combi = modulo(combi, m)
  end function combi

  subroutine calc_fact()
    integer(8) :: i
    fact(0) = 1
    ifact(0) = inv(1_8)
    do i = 1, n
      fact(i) = modulo(fact(i - 1) * i, m)
      ifact(i) = inv(fact(i))
    end do
  end subroutine calc_fact

  function inv(n)
    integer(8), intent(in) :: n
    integer(8) :: inv
    inv = pow(n, m - 2)
  end function inv

  recursive function pow(n, p) result(res)
    integer(8), intent(in) :: n, p
    integer(8) :: res
    if (p == 0) then
      res = 1
    else if (modulo(p, 2_8) == 0) then
      res = modulo(pow(n, p / 2) ** 2, m)
    else
      res = modulo(pow(n, p - 1) * n, m)
    end if
  end function pow
end program main

; ModuleID = 'FIRModule'
source_filename = "FIRModule"
target datalayout = "e-m:e-p270:32:32-p271:32:32-p272:64:64-i64:64-f80:128-n8:16:32:64-S128"
target triple = "x86_64-pc-linux-gnu"

@_QFEa.0 = internal unnamed_addr global ptr null
@_QFEa.8 = internal unnamed_addr global i64 0
@_QFEfact.0 = internal unnamed_addr global ptr null
@_QFEfact.7 = internal unnamed_addr global i64 1
@_QFEifact.0 = internal unnamed_addr global ptr null
@_QFEifact.7 = internal unnamed_addr global i64 1
@_QQcl.1caebacab38ce13f4cbbbce7bfba5aa8 = linkonce constant [73 x i8] c"/home/ubuntu/Downloads/TheStack_IR/Collated/Fortran/Source_2646/source.f\00"
@_QQcl.28693029 = linkonce constant [4 x i8] c"(i0)"
@_QQEnvironmentDefaults = local_unnamed_addr constant ptr null

; Function Attrs: mustprogress nofree nounwind willreturn allockind("alloc,uninitialized") allocsize(0) memory(inaccessiblemem: readwrite)
declare noalias noundef ptr @malloc(i64 noundef) local_unnamed_addr #0

define void @_QQmain() local_unnamed_addr {
  %1 = alloca i64, align 8
  %2 = alloca i64, align 8
  %3 = alloca i64, align 8
  %4 = alloca { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, align 8
  %5 = alloca { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, align 8
  %6 = alloca i64, align 8
  %7 = alloca i64, align 8
  %8 = alloca { ptr }, align 8
  store ptr %7, ptr %8, align 8, !tbaa !3
  %9 = call ptr @_FortranAioBeginExternalListInput(i32 -1, ptr nonnull @_QQcl.1caebacab38ce13f4cbbbce7bfba5aa8, i32 8)
  %10 = call i1 @_FortranAioInputInteger(ptr %9, ptr nonnull %7, i32 8)
  %11 = call i1 @_FortranAioInputInteger(ptr %9, ptr nonnull %6, i32 8)
  %12 = call i32 @_FortranAioEndIoStatement(ptr %9)
  %13 = load i64, ptr %7, align 8, !tbaa !3
  %14 = call i64 @llvm.smax.i64(i64 %13, i64 0)
  %15 = shl i64 %14, 3
  %16 = call ptr @malloc(i64 %15)
  store ptr %16, ptr @_QFEa.0, align 8, !tbaa !7
  store i64 %14, ptr @_QFEa.8, align 8, !tbaa !7
  %17 = add i64 %13, 1
  %18 = call i64 @llvm.smax.i64(i64 %17, i64 0)
  %19 = shl i64 %18, 3
  %20 = call ptr @malloc(i64 %19)
  %21 = icmp ugt i64 %13, 9223372036854775806
  %22 = zext i1 %21 to i64
  store ptr %20, ptr @_QFEfact.0, align 8, !tbaa !7
  store i64 %22, ptr @_QFEfact.7, align 8, !tbaa !7
  %23 = call ptr @malloc(i64 %19)
  store ptr %23, ptr @_QFEifact.0, align 8, !tbaa !7
  store i64 %22, ptr @_QFEifact.7, align 8, !tbaa !7
  %24 = call ptr @_FortranAioBeginExternalListInput(i32 -1, ptr nonnull @_QQcl.1caebacab38ce13f4cbbbce7bfba5aa8, i32 10)
  %.unpack = load ptr, ptr @_QFEa.0, align 8, !tbaa !7
  %.unpack271.unpack.unpack273 = load i64, ptr @_QFEa.8, align 8, !tbaa !7
  store ptr %.unpack, ptr %5, align 8, !tbaa !7
  %.fca.1.gep = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %5, i64 0, i32 1
  store i64 8, ptr %.fca.1.gep, align 8, !tbaa !7
  %.fca.2.gep = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %5, i64 0, i32 2
  store i32 20180515, ptr %.fca.2.gep, align 8, !tbaa !7
  %.fca.3.gep = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %5, i64 0, i32 3
  store i8 1, ptr %.fca.3.gep, align 4, !tbaa !7
  %.fca.4.gep = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %5, i64 0, i32 4
  store i8 10, ptr %.fca.4.gep, align 1, !tbaa !7
  %.fca.5.gep = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %5, i64 0, i32 5
  store i8 0, ptr %.fca.5.gep, align 2, !tbaa !7
  %.fca.6.gep = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %5, i64 0, i32 6
  store i8 0, ptr %.fca.6.gep, align 1, !tbaa !7
  %.fca.7.0.0.gep = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %5, i64 0, i32 7, i64 0, i64 0
  store i64 1, ptr %.fca.7.0.0.gep, align 8, !tbaa !7
  %.fca.7.0.1.gep = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %5, i64 0, i32 7, i64 0, i64 1
  store i64 %.unpack271.unpack.unpack273, ptr %.fca.7.0.1.gep, align 8, !tbaa !7
  %.fca.7.0.2.gep = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %5, i64 0, i32 7, i64 0, i64 2
  store i64 8, ptr %.fca.7.0.2.gep, align 8, !tbaa !7
  %25 = call i1 @_FortranAioInputDescriptor(ptr %24, ptr nonnull %5)
  %26 = call i32 @_FortranAioEndIoStatement(ptr %24)
  %.unpack276 = load ptr, ptr @_QFEa.0, align 8, !tbaa !7
  %.unpack283.unpack.unpack285 = load i64, ptr @_QFEa.8, align 8, !tbaa !7
  store ptr %.unpack276, ptr %4, align 8, !tbaa !7
  %.fca.1.gep159 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %4, i64 0, i32 1
  store i64 8, ptr %.fca.1.gep159, align 8, !tbaa !7
  %.fca.2.gep161 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %4, i64 0, i32 2
  store i32 20180515, ptr %.fca.2.gep161, align 8, !tbaa !7
  %.fca.3.gep163 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %4, i64 0, i32 3
  store i8 1, ptr %.fca.3.gep163, align 4, !tbaa !7
  %.fca.4.gep165 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %4, i64 0, i32 4
  store i8 10, ptr %.fca.4.gep165, align 1, !tbaa !7
  %.fca.5.gep167 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %4, i64 0, i32 5
  store i8 0, ptr %.fca.5.gep167, align 2, !tbaa !7
  %.fca.6.gep169 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %4, i64 0, i32 6
  store i8 0, ptr %.fca.6.gep169, align 1, !tbaa !7
  %.fca.7.0.0.gep171 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %4, i64 0, i32 7, i64 0, i64 0
  store i64 1, ptr %.fca.7.0.0.gep171, align 8, !tbaa !7
  %.fca.7.0.1.gep173 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %4, i64 0, i32 7, i64 0, i64 1
  store i64 %.unpack283.unpack.unpack285, ptr %.fca.7.0.1.gep173, align 8, !tbaa !7
  %.fca.7.0.2.gep175 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %4, i64 0, i32 7, i64 0, i64 2
  store i64 8, ptr %.fca.7.0.2.gep175, align 8, !tbaa !7
  call void @_QFPquicksort(ptr nonnull %4, ptr nonnull %8)
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %3)
  %27 = load ptr, ptr %8, align 8, !tbaa !3
  %.unpack.i = load ptr, ptr @_QFEfact.0, align 8, !tbaa !7
  %.unpack90.unpack.unpack.i = load i64, ptr @_QFEfact.7, align 8, !tbaa !7
  %28 = sub i64 0, %.unpack90.unpack.unpack.i
  %29 = getelementptr i64, ptr %.unpack.i, i64 %28
  store i64 1, ptr %29, align 8, !tbaa !3
  store i64 1, ptr %3, align 8, !tbaa !3
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %2)
  store i64 1000000005, ptr %2, align 8, !tbaa !3
  %30 = call i64 @_QFPpow(ptr nonnull %3, ptr nonnull %2, ptr nonnull %8)
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %2)
  %.unpack95.i = load ptr, ptr @_QFEifact.0, align 8, !tbaa !7
  %.unpack102.unpack.unpack.i = load i64, ptr @_QFEifact.7, align 8, !tbaa !7
  %31 = sub i64 0, %.unpack102.unpack.unpack.i
  %32 = getelementptr i64, ptr %.unpack95.i, i64 %31
  store i64 %30, ptr %32, align 8, !tbaa !3
  %33 = load i64, ptr %27, align 8, !tbaa !3
  %34 = icmp sgt i64 %33, 0
  br i1 %34, label %.lr.ph.i, label %_QFPcalc_fact.exit

.lr.ph.i:                                         ; preds = %0
  %35 = xor i64 %.unpack90.unpack.unpack.i, -1
  br label %36

36:                                               ; preds = %36, %.lr.ph.i
  %37 = phi i64 [ 1, %.lr.ph.i ], [ %53, %36 ]
  %38 = add i64 %37, %35
  %39 = getelementptr i64, ptr %.unpack.i, i64 %38
  %40 = load i64, ptr %39, align 8, !tbaa !3
  %41 = mul i64 %40, %37
  %42 = srem i64 %41, 1000000007
  %43 = icmp slt i64 %41, 0
  %44 = icmp ne i64 %42, 0
  %45 = and i1 %43, %44
  %46 = add nsw i64 %42, 1000000007
  %47 = select i1 %45, i64 %46, i64 %42
  %48 = sub i64 %37, %.unpack90.unpack.unpack.i
  %49 = getelementptr i64, ptr %.unpack.i, i64 %48
  store i64 %47, ptr %49, align 8, !tbaa !3
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %1)
  store i64 1000000005, ptr %1, align 8, !tbaa !3
  %50 = call i64 @_QFPpow(ptr nonnull %49, ptr nonnull %1, ptr nonnull %8)
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %1)
  %51 = sub i64 %37, %.unpack102.unpack.unpack.i
  %52 = getelementptr i64, ptr %.unpack95.i, i64 %51
  store i64 %50, ptr %52, align 8, !tbaa !3
  %53 = add nuw i64 %37, 1
  %exitcond.not.i = icmp eq i64 %37, %33
  br i1 %exitcond.not.i, label %_QFPcalc_fact.exit, label %36

_QFPcalc_fact.exit:                               ; preds = %36, %0
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %3)
  %54 = load i64, ptr %7, align 8, !tbaa !3
  %55 = icmp sgt i64 %54, 0
  br i1 %55, label %.lr.ph, label %._crit_edge

.lr.ph:                                           ; preds = %_QFPcalc_fact.exit
  %.unpack300 = load ptr, ptr @_QFEa.0, align 8, !tbaa !7
  %56 = load i64, ptr %6, align 8, !tbaa !3
  %57 = add i64 %56, -1
  %58 = sub i64 %57, %.unpack102.unpack.unpack.i
  %59 = getelementptr i64, ptr %.unpack95.i, i64 %58
  %.neg323 = sub i64 1, %56
  br label %64

.preheader:                                       ; preds = %_QFPcombi.exit
  br i1 %55, label %.lr.ph327, label %._crit_edge

.lr.ph327:                                        ; preds = %.preheader
  %.unpack288 = load ptr, ptr @_QFEa.0, align 8, !tbaa !7
  %60 = load i64, ptr %6, align 8, !tbaa !3
  %61 = add i64 %60, -1
  %62 = sub i64 %61, %.unpack102.unpack.unpack.i
  %63 = getelementptr i64, ptr %.unpack95.i, i64 %62
  %.neg = sub i64 1, %60
  br label %107

64:                                               ; preds = %.lr.ph, %_QFPcombi.exit
  %65 = phi i64 [ 1, %.lr.ph ], [ %106, %_QFPcombi.exit ]
  %.0325 = phi i64 [ 0, %.lr.ph ], [ %105, %_QFPcombi.exit ]
  %66 = add i64 %65, -1
  %67 = getelementptr i64, ptr %.unpack300, i64 %66
  %68 = load i64, ptr %67, align 8, !tbaa !3
  %.not = icmp sgt i64 %65, %57
  br i1 %.not, label %69, label %_QFPcombi.exit

69:                                               ; preds = %64
  %70 = add nsw i64 %65, -1
  %71 = sub i64 %70, %.unpack90.unpack.unpack.i
  %72 = getelementptr i64, ptr %.unpack.i, i64 %71
  %73 = load i64, ptr %72, align 8, !tbaa !3
  %74 = load i64, ptr %59, align 8, !tbaa !3
  %75 = mul i64 %74, %73
  %76 = srem i64 %75, 1000000007
  %77 = icmp slt i64 %75, 0
  %78 = icmp ne i64 %76, 0
  %79 = and i1 %77, %78
  %80 = add nsw i64 %76, 1000000007
  %81 = select i1 %79, i64 %80, i64 %76
  %.neg324 = sub i64 %70, %.unpack102.unpack.unpack.i
  %82 = add i64 %.neg324, %.neg323
  %83 = getelementptr i64, ptr %.unpack95.i, i64 %82
  %84 = load i64, ptr %83, align 8, !tbaa !3
  %85 = mul i64 %81, %84
  %86 = srem i64 %85, 1000000007
  %87 = icmp slt i64 %85, 0
  %88 = icmp ne i64 %86, 0
  %89 = and i1 %87, %88
  %90 = add nsw i64 %86, 1000000007
  %91 = select i1 %89, i64 %90, i64 %86
  br label %_QFPcombi.exit

_QFPcombi.exit:                                   ; preds = %64, %69
  %.0.i = phi i64 [ %91, %69 ], [ 0, %64 ]
  %92 = mul i64 %.0.i, %68
  %93 = srem i64 %92, 1000000007
  %94 = icmp slt i64 %92, 0
  %95 = icmp ne i64 %93, 0
  %96 = and i1 %94, %95
  %97 = add nsw i64 %93, 1000000007
  %98 = select i1 %96, i64 %97, i64 %93
  %99 = add nsw i64 %98, %.0325
  %100 = srem i64 %99, 1000000007
  %101 = icmp slt i64 %99, 0
  %102 = icmp ne i64 %100, 0
  %103 = and i1 %101, %102
  %104 = add nsw i64 %100, 1000000007
  %105 = select i1 %103, i64 %104, i64 %100
  %106 = add nuw i64 %65, 1
  %exitcond.not = icmp eq i64 %65, %54
  br i1 %exitcond.not, label %.preheader, label %64

107:                                              ; preds = %.lr.ph327, %_QFPcombi.exit318
  %108 = phi i64 [ 1, %.lr.ph327 ], [ %150, %_QFPcombi.exit318 ]
  %.1326 = phi i64 [ %105, %.lr.ph327 ], [ %149, %_QFPcombi.exit318 ]
  %109 = add i64 %108, -1
  %110 = getelementptr i64, ptr %.unpack288, i64 %109
  %111 = load i64, ptr %110, align 8, !tbaa !3
  %112 = sub nsw i64 %54, %108
  %113 = icmp slt i64 %112, %61
  br i1 %113, label %_QFPcombi.exit318, label %114

114:                                              ; preds = %107
  %115 = sub i64 %112, %.unpack90.unpack.unpack.i
  %116 = getelementptr i64, ptr %.unpack.i, i64 %115
  %117 = load i64, ptr %116, align 8, !tbaa !3
  %118 = load i64, ptr %63, align 8, !tbaa !3
  %119 = mul i64 %118, %117
  %120 = srem i64 %119, 1000000007
  %121 = icmp slt i64 %119, 0
  %122 = icmp ne i64 %120, 0
  %123 = and i1 %121, %122
  %124 = add nsw i64 %120, 1000000007
  %125 = select i1 %123, i64 %124, i64 %120
  %.neg322 = sub i64 %112, %.unpack102.unpack.unpack.i
  %126 = add i64 %.neg322, %.neg
  %127 = getelementptr i64, ptr %.unpack95.i, i64 %126
  %128 = load i64, ptr %127, align 8, !tbaa !3
  %129 = mul i64 %125, %128
  %130 = srem i64 %129, 1000000007
  %131 = icmp slt i64 %129, 0
  %132 = icmp ne i64 %130, 0
  %133 = and i1 %131, %132
  %134 = add nsw i64 %130, 1000000007
  %135 = select i1 %133, i64 %134, i64 %130
  br label %_QFPcombi.exit318

_QFPcombi.exit318:                                ; preds = %107, %114
  %.0.i317 = phi i64 [ %135, %114 ], [ 0, %107 ]
  %136 = mul i64 %.0.i317, %111
  %137 = srem i64 %136, 1000000007
  %138 = icmp slt i64 %136, 0
  %139 = icmp ne i64 %137, 0
  %140 = and i1 %138, %139
  %141 = add nsw i64 %137, 1000000007
  %142 = select i1 %140, i64 %141, i64 %137
  %143 = sub nsw i64 %.1326, %142
  %144 = srem i64 %143, 1000000007
  %145 = icmp slt i64 %143, 0
  %146 = icmp ne i64 %144, 0
  %147 = and i1 %145, %146
  %148 = add nsw i64 %144, 1000000007
  %149 = select i1 %147, i64 %148, i64 %144
  %150 = add nuw i64 %108, 1
  %exitcond330.not = icmp eq i64 %54, %108
  br i1 %exitcond330.not, label %._crit_edge, label %107

._crit_edge:                                      ; preds = %_QFPcombi.exit318, %_QFPcalc_fact.exit, %.preheader
  %.1.lcssa = phi i64 [ %105, %.preheader ], [ 0, %_QFPcalc_fact.exit ], [ %149, %_QFPcombi.exit318 ]
  %151 = call ptr @_FortranAioBeginExternalFormattedOutput(ptr nonnull @_QQcl.28693029, i64 4, ptr null, i32 -1, ptr nonnull @_QQcl.1caebacab38ce13f4cbbbce7bfba5aa8, i32 28)
  %152 = call i1 @_FortranAioOutputInteger64(ptr %151, i64 %.1.lcssa)
  %153 = call i32 @_FortranAioEndIoStatement(ptr %151)
  ret void
}

define void @_QFPquicksort(ptr nocapture readonly %0, ptr nest %1) local_unnamed_addr {
  %3 = alloca { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, align 8
  %4 = alloca { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, align 8
  %5 = alloca { ptr, i64, i32, i8, i8, i8, i8 }, align 8
  %6 = alloca double, align 8
  %7 = getelementptr { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %0, i64 0, i32 7, i64 0, i64 1
  %8 = load i64, ptr %7, align 8, !tbaa !7
  %9 = icmp slt i64 %8, 2
  br i1 %9, label %53, label %10

10:                                               ; preds = %2
  store ptr %6, ptr %5, align 8, !tbaa !7
  %.fca.1.gep = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8 }, ptr %5, i64 0, i32 1
  store i64 8, ptr %.fca.1.gep, align 8, !tbaa !7
  %.fca.2.gep = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8 }, ptr %5, i64 0, i32 2
  store i32 20180515, ptr %.fca.2.gep, align 8, !tbaa !7
  %.fca.3.gep = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8 }, ptr %5, i64 0, i32 3
  store i8 0, ptr %.fca.3.gep, align 4, !tbaa !7
  %.fca.4.gep = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8 }, ptr %5, i64 0, i32 4
  store i8 28, ptr %.fca.4.gep, align 1, !tbaa !7
  %.fca.5.gep = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8 }, ptr %5, i64 0, i32 5
  store i8 0, ptr %.fca.5.gep, align 2, !tbaa !7
  %.fca.6.gep = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8 }, ptr %5, i64 0, i32 6
  store i8 0, ptr %.fca.6.gep, align 1, !tbaa !7
  %11 = call {} @_FortranARandomNumber(ptr nonnull %5, ptr nonnull @_QQcl.1caebacab38ce13f4cbbbce7bfba5aa8, i32 39)
  %12 = load double, ptr %6, align 8, !tbaa !3
  %13 = sitofp i64 %8 to double
  %14 = fmul contract double %12, %13
  %15 = call contract double @llvm.floor.f64(double %14)
  %16 = fptosi double %15 to i32
  %17 = add i32 %16, 1
  %18 = sext i32 %17 to i64
  %19 = call i64 @llvm.smin.i64(i64 %8, i64 %18)
  %20 = call i64 @llvm.smax.i64(i64 %19, i64 1)
  %21 = add nsw i64 %20, -1
  %22 = load ptr, ptr %0, align 8, !tbaa !7
  %23 = getelementptr { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %0, i64 0, i32 7, i64 0, i64 2
  %24 = load i64, ptr %23, align 8, !tbaa !7
  %25 = mul i64 %21, %24
  %26 = getelementptr i8, ptr %22, i64 %25
  %27 = load i64, ptr %26, align 8, !tbaa !3
  br label %.outer

.outer:                                           ; preds = %51, %10
  %.058.ph = phi i64 [ %37, %51 ], [ %8, %10 ]
  %.0.ph = phi i64 [ %34, %51 ], [ 1, %10 ]
  br label %28

28:                                               ; preds = %28, %.outer
  %.0 = phi i64 [ %34, %28 ], [ %.0.ph, %.outer ]
  %29 = add i64 %.0, -1
  %30 = mul i64 %29, %24
  %31 = getelementptr i8, ptr %22, i64 %30
  %32 = load i64, ptr %31, align 8, !tbaa !3
  %33 = icmp sgt i64 %27, %32
  %34 = add i64 %.0, 1
  br i1 %33, label %28, label %.preheader

.preheader:                                       ; preds = %28
  %35 = getelementptr i8, ptr %22, i64 %30
  br label %36

36:                                               ; preds = %.preheader, %36
  %.1 = phi i64 [ %37, %36 ], [ %.058.ph, %.preheader ]
  %37 = add i64 %.1, -1
  %38 = mul i64 %37, %24
  %39 = getelementptr i8, ptr %22, i64 %38
  %40 = load i64, ptr %39, align 8, !tbaa !3
  %41 = icmp sgt i64 %40, %27
  br i1 %41, label %36, label %42

42:                                               ; preds = %36
  %.not = icmp slt i64 %.0, %.1
  br i1 %.not, label %51, label %43

43:                                               ; preds = %42
  %44 = call i64 @llvm.smax.i64(i64 %29, i64 0)
  store ptr %22, ptr %4, align 8, !tbaa !7
  %.fca.1.gep27 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %4, i64 0, i32 1
  store i64 8, ptr %.fca.1.gep27, align 8, !tbaa !7
  %.fca.2.gep29 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %4, i64 0, i32 2
  store i32 20180515, ptr %.fca.2.gep29, align 8, !tbaa !7
  %.fca.3.gep31 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %4, i64 0, i32 3
  store i8 1, ptr %.fca.3.gep31, align 4, !tbaa !7
  %.fca.4.gep33 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %4, i64 0, i32 4
  store i8 10, ptr %.fca.4.gep33, align 1, !tbaa !7
  %.fca.5.gep35 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %4, i64 0, i32 5
  store i8 0, ptr %.fca.5.gep35, align 2, !tbaa !7
  %.fca.6.gep37 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %4, i64 0, i32 6
  store i8 0, ptr %.fca.6.gep37, align 1, !tbaa !7
  %.fca.7.0.0.gep = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %4, i64 0, i32 7, i64 0, i64 0
  store i64 1, ptr %.fca.7.0.0.gep, align 8, !tbaa !7
  %.fca.7.0.1.gep = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %4, i64 0, i32 7, i64 0, i64 1
  store i64 %44, ptr %.fca.7.0.1.gep, align 8, !tbaa !7
  %.fca.7.0.2.gep = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %4, i64 0, i32 7, i64 0, i64 2
  store i64 %24, ptr %.fca.7.0.2.gep, align 8, !tbaa !7
  call void @_QFPquicksort(ptr nonnull %4, ptr %1)
  %45 = load i64, ptr %23, align 8, !tbaa !7
  %46 = load ptr, ptr %0, align 8, !tbaa !7
  %47 = mul i64 %45, %.1
  %48 = getelementptr i8, ptr %46, i64 %47
  %49 = sub i64 %8, %.1
  %50 = call i64 @llvm.smax.i64(i64 %49, i64 0)
  store ptr %48, ptr %3, align 8, !tbaa !7
  %.fca.1.gep41 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %3, i64 0, i32 1
  store i64 8, ptr %.fca.1.gep41, align 8, !tbaa !7
  %.fca.2.gep43 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %3, i64 0, i32 2
  store i32 20180515, ptr %.fca.2.gep43, align 8, !tbaa !7
  %.fca.3.gep45 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %3, i64 0, i32 3
  store i8 1, ptr %.fca.3.gep45, align 4, !tbaa !7
  %.fca.4.gep47 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %3, i64 0, i32 4
  store i8 10, ptr %.fca.4.gep47, align 1, !tbaa !7
  %.fca.5.gep49 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %3, i64 0, i32 5
  store i8 0, ptr %.fca.5.gep49, align 2, !tbaa !7
  %.fca.6.gep51 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %3, i64 0, i32 6
  store i8 0, ptr %.fca.6.gep51, align 1, !tbaa !7
  %.fca.7.0.0.gep53 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %3, i64 0, i32 7, i64 0, i64 0
  store i64 1, ptr %.fca.7.0.0.gep53, align 8, !tbaa !7
  %.fca.7.0.1.gep55 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %3, i64 0, i32 7, i64 0, i64 1
  store i64 %50, ptr %.fca.7.0.1.gep55, align 8, !tbaa !7
  %.fca.7.0.2.gep57 = getelementptr inbounds { ptr, i64, i32, i8, i8, i8, i8, [1 x [3 x i64]] }, ptr %3, i64 0, i32 7, i64 0, i64 2
  store i64 %45, ptr %.fca.7.0.2.gep57, align 8, !tbaa !7
  call void @_QFPquicksort(ptr nonnull %3, ptr %1)
  br label %53

51:                                               ; preds = %42
  %52 = getelementptr i8, ptr %22, i64 %38
  store i64 %40, ptr %35, align 8, !tbaa !3
  store i64 %32, ptr %52, align 8, !tbaa !3
  br label %.outer

53:                                               ; preds = %43, %2
  ret void
}

; Function Attrs: mustprogress nofree norecurse nosync nounwind willreturn memory(read, inaccessiblemem: none)
define i64 @_QFPcombi(ptr nocapture readonly %0, ptr nocapture readonly %1, ptr nest nocapture readnone %2) local_unnamed_addr #1 {
  %4 = load i64, ptr %0, align 8, !tbaa !3
  %5 = load i64, ptr %1, align 8, !tbaa !3
  %6 = icmp slt i64 %4, %5
  br i1 %6, label %32, label %7

7:                                                ; preds = %3
  %.unpack = load ptr, ptr @_QFEfact.0, align 8, !tbaa !7
  %.unpack53.unpack.unpack = load i64, ptr @_QFEfact.7, align 8, !tbaa !7
  %8 = sub i64 %4, %.unpack53.unpack.unpack
  %9 = getelementptr i64, ptr %.unpack, i64 %8
  %10 = load i64, ptr %9, align 8, !tbaa !3
  %.unpack58 = load ptr, ptr @_QFEifact.0, align 8, !tbaa !7
  %.unpack65.unpack.unpack = load i64, ptr @_QFEifact.7, align 8, !tbaa !7
  %11 = sub i64 %5, %.unpack65.unpack.unpack
  %12 = getelementptr i64, ptr %.unpack58, i64 %11
  %13 = load i64, ptr %12, align 8, !tbaa !3
  %14 = mul i64 %13, %10
  %15 = srem i64 %14, 1000000007
  %16 = icmp slt i64 %14, 0
  %17 = icmp ne i64 %15, 0
  %18 = and i1 %16, %17
  %19 = add nsw i64 %15, 1000000007
  %20 = select i1 %18, i64 %19, i64 %15
  %21 = add i64 %5, %.unpack65.unpack.unpack
  %22 = sub i64 %4, %21
  %23 = getelementptr i64, ptr %.unpack58, i64 %22
  %24 = load i64, ptr %23, align 8, !tbaa !3
  %25 = mul i64 %20, %24
  %26 = srem i64 %25, 1000000007
  %27 = icmp slt i64 %25, 0
  %28 = icmp ne i64 %26, 0
  %29 = and i1 %27, %28
  %30 = add nsw i64 %26, 1000000007
  %31 = select i1 %29, i64 %30, i64 %26
  br label %32

32:                                               ; preds = %3, %7
  %.0 = phi i64 [ %31, %7 ], [ 0, %3 ]
  ret i64 %.0
}

; Function Attrs: nofree nosync nounwind memory(readwrite, inaccessiblemem: none)
define void @_QFPcalc_fact(ptr nest %0) local_unnamed_addr #2 {
  %2 = alloca i64, align 8
  %3 = alloca i64, align 8
  %4 = alloca i64, align 8
  %5 = load ptr, ptr %0, align 8, !tbaa !3
  %.unpack = load ptr, ptr @_QFEfact.0, align 8, !tbaa !7
  %.unpack90.unpack.unpack = load i64, ptr @_QFEfact.7, align 8, !tbaa !7
  %6 = sub i64 0, %.unpack90.unpack.unpack
  %7 = getelementptr i64, ptr %.unpack, i64 %6
  store i64 1, ptr %7, align 8, !tbaa !3
  store i64 1, ptr %4, align 8, !tbaa !3
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %3)
  store i64 1000000005, ptr %3, align 8, !tbaa !3
  %8 = call i64 @_QFPpow(ptr nonnull %4, ptr nonnull %3, ptr nonnull %0)
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %3)
  %.unpack95 = load ptr, ptr @_QFEifact.0, align 8, !tbaa !7
  %.unpack102.unpack.unpack = load i64, ptr @_QFEifact.7, align 8, !tbaa !7
  %9 = sub i64 0, %.unpack102.unpack.unpack
  %10 = getelementptr i64, ptr %.unpack95, i64 %9
  store i64 %8, ptr %10, align 8, !tbaa !3
  %11 = load i64, ptr %5, align 8, !tbaa !3
  %12 = icmp sgt i64 %11, 0
  br i1 %12, label %.lr.ph, label %._crit_edge

.lr.ph:                                           ; preds = %1
  %13 = xor i64 %.unpack90.unpack.unpack, -1
  br label %14

14:                                               ; preds = %.lr.ph, %14
  %15 = phi i64 [ 1, %.lr.ph ], [ %31, %14 ]
  %16 = add i64 %15, %13
  %17 = getelementptr i64, ptr %.unpack, i64 %16
  %18 = load i64, ptr %17, align 8, !tbaa !3
  %19 = mul i64 %18, %15
  %20 = srem i64 %19, 1000000007
  %21 = icmp slt i64 %19, 0
  %22 = icmp ne i64 %20, 0
  %23 = and i1 %21, %22
  %24 = add nsw i64 %20, 1000000007
  %25 = select i1 %23, i64 %24, i64 %20
  %26 = sub i64 %15, %.unpack90.unpack.unpack
  %27 = getelementptr i64, ptr %.unpack, i64 %26
  store i64 %25, ptr %27, align 8, !tbaa !3
  call void @llvm.lifetime.start.p0(i64 8, ptr nonnull %2)
  store i64 1000000005, ptr %2, align 8, !tbaa !3
  %28 = call i64 @_QFPpow(ptr nonnull %27, ptr nonnull %2, ptr nonnull %0)
  call void @llvm.lifetime.end.p0(i64 8, ptr nonnull %2)
  %29 = sub i64 %15, %.unpack102.unpack.unpack
  %30 = getelementptr i64, ptr %.unpack95, i64 %29
  store i64 %28, ptr %30, align 8, !tbaa !3
  %31 = add nuw i64 %15, 1
  %exitcond.not = icmp eq i64 %15, %11
  br i1 %exitcond.not, label %._crit_edge, label %14

._crit_edge:                                      ; preds = %14, %1
  ret void
}

; Function Attrs: nofree nosync nounwind memory(argmem: read)
define i64 @_QFPinv(ptr readonly %0, ptr nest readonly %1) local_unnamed_addr #3 {
  %3 = alloca i64, align 8
  store i64 1000000005, ptr %3, align 8, !tbaa !3
  %4 = call i64 @_QFPpow(ptr %0, ptr nonnull %3, ptr %1)
  ret i64 %4
}

; Function Attrs: nofree nosync nounwind memory(argmem: read)
define i64 @_QFPpow(ptr %0, ptr nocapture readonly %1, ptr nest %2) local_unnamed_addr #3 {
  %4 = alloca i64, align 8
  %5 = alloca i64, align 8
  %6 = load i64, ptr %1, align 8, !tbaa !3
  %7 = icmp eq i64 %6, 0
  br i1 %7, label %37, label %8

8:                                                ; preds = %3
  %9 = srem i64 %6, 2
  %10 = sdiv i64 %6, 2
  %11 = icmp slt i64 %6, 0
  %12 = icmp ne i64 %9, 0
  %13 = and i1 %11, %12
  %14 = add nsw i64 %9, 2
  %15 = select i1 %13, i64 %14, i64 %9
  %16 = icmp eq i64 %15, 0
  br i1 %16, label %17, label %26

17:                                               ; preds = %8
  store i64 %10, ptr %5, align 8, !tbaa !3
  %18 = call i64 @_QFPpow(ptr %0, ptr nonnull %5, ptr %2)
  %19 = mul i64 %18, %18
  %20 = srem i64 %19, 1000000007
  %21 = icmp slt i64 %19, 0
  %22 = icmp ne i64 %20, 0
  %23 = and i1 %21, %22
  %24 = add nsw i64 %20, 1000000007
  %25 = select i1 %23, i64 %24, i64 %20
  br label %37

26:                                               ; preds = %8
  %27 = add i64 %6, -1
  store i64 %27, ptr %4, align 8, !tbaa !3
  %28 = call i64 @_QFPpow(ptr %0, ptr nonnull %4, ptr %2)
  %29 = load i64, ptr %0, align 8, !tbaa !3
  %30 = mul i64 %29, %28
  %31 = srem i64 %30, 1000000007
  %32 = icmp slt i64 %30, 0
  %33 = icmp ne i64 %31, 0
  %34 = and i1 %32, %33
  %35 = add nsw i64 %31, 1000000007
  %36 = select i1 %34, i64 %35, i64 %31
  br label %37

37:                                               ; preds = %3, %17, %26
  %.0 = phi i64 [ %25, %17 ], [ %36, %26 ], [ 1, %3 ]
  ret i64 %.0
}

declare ptr @_FortranAioBeginExternalListInput(i32, ptr, i32) local_unnamed_addr

declare i1 @_FortranAioInputInteger(ptr, ptr, i32) local_unnamed_addr

declare i32 @_FortranAioEndIoStatement(ptr) local_unnamed_addr

declare i1 @_FortranAioInputDescriptor(ptr, ptr) local_unnamed_addr

declare ptr @_FortranAioBeginExternalFormattedOutput(ptr, i64, ptr, i32, ptr, i32) local_unnamed_addr

declare i1 @_FortranAioOutputInteger64(ptr, i64) local_unnamed_addr

declare {} @_FortranARandomNumber(ptr, ptr, i32) local_unnamed_addr

; Function Attrs: mustprogress nocallback nofree nosync nounwind speculatable willreturn memory(none)
declare double @llvm.floor.f64(double) #4

; Function Attrs: nocallback nofree nosync nounwind speculatable willreturn memory(none)
declare i64 @llvm.smax.i64(i64, i64) #5

; Function Attrs: nocallback nofree nosync nounwind speculatable willreturn memory(none)
declare i64 @llvm.smin.i64(i64, i64) #5

; Function Attrs: nocallback nofree nosync nounwind willreturn memory(argmem: readwrite)
declare void @llvm.lifetime.start.p0(i64 immarg, ptr nocapture) #6

; Function Attrs: nocallback nofree nosync nounwind willreturn memory(argmem: readwrite)
declare void @llvm.lifetime.end.p0(i64 immarg, ptr nocapture) #6

attributes #0 = { mustprogress nofree nounwind willreturn allockind("alloc,uninitialized") allocsize(0) memory(inaccessiblemem: readwrite) "alloc-family"="malloc" }
attributes #1 = { mustprogress nofree norecurse nosync nounwind willreturn memory(read, inaccessiblemem: none) }
attributes #2 = { nofree nosync nounwind memory(readwrite, inaccessiblemem: none) }
attributes #3 = { nofree nosync nounwind memory(argmem: read) }
attributes #4 = { mustprogress nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #5 = { nocallback nofree nosync nounwind speculatable willreturn memory(none) }
attributes #6 = { nocallback nofree nosync nounwind willreturn memory(argmem: readwrite) }

!llvm.module.flags = !{!0, !1, !2}

!0 = !{i32 2, !"Debug Info Version", i32 3}
!1 = !{i32 8, !"PIC Level", i32 2}
!2 = !{i32 7, !"PIE Level", i32 2}
!3 = !{!4, !4, i64 0}
!4 = !{!"any data access", !5, i64 0}
!5 = !{!"any access", !6, i64 0}
!6 = !{!"Flang Type TBAA Root"}
!7 = !{!8, !8, i64 0}
!8 = !{!"descriptor member", !5, i64 0}
