Pagini recente » Cod sursa (job #2901560) | Cod sursa (job #145982) | Cod sursa (job #3234853) | Cod sursa (job #2701071) | Cod sursa (job #2872813)
#include <bits/stdc++.h>
#define INF 0x3f3f3f3f
using namespace std;
ifstream fin("robotei.in");
ofstream fout("robotei.out");
struct point {
int x, y;
};
const int kN = 1e3;
const int kM = 7e5;
int n, m, X, Y, modX, modY, offX, offY, len, cntX1[kN], cntY1[kN], cntX2[kN], cntY2[kN], dp[kN][kN];
vector<point> g[kN][kN];
point q[kN * kN];
int64_t sol[kM];
void findCycle(int x, int y) {
len += 1;
dp[x][y] = 1;
int xv = (x * x + offX) % modX;
int yv = (y * y + offY) % modY;
if (xv == X && yv == Y) {
return;
}
if (dp[xv][yv]) {
len = 0;
return;
}
findCycle(xv, yv);
}
void testCase() {
fin >> n >> m >> X >> Y >> modX >> modY >> offX >> offY;
if (X >= modX || Y >= modY) {
fout << "0 " << (int64_t)n * n - 1 << '\n' << "1 1\n";
return;
}
for (int i = 0; i < n; ++i) {
cntX1[((int64_t)i * i + offX) % modX] += 1;
cntY1[((int64_t)i * i + offY) % modY] += 1;
if (i >= modX) {
cntX2[((int64_t)i * i + offX) % modX] += 1;
}
if (i >= modY) {
cntY2[((int64_t)i * i + offY) % modY] += 1;
}
}
findCycle(X, Y);
for (int x = 0; x < modX; ++x) {
for (int y = 0; y < modY; ++y) {
int xv = (x * x + offX) % modX;
int yv = (y * y + offY) % modY;
g[xv][yv].push_back({x, y});
dp[x][y] = INF;
}
}
int l = 0, r = -1;
dp[X][Y] = 0;
q[++r] = {X, Y};
while (l <= r) {
point curr = q[l++];
for (const point &vec : g[curr.x][curr.y]) {
if (dp[vec.x][vec.y] > dp[curr.x][curr.y] + 1) {
dp[vec.x][vec.y] = dp[curr.x][curr.y] + 1;
q[++r] = vec;
}
}
}
/* for (int i = 0; i < modX; ++i) {
for (int j = 0; j < modY; ++j) {
for (int d = dp[i][j]; d <= dp[i][j] + 1; ++d) {
int visits = (d <= m), rem = max(0, m - d);
if (len) {
visits += rem / len;
}
if (d == dp[i][j]) {
sol[visits] += 1;
} else {
sol[visits] += (int64_t)cntX2[i] * cntY1[j] + (int64_t)cntX1[i] * cntY2[j] - (int64_t)cntX2[i] * cntY2[j];
}
}
}
} */
for (int i = 0; i <= m; ++i) {
if (sol[i]) {
fout << i << ' ' << sol[i] << '\n';
}
}
}
int main() {
int tests = 1;
for (int tc = 0; tc < tests; ++tc) {
testCase();
}
fin.close();
fout.close();
return 0;
}