#include<bits/stdc++.h>
using namespace std;
int a[500005];
long long S[500005];
long long f[500005];
int R[500005];
int sm[15000005],ls[15000005],rs[15000005],tot=0;
void fuzhi(int l,int r,int o,int ll,int rr,int z){
	if(l>=ll&&r<=rr){
		sm[o]=min(sm[o],z);
		return;
	}
	int mid=(l+r)>>1;
	if(mid>=ll){
		if(!ls[o])ls[o]=++tot,sm[tot]=1e9;
		fuzhi(l,mid,ls[o],ll,rr,z);
	}
	if(mid<rr){
		if(!rs[o])rs[o]=++tot,sm[tot]=1e9;
		fuzhi(mid+1,r,rs[o],ll,rr,z);
	}
}
int qiu(int l,int r,int o,int x){
	if(!o)return 1e9;
	if(l==r)return sm[o];
	int mid=(l+r)>>1,ans=sm[o];
	if(x<=mid)ans=min(ans,qiu(l,mid,ls[o],x));
	else ans=min(ans,qiu(mid+1,r,rs[o],x));
	return ans;
}
int n,m,k,c;
void upd(long long &x,long long y){x=min(x,y);}
int main(){
	int T;
	scanf("%d",&T);
	while(T--){
		scanf("%d%d%d%d",&n,&m,&k,&c);
		for(int i=1;i<=n;++i)scanf("%d",&a[i]),S[i]=S[i-1]+a[i];
		if(m==1){
			long long ans=0;
			for(int i=1;i<=n;++i){
				ans+=1ll*a[i]/k*c;
				ans+=min(c,a[i]%k);
			}
			printf("%lld\n",ans);
			continue;
		}
		for(int i=0;i<=n;++i)f[i]=1e18;
		f[n]=0;
		for(int i=1;i<=tot;++i)sm[i]=ls[i]=rs[i]=0;
		tot=1;sm[1]=n;
		for(int i=n;i>=1;--i){
			if(i+m<=n){
				int j=i+m;
				long long z=k+S[j-m]-S[j-1]-1;
				if(z>=0){
					int L=(S[j-m]-z)%k,R=S[j-m]%k;
					if(L<0)L+=k;
					if(R<0)R+=k;
					if(L<=R){
						fuzhi(0,k-1,1,L,R,j-1);
					}else{
						fuzhi(0,k-1,1,L,k-1,j-1);
						fuzhi(0,k-1,1,0,R,j-1);
					}
				}
			}
			R[i-1]=qiu(0,k-1,1,S[i-1]%k);
		}
		f[0] = 0;
		for(int i=1;i<=n;i++)f[i]=1e18;
		// for(int i=0;i<n;i++)std::cerr << R[i] << std::endl;
		for(int i=0;i<n;i++){
			int j=R[i];long long s;
			// std::cerr << i << " " << j << std::endl;
			s=S[j]-S[i];
			// upd(f[j],f[i]+s/k*c+s%k);
			upd(f[j],f[i]+(s+k-1)/k*c);
			s=S[j-m+1]-S[i];
			upd(f[j-m+2],f[i]+s/k*c+s%k);
		}
		// std::cerr << f[5] << std::endl;
		cout<<f[n]<<endl;
	}
	return 0;
}