題意：給我們n個數，然後給我們一個數字K，再給我們兩個區間[L,R],[U,V]，問我們從這兩個區間裡分別取一個數，有多少種取法，可以使得所取的兩個數和為K。

分析：莫隊演算法+容斥定理。因為莫隊演算法是知道[l,r]的值後可以求[l’,r’]的值，因為每次都會給兩個區間，所以我們需要將兩個不相關的區間想辦法讓他們產生連續。對於[x1,y1][x2,y2]我們這時候可以通過容斥定理將其變為F[x1,y1,x2,y2]=F[x1,y2]+F[y1+1,x2-1]-F[y1+1,y2]-F[x1,x2-1]將一次求解變成4個區間分別求解，每次在一個區間中求解有多少種取法可以使得兩個數和為k。

AC程式碼如下：

#include <iostream>
#include <stdio.h>
#include <math.h>
#include <string.h>
#include <algorithm>
#include <iostream>
#include <stdlib.h>
#include <queue>
#include <map>
#include <vector>

using namespace std;

const int maxn = 30000+10;

int N,K,a[maxn],m,belong[maxn],ans;
int
 
 sum[2*maxn];
struct Query
{
    int x1,y1,x2,y2,ans;
}q[maxn];///儲存查詢

struct Query1
{
    int l,r,ans,flag,fa;

}arr[maxn*4];///儲存區間

bool cmp1(Query1 a, Query1 b)///分塊排序
{
    if(belong[a.l]==belong[b.l])
        return a.r<b.r;
    return belong[a.l]<belong[b.l];
}

void updata(int pos, int flag)
{
    sum[a[pos]]+=flag;
    if
 
(K>a[pos])
        ans+=sum[K-a[pos]]*flag;
}

int main()
{
    while(scanf("%d",&N)!=EOF)
    {
        scanf("%d",&K);
        for(int i=1; i<=N; i++)
            scanf("%d",&a[i]);

        int block =sqrt(N);
        for(int i=1; i<=N; i++)
            belong[i]=(i-1)/block+1;

        scanf("%d",&m);
        for(int i=1; i<=m; i++)
        {
            scanf("%d %d %d %d",&q[i].x1,&q[i].y1,&q[i].x2,&q[i].y2);
            q[i].ans=0;
            for(int j=0; j<4; j++)
            {
                arr[i+j*m].fa=i;
                arr[i+j*m].ans=0;
            }///初始化各區間

            arr[i].flag=1,arr[i].l=q[i].x1,arr[i].r=q[i].y2;///[x1,y2];

            arr[i+m].flag=1,arr[i+m].l=q[i].y1+1,arr[i+m].r=q[i].x2-1;///[y1+1,x2-1]

            arr[i+m*2].flag=0,arr[i+m*2].l=q[i].y1+1,arr[i+m*2].r=q[i].y2;///[y1+1,y2]

            arr[i+m*3].flag=0,arr[i+m*3].l=q[i].x1,arr[i+m*3].r=q[i].x2-1;///[x1,x2-1]
        }
        sort(arr+1,arr+4*m+1,cmp1);

        ans=0;
        int l=1,r=0;
       memset(sum,0,sizeof(sum));
        sum[a[1]]=1;

        for(int i=1; i<=4*m; i++)
        {
            while(l<arr[i].l)
            {
                updata(l,-1);
                l++;
            }

            while(l>arr[i].l)
            {
                l--;
                updata(l,1);
            }

            while(r<arr[i].r)
            {
                r++;
                updata(r,1);
            }

            while(r>arr[i].r)
            {
                updata(r,-1);
                r--;
            }
            arr[i].ans=ans;
        }
        for(int i=1; i<=4*m; i++)
        {
            int x=arr[i].fa;
            int y=arr[i].ans;

            if(arr[i].flag) q[x].ans+=y;
            else
                q[x].ans-=y;
        }
        for(int i=1; i<=m; i++)
            printf("%d\n",q[i].ans);
    }
    return 0;
}