区间【NOIP2017提高A组模拟7.10】
2017-07-15 11:00
323 查看
题目
输入
样例输入
sample1:4 2 10
5 1 1 10
sample2:
1000 97 96998351
41 1668 505 2333
样例输出
sample1:
4
sample2:
1749769
数据范围
思路
当时想着,一个新区间不就是加上后面一个与减去第一个嘛,求出逆元后,扫一遍就可以了。←_←然而模P不一定有逆元。
这就很尴尬了。。。。。
所以就应该分块!
解法
分块!每k个分一块(最后一块不一定有k个)。
对于一个区间,它要么是一个块,要么是横跨两个块。所以只需求出每一个块的前缀与后缀,扫一遍计算即可。
时间O(n)。
代码
#include<cstdio> #include<algorithm> #include<cstdlib> #include<cstring> #include<cmath> #define fo(i,a,b) for(ll i=a;i<=b;i++) #define down(i,a,b) for(ll i=a;i>=b;i--) #define ll long long using namespace std; const int maxn=1e7*2+5; ll ans,n,k,mo,A,B,C,D,block; int a[maxn],sl[maxn],sr[maxn]; int main() { //freopen("T.in","r",stdin); scanf("%lld%lld%lld",&n,&k,&mo); scanf("%lld%lld%lld%lld",&A,&B,&C,&D); a[1]=A; fo(i,2,n) a[i]=((ll)a[i-1]*B+C)%D; block=n/k; fo(i,1,block*k) { if ((i-1)%k==0) sl[i]=(ll)a[i]%mo; else sl[i]=(ll)sl[i-1]*(ll)a[i]%mo; } down(i,block*k,1){ if (i%k==0) sr[i]=(ll)a[i]%mo; else sr[i]=(ll)sr[i+1]*(ll)a[i]%mo; } if (n%k!=0){ fo(i,block*k+1,n) { sl[i]=sr[i]=a[i]%mo; if (i!=block*k+1) sl[i]=(ll)sl[i-1]*(ll)a[i]%mo; } } fo(i,1,block*k){ ll now; if (i+k-1>n) break; if ((i-1)%k==0) now=(ll)sr[i]%mo; else now=(ll)sr[i]*(ll)sl[i+k-1]%mo; ans=ans^now; } printf("%lld\n",ans); }
相关文章推荐
- JZOJ 100035【NOIP2017提高A组模拟7.10】区间
- JZOJ 100035. 【NOIP2017提高A组模拟7.10】区间
- 区间【NOIP2017提高A组模拟9.16】
- 【NOIP2017提高组模拟6.27】C
- 【NOIP2017提高组模拟12.24】C
- 【NOIP2017提高A组模拟10.5】Ping
- 【NOIP2017提高组模拟12.10】幻魔皇
- A【NOIP2017提高组模拟12.18】
- JZOJ5347【NOIP2017提高A组模拟9.5】遥远的金字塔 斜率优化 DP
- NOIP2017提高组模拟22-雅礼国庆10.5
- 【NOIP2017提高A组模拟8.24】早苗
- 【NOIP2017提高A组模拟8.25】夜莺与玫瑰
- B【NOIP2017提高组模拟12.18】
- 【NOIP2017提高组模拟12.17】环
- 【NOIP2017提高组模拟12.24】B
- 【NOIP2017提高A组模拟9.12】Shorten Diameter
- 【NOIP2017提高A组模拟9.23】碎
- 【NOIP2017提高A组模拟10.10】Graph
- 【NOIP2017提高A组模拟8.22】密码
- 【NOIP2017提高组模拟12.10】幻魔皇