2018.08.29 NOIP模擬 pmatrix（線性篩）

【問題描述】

根據哥德巴赫猜想（每個不小于 6 的偶數都可以表示為兩個奇素數之和），定義
哥德巴赫矩陣 A 如下：對于正整數對(i,j)，若 i+j 為偶數且 i,j 均為奇素數，則 Ai,j = 1，
否則 Ai,j = 0。現在有若干詢問(x1,y1,x2,y2)，你需要回答下列式子的值

【輸入】

第一行一個整數 m
接下來 m 行，每行四個整數 x1 y1 x2 y2，表示一個詢問

【輸出】

m 行，每行一個整數，表示對應詢問的答案

【輸入樣例】

1
1 1 3 5

【輸出樣例】

2
30%的數據保證 x2, y2, m ≤ 100
100%的數據保證 1 ≤ x1 ≤ x2 ≤ 10^6; 1 ≤ y1 ≤ y2 ≤ 10^6; m ≤ 1000

---

簽到題。
一開始題讀錯了。。。
差點忘記開long long。。。
就是線性篩之后維護一個素數出現次數的前綴和。
代碼：

#include<bits/stdc++.h>
#define N 1000005
#define ll long long
using namespace std;
inline int read(){int ans=0;char ch=getchar();while(!isdigit(ch))ch=getchar();while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch^48),ch=getchar();return ans;
}
int pri[N],tot=0;
ll cnt[N];
bool vis[N];
inline void init(int len){vis[1]=true;for(int i=2;i<=len;++i){if(!vis[i])pri[++tot]=i;for(int j=1;j<=tot;++j){if(pri[j]*i>len)break;vis[pri[j]*i]=1;if(i%pri[j]==0)break;}}for(int i=3;i<=len;++i)cnt[i]=vis[i]?cnt[i-1]:cnt[i-1]+1ll;
}
int main(){
//  freopen("pmatrix.in","r",stdin);
//  freopen("pmatrix.out","w",stdout);init(1000000);int m=read();while(m--){int a=read()-1,c=read()-1,b=read(),d=read();cout<<(cnt[b]-cnt[a])*(cnt[d]-cnt[c])<<'\n';}return 0;
}