結果
問題 | No.896 友達以上恋人未満 |
ユーザー |
![]() |
提出日時 | 2019-09-27 22:01:20 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2,058 ms / 3,500 ms |
コード長 | 1,413 bytes |
コンパイル時間 | 2,293 ms |
コンパイル使用メモリ | 195,480 KB |
最終ジャッジ日時 | 2025-01-07 19:24:26 |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 7 |
ソースコード
#include<bits/stdc++.h>using namespace std;using Int = long long;template<typename T1,typename T2> inline void chmin(T1 &a,T2 b){if(a>b) a=b;}template<typename T1,typename T2> inline void chmax(T1 &a,T2 b){if(a<b) a=b;}struct FastIO{FastIO(){cin.tie(0);ios::sync_with_stdio(0);}}fastio_beet;//INSERT ABOVE HEREusing ll = long long;const int sz = 1<<24;ll zs[sz]={};signed main(){ll m,n,mulX,addX,mulY,addY,MOD;cin>>m>>n>>mulX>>addX>>mulY>>addY>>MOD;vector<ll> xs(m),ys(m),as(m),bs(m);for(int i=0;i<m;i++) cin>>xs[i];for(int i=0;i<m;i++) cin>>ys[i];for(int i=0;i<m;i++) cin>>as[i];for(int i=0;i<m;i++) cin>>bs[i];for(int i=0;i<m;i++){zs[xs[i]]+=ys[i];}ll px,py,pa,pb;px=xs[m-1],py=ys[m-1],pa=as[m-1],pb=bs[m-1];for(int i=m;i<n;i++){ll nx=(px*mulX+addX)&(MOD-1);ll ny=(py*mulY+addY)&(MOD-1);zs[nx]+=ny;px=nx;py=ny;}for(int i=1;i<sz;i++)for(int j=i+i;j<sz;j+=i)zs[i]+=zs[j];ll sum=0;for(int i=0;i<m;i++){ll res=zs[as[i]];if(as[i]*bs[i]<MOD) res-=zs[as[i]*bs[i]];cout<<res<<"\n";sum^=res;}px=xs[m-1],py=ys[m-1],pa=as[m-1],pb=bs[m-1];for(int i=m;i<n;i++){ll na=((pa*mulX+addX+MOD-1)&(MOD-1))+1;ll nb=((pb*mulY+addY+MOD-1)&(MOD-1))+1;ll res=zs[na];if(na*nb<MOD) res-=zs[na*nb];sum^=res;pa=na;pb=nb;}cout<<sum<<endl;return 0;}