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;
}

转载于:https://www.cnblogs.com/ldxcaicai/p/9738339.html

2018.08.29 NOIP模拟 pmatrix（线性筛）相关推荐

2018.11.03 NOIP模拟地球发动机（线性dp）
传送门考试5分钟写完. 如果没这题今天多半爆零了(汗直接二分出合法的转移范围. 然后用后面的状态更新前面的就可以了. 代码
2018.10.20 NOIP模拟蛋糕（线段树+贪心/lis）
传送门听说是最长反链衍生出的对偶定理就能秒了. 本蒟蒻直接用线段树模拟维护的. 对于第一维排序. 维护第二维的偏序关系可以借助线段树/树状数组维护逆序对的思想建立权值线段树贪心求解. 代码转载于: ...
2018.10.26 NOIP模拟瓶子（dp/贪心）
传送门正解是dp并不想去想了. 自己yy了一个贪心拿了95pts95pts95pts,唯一没过的点还只有一个地方错了,面向数据变成之后过啦! 所以我讲讲如何贪心. 考虑到最后都只会合并成一种颜色,所 ...
EZ 2018 07 06 NOIP模拟赛
又是慈溪那边给的题目,这次终于没有像上次那样尴尬了, T1拿到了较高的暴力分,T2没写炸,然后T3写了一个优雅的暴力就203pts,Rank3了. 听说其它学校的分数普遍100+,那我们学校还不是强到 ...
2018/3/18 noip模拟赛 20分
T1 dp,特别裸特别简单,我放弃了写了个dfs. T2 树归,特别裸特别简单,我不会写. T3 贪心二分不知道什么玩意儿反正不会写就对了. 我是个智障转载于:https://www.cnblogs ...
2018.11.05 NOIP模拟规避（最短路计数）
传送门正难则反. 考虑计算两人相遇的方案数. 先正反跑一遍最短路计数. 然后对于一条在最短路上的边(u,v)(u,v)(u,v),如果(dis(s,u)*2<total&&di ...
2018.10.24 NOIP模拟小 C 的序列（链表+数论）
传送门考虑到a[l],gcd(a[l],a[l+1]),gcd(a[l],a[l+1],a[l+2])....gcd(a[l]...a[r])a[l],gcd(a[l],a[l+1]),gcd(a[ ...
2018.10.01 NOIP模拟卡牌游戏（贪心）
传送门简单贪心题. 然而考试的时候失了智少讨论了一种情况导致gg. 实际上用到了二分图匹配的思想,L每次找到刚好比当前的牌小一点的出出去,看能匹配几个. 如何处理? 我们先考虑第一种比分策略. 我们 ...
2018.10.04 NOIP模拟航班（tarjan+树形dp）
传送门考场上自己yy了一个双连通只有40分. 然后换根dp求最长路就行了. 代码转载于:https://www.cnblogs.com/ldxcaicai/p/10084917.html