結果

問題 No.896 友達以上恋人未満
ユーザー milanis48663220
提出日時 2019-09-28 00:57:16
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 1,403 ms / 3,500 ms
コード長 1,357 bytes
コンパイル時間 591 ms
コンパイル使用メモリ 56,740 KB
実行使用メモリ 173,668 KB
最終ジャッジ日時 2024-09-25 04:00:02
合計ジャッジ時間 4,553 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 7
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>

using namespace std;

int z[(1<<24)+1];
long cnt[(1<<24)+1];
long x[1000], y[1000], a[1000], b[1000];

int main(){
    long M, N, mx, my, ax, ay, MOD;
    cin >> M >> N >> mx >> ax >> my >> ay >> MOD;
    for(long i = 0; i < M; i++){
        cin >> x[i];
    }
    for(long i = 0; i < M; i++){
        cin >> y[i];
    }
    for(long i = 0; i < M; i++){
        cin >> a[i];
    }
    for(long i = 0; i < M; i++){
        cin >> b[i];
        z[x[i]] += y[i];
    }
    long xp = x[M-1];
    long yp = y[M-1];
    for(long i = M; i < N; i++){
        long X, Y;
        X = xp*mx+ax;
        Y = yp*my+ay;
        X &= (MOD-1);
        Y &= (MOD-1);
        z[X] += Y;
        xp = X;
        yp = Y;
    }
    for(long i = 1; i < MOD; i++){
        for(long j = 1; i*j < MOD; j++){
            cnt[i] += z[i*j];
        }
    }
    long xo = 0;
    for(long i = 0; i < M; i++){
        long ans = cnt[a[i]];
        if(a[i]*b[i] < MOD) ans -= cnt[a[i]*b[i]];
        cout << ans << endl;
        xo ^= ans;
    }
    long ap = a[M-1], bp = b[M-1];
    for(long i = M; i < N; i++){
        long A, B;
        A = ((ap*mx+ax+MOD-1)&(MOD-1))+1;
        B = ((bp*my+ay+MOD-1)&(MOD-1))+1;
        long ans = cnt[A];
        if(A*B < MOD) ans -= cnt[A*B];
        xo ^= ans;
        ap = A;
        bp = B;
    }
    cout << xo << endl;
}
0