結果

問題 No.896 友達以上恋人未満
コンテスト
ユーザー chocorusk
提出日時 2019-09-28 10:50:43
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 1,451 bytes
コンパイル時間 1,004 ms
コンパイル使用メモリ 100,084 KB
実行使用メモリ 134,400 KB
最終ジャッジ日時 2024-10-01 19:25:04
合計ジャッジ時間 6,320 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 1 WA * 6
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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];
int main()
{
	cin>>m>>n>>mulx>>addx>>muly>>addy>>mod;
	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=1; i<mod; i++){
		z[i]+=z[0];
		for(int j=2*i; j<mod; j+=i){
			z[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;
}
0