Problem: problem

You are playing with J using the Midakuji. The Amida lottery consists of n vertical bars and m horizontal bars. The vertical bars are numbered from 1 to n in order from the left, and the vertical bar i The positive integer si is written at the bottom of.

<image>

Figure 3-1 Example of Amidakuji (n = 4, m = 5, s1 = 20, s2 = 80, s3 = 100, s4 = 50)

The integer written at the bottom of the vertical bar i, which is reached by following the path from the top of the vertical bar i, is the score when the vertical bar i is selected. For example, in Fig. 3-1 the vertical bar 1 is selected. And the score is 80 points, and if the vertical bar 2 is selected, the score is 100 points.

<image>

Figure 3-2 Example of how to follow the road

Mr. J decides to choose a series of k books from vertical bar 1 to vertical bar k. The total score when selecting those k vertical bars is J's score. However, you are in the Amidakuji You can select one horizontal bar and delete that horizontal bar from the Amidakuji. (You do not have to delete it.) If you delete one horizontal bar, the vertical bar in the deleted Amidakuji will be displayed vertically. J's score is the total of the points when selecting k consecutive vertical bars from bar 1 to vertical bar k.

Given the shape of the Amida lottery and the number of vertical bars k that J chooses as input, create a program that finds the minimum value of J's score.



input

The input consists of multiple datasets. Each dataset is given in the following format.

On the first line, four integers n, m, h, k are written separated by blanks. N (2 ≤ n ≤ 1000) is the number of vertical bars, and m (1 ≤ m ≤ 100000) is the horizontal bar. H (2 ≤ h ≤ 1000) represents the length of the vertical bar, and k (1 ≤ k ≤ n) represents the number of vertical bars you choose.

The following n lines contain the score written at the bottom of the vertical bar. The first line (1 ≤ i ≤ n) contains the positive integer si. Also, s1 + s2 + ... + sn ≤ 2000000000 = 2 × 109.

The following m lines indicate the position of the horizontal bar. The horizontal bars are numbered from 1 to m. The first line (1 ≤ i ≤ m) is the horizontal bar i. The two integers ai, bi (1 ≤ ai ≤ n --1, 1 ≤ bi ≤ h --1) representing the position are separated by blanks, and the horizontal bar i connects the vertical bar ai and the vertical bar ai + 1. , Indicates that the distance from the top edge of the horizontal bar i is bi. However, no two horizontal bars share an endpoint.

Of the scoring data, 20% of the points will have the lowest score for Mr. J if the horizontal bar is not deleted. Also, 30% of the points will satisfy n ≤ 20, m ≤ 30, h ≤ 10. 60% of the points are m ≤ 1000.

The end of the input is indicated by a line containing four zeros. The number of datasets does not exceed 10.

output

For each dataset, output the minimum score of Mr. J on one line.

Examples

Input

4 5 7 2
20
80
100
50
1 1
2 6
2 3
1 5
3 1
2 2 5 1
10
20
1 1
1 3
0 0 0 0


Output

100
10


Input

None


Output

None

Solution (in CPP):

Solution:
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <climits>
#include <iostream>
#include <vector>
#include <string>
#include <sstream>
#include <algorithm>
#include <utility>
#include <set>
#include <map>
#include <stack>
#include <queue>
#include <deque>
#include <functional>

using namespace std;

#define int long long
#define MP make_pair
#define PB push_back
#define FI first
#define SE second
typedef pair<int, pair<int, int> > PPII;
typedef pair<int, int> PII;
static const int INF = 1ll<<60;

int N,M,H,K;
vector<int> V[1005];
int a[1005];
int flg[1005][1005];
vector<int> pass[1005][1005];
int vi[1005];
int s;
int maxd;
int temp=0;
int ans=0;

void dfs(int n,int m,bool down){
    //cout<<n<<" "<<V[n][m]<<" "<<s<<endl;
    pass[n][V[n][m]].PB(s);
    if(V[n][m]==maxd+1){
        vi[s]=a[n];
        return ;
    }
    if(down==true&&m!=0){
        int it = lower_bound(V[n+1].begin(),V[n+1].end(),V[n][m])-V[n+1].begin();
        if(it==V[n+1].size()||flg[n][V[n][m]]==0||V[n+1][it]!=V[n][m]){
            it = lower_bound(V[n-1].begin(),V[n-1].end(),V[n][m])-V[n-1].begin();
            dfs(n-1,it,false);
        }
        else dfs(n+1,it,false);
    }
    else {
        int it = upper_bound(V[n].begin(),V[n].end(),V[n][m])-V[n].begin();
        dfs(n,it,true);
    }
}

signed main(){
    while(1){
        cin>>N>>M>>H>>K;
        if(N==0)break;
        for(int i=0;i<N;++i){
            cin>>a[i];
        }
        for(int i=0;i<M;++i){
            int A,B;
            cin>>A>>B;
            A--;
            maxd=max(maxd,B);
            V[A].PB(B);
            V[A+1].PB(B);
            flg[A][B]=1;
        }
        for(int i=0;i<N;++i){
            V[i].PB(0);
            V[i].PB(maxd+1);
        }
        for(int i=0;i<N;++i)sort(V[i].begin(),V[i].end());
        for(int i=0;i<N;++i){
            s=i;
            int it = lower_bound(V[i].begin(),V[i].end(),0)-V[i].begin();
            dfs(i,it,true);
        }
        for(int i=0;i<K;++i)temp+=vi[i];
        ans=temp;
        for(int i=0;i<N;++i){
            for(int j=0;j<V[i].size();++j){
                int it = lower_bound(V[i+1].begin(),V[i+1].end(),V[i][j])-V[i+1].begin();
                if(it>=V[i+1].size())continue;
                if(V[i+1][it]!=V[i][j])continue;
                //cout<<V[i+1][it]<<" "<<V[i][j]<<endl;
                if(pass[i][V[i][j]].size()!=2||pass[i+1][V[i][j]].size()!=2)continue;
                if(min(pass[i][V[i][j]][0],pass[i][V[i][j]][1])!=min(pass[i+1][V[i][j]][0],pass[i+1][V[i][j]][1]))continue;
                if(max(pass[i][V[i][j]][0],pass[i][V[i][j]][1])!=max(pass[i+1][V[i][j]][0],pass[i+1][V[i][j]][1]))continue;
                int p1=min(pass[i][V[i][j]][0],pass[i][V[i][j]][1]),p2=max(pass[i+1][V[i][j]][0],pass[i+1][V[i][j]][1]);
                //cout<<p1<<" "<<p2<<endl;
                if(p1>=K&&p2>=K)continue;
                if(p1<K&&p2<K)continue;
                //cout<<p1<<" "<<p2<<" "<<temp-vi[min(p1,p2)]+vi[max(p1,p2)]<<endl;
                if(p1>=K||p2>=K)ans=min(ans,temp-vi[min(p1,p2)]+vi[max(p1,p2)]);
            }
        }
        cout<<ans<<endl;
        ans=0;
        temp=0;
        maxd=0;
        for(int i=0;i<N;++i){
            for(int j=0;j<V[i].size();++j){
                flg[i][V[i][j]]=0;
                pass[i][V[i][j]].clear();
            }
            V[i].clear();
            vi[i]=0;
        }
    }
}