結果

問題 No.896 友達以上恋人未満
コンテスト
ユーザー Lepton_s
提出日時 2018-03-02 18:21:46
言語 C++17
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=c++17 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
MLE  
(最新)
AC  
(最初)
実行時間 -
コード長 1,183 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 698 ms
コンパイル使用メモリ 81,308 KB
実行使用メモリ 213,108 KB
最終ジャッジ日時 2026-06-06 22:18:35
合計ジャッジ時間 6,005 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 5 MLE * 2
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

// naive2
#include <cstdio>
#include <cassert>
#include <iostream>
#include <cstring>
using namespace std;
typedef long long ll;
#define N (1<<24)
#define M 1000
int m, n, mx, ax, my, ay, mod;
ll x[M], y[M], z[N], a[M], b[M], dp[N];

ll solve2(ll a){
	if(a>mod) return 0;
	if(dp[a]>=0) return dp[a];
	ll res = 0;
	for(ll i = a; i < mod; i+=a){
		res += z[i];
	}
	return dp[a] = res;
}

int main(){
	scanf("%d%d%d%d%d%d%d",&m,&n,&mx,&ax,&my,&ay,&mod);
	for(int i = 0; i < m; i++) scanf("%d", x+i);
	for(int i = 0; i < m; i++) scanf("%d", y+i);
	for(int i = 0; i < m; i++) z[x[i]] += y[i];
	ll xx = x[m-1], yy = y[m-1];
	for(int i = m; i < n; i++){
		xx = (xx*mx+(ll)ax)&(mod-1);
		yy = (yy*my+(ll)ay)&(mod-1);
		z[xx] += yy;
	}
	
	for(int i = 0; i < m; i++) scanf("%d", a+i);
	for(int i = 0; i < m; i++) scanf("%d", b+i);

	memset(dp, -1, sizeof(dp));
	ll res = 0;
	for(int i = 0; i < m; i++){
		ll r = solve2(a[i])-solve2(a[i]*b[i]);
		printf("%lld\n", r);
		res^=r;
	}
	ll aa = a[m-1], bb = b[m-1];
	for(int i = m; i < n; i++){
		aa = ((aa*mx+(ll)ax+mod-1)&(mod-1))+1;
		bb = ((bb*my+(ll)ay+mod-1)&(mod-1))+1;
		res^=solve2(aa)-solve2(aa*bb);
	}
	printf("%lld\n", res);
	return 0;
}
0