結果
問題 | No.896 友達以上恋人未満 |
ユーザー | chocorusk |
提出日時 | 2019-09-28 11:10:54 |
言語 | C++11 (gcc 11.4.0) |
結果 |
AC
|
実行時間 | 807 ms / 3,500 ms |
コード長 | 1,714 bytes |
コンパイル時間 | 900 ms |
コンパイル使用メモリ | 102,128 KB |
実行使用メモリ | 99,092 KB |
最終ジャッジ日時 | 2024-10-01 19:38:22 |
合計ジャッジ時間 | 5,785 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 172 ms
20,888 KB |
testcase_01 | AC | 155 ms
21,712 KB |
testcase_02 | AC | 174 ms
20,216 KB |
testcase_03 | AC | 162 ms
20,120 KB |
testcase_04 | AC | 167 ms
20,116 KB |
testcase_05 | AC | 399 ms
52,840 KB |
testcase_06 | AC | 565 ms
52,548 KB |
testcase_07 | AC | 373 ms
52,832 KB |
testcase_08 | AC | 373 ms
53,728 KB |
testcase_09 | AC | 575 ms
85,880 KB |
testcase_10 | AC | 807 ms
99,092 KB |
ソースコード
#include <cstdio> #include <cstring> #include <iostream> #include <string> #include <cmath> #include <bitset> #include <vector> #include <map> #include <set> #include <queue> #include <deque> #include <algorithm> #include <complex> #include <unordered_map> #include <unordered_set> #include <random> #include <cassert> #include <fstream> #include <utility> #include <functional> #define popcount __builtin_popcount using namespace std; typedef long long int ll; typedef pair<int, int> P; int m, n, mulx, addx, muly, addy, mod; ll z[1<<24]; int a[1010], b[1010], x[1010], y[1010]; const int MAX=1<<24; bool isprime[MAX]; void sieve(){ for(ll i=3; i<MAX; i+=2){ isprime[i]=1; } isprime[2]=1; for(ll i=3; i<MAX; i++){ if(isprime[i]){ for(ll j=2*i; j<MAX; j+=i){ isprime[j]=0; } } } return; } int main() { cin>>m>>n>>mulx>>addx>>muly>>addy>>mod; sieve(); for(int i=1; i<=m; i++) cin>>x[i]; for(int i=1; i<=m; i++) cin>>y[i]; for(int i=1; i<=m; i++) cin>>a[i]; for(int i=1; i<=m; i++) cin>>b[i]; for(int i=1; i<=m; i++){ z[x[i]]+=y[i]; } ll x1=x[m], y1=y[m]; for(int i=m+1; i<=n; i++){ x1=((x1*mulx+addx)&(mod-1)); y1=((y1*muly+addy)&(mod-1)); z[x1]+=y1; } for(int i=2; i<mod; i++){ if(!isprime[i]) continue; for(int j=(mod-1)/i*i; j>0; j-=i){ z[j/i]+=z[j]; } } ll ans=0; for(int i=1; i<=m; i++){ ll ans1=0; if(a[i]<mod) ans1+=z[a[i]]; if((ll)a[i]*b[i]<mod) ans1-=z[a[i]*b[i]]; cout<<ans1<<endl; ans^=ans1; } ll a1=a[m], b1=b[m]; for(int i=m+1; i<=n; i++){ a1=((a1*mulx+addx+mod-1)&(mod-1))+1; b1=((b1*muly+addy+mod-1)&(mod-1))+1; ll ans1=0; if(a1<mod) ans1+=z[a1]; if(a1*b1<mod) ans1-=z[a1*b1]; ans^=ans1; } cout<<ans<<endl; return 0; }