R=range
J=lambda:map(int,input().split())
O=print
N,*_=J()
A=list(J())
A=sorted([A[i],i+1]for i in R(N*2))
O(sum(A[N+i][0]-A[i][0]for i in R(N)))
for i in R(N):O(A[i][1],A[N+i][1])