// 以上帝的名义
// 候选硕士
#include <bits/stdc++.h>
#ifdef local
#include "algo/debug.h"
#else
#define dbg(x...) 0
#endif
using namespace std ;
using ll = long long ;
#define int ll
bool isPrime(int x) {
if (x < 2) return 0 ;
for (int i = 2 ; i * i <= x ; i++) {
if (x % i == 0) {
return 0 ;
}
}
return 1 ;
}
struct Fenwick {
vector<ll> t ;
int n;
Fenwick(int len) {
n = len ;
t.assign(n + 1, 0) ;
}
void upd(int i, ll x) {
while (i <= n) {
t[i] += x;
i += (i & (-i)) ;
}
}
ll get(int i) {
ll sum = 0 ;
while (i) {
sum += t[i] ;
i -= (i & (-i)) ;
}
return sum ;
}
ll get(int l, int r) {
return get(r) - get(l - 1) ;
}
};
int32_t main() {
cin.tie(0)->sync_with_stdio(false) ;
ll n, k ; cin >> n >> k ;
struct rect {
int x1, y1, x2, y2;
void read() {
cin >> x1 >> y1 >> x2 >> y2 ;
}
};
vector<rect> a(k) ;
for (int i = 0 ; i < k ; i++) a[i].read() ;
vector<int> d;
for (int i = 1 ; i * i <= n ; i++) {
if (n % i) continue ;
d.push_back(i) ;
if (n / i != i) d.push_back(n / i) ;
}
auto calc = [&](ll del) -> ll{
vector<vector<ll>> cnt(del + 1, vector<ll>(2, 0)) ;
int block = n / del ;
vector<Fenwick> f(2, Fenwick(n + 2)) ;
for (int i = 0 ; i < k ; i++) {
int x1 = a[i].x1 , y1 = a[i].y1 , x2 = a[i].x2, y2 = a[i].y2 ;
ll leny = (y2 - y1 + 1) ;
ll lenx = (x2 - x1 + 1) ;
ll evy = ((y1 & 1) ? (leny / 2) : (leny + 1) / 2) , ody = leny - evy ;
f[0].upd(x1, evy) ;
f[0].upd(x2 + 1, -evy) ;
f[1].upd(x1, ody) ;
f[1].upd(x2 + 1, -ody) ;
}
for (int i = 1 ; i <= n ; i++) {
// dbg(f[0].get(i), f[1].get(i)) ;
cnt[(i - 1) / block][1] += f[0].get(i) ;
cnt[(i - 1) / block][0] += f[1].get(i) ;
}
// dbg(del, cnt) ;
ll ret = 0 , sum = 0 ;
for (int i = 0 ; i < del ; i++) {
// for (int j = 0 ; j < del ; j++) {
if (i % 2 == 0) {
ll c = ((del + 1) / 2) * (block * block) ;
ll cur = c - cnt[i][0] + cnt[i][1] ;
sum += cur ;
} else {
ll c = (del / 2) * (block * block) ;
ll cur = c - cnt[i][1] + cnt[i][0] ;
sum += cur ;
}
// }
}
ret = sum ;
sum = 0 ;
for (int i = 0 ; i < del ; i++) {
// for (int j = 0 ; j < del ; j++) {
if (i % 2 != 0) {
ll c = ((del + 1) / 2) * (block * block) ;
ll cur = c - cnt[i][0] + cnt[i][1] ;
sum += cur ;
} else {
ll c = (del / 2) * (block * block) ;
ll cur = c - cnt[i][1] + cnt[i][0] ;
sum += cur ;
}
// }
}
ret = min(ret, sum) ;
return ret ;
};
ll ret = LLONG_MAX ;
for (int i : d) {
if (i > 1) {
ret = min(ret, calc(i)) ;
}
}
cout << ret ;
return 0 ;
}
// 希望白银
Compilation message
chessboard.cpp: In lambda function:
chessboard.cpp:76:16: warning: unused variable 'lenx' [-Wunused-variable]
76 | ll lenx = (x2 - x1 + 1) ;
| ^~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
348 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
3 |
Correct |
0 ms |
348 KB |
Output is correct |
4 |
Correct |
0 ms |
348 KB |
Output is correct |
5 |
Correct |
0 ms |
348 KB |
Output is correct |
6 |
Correct |
0 ms |
348 KB |
Output is correct |
7 |
Correct |
0 ms |
348 KB |
Output is correct |
8 |
Correct |
0 ms |
452 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
24 ms |
8280 KB |
Output is correct |
2 |
Correct |
9 ms |
2396 KB |
Output is correct |
3 |
Correct |
18 ms |
8796 KB |
Output is correct |
4 |
Correct |
13 ms |
2776 KB |
Output is correct |
5 |
Correct |
21 ms |
7552 KB |
Output is correct |
6 |
Correct |
17 ms |
5980 KB |
Output is correct |
7 |
Correct |
4 ms |
2600 KB |
Output is correct |
8 |
Correct |
19 ms |
5980 KB |
Output is correct |
9 |
Correct |
41 ms |
6720 KB |
Output is correct |
10 |
Correct |
18 ms |
5720 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
0 ms |
348 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
0 ms |
348 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
24 ms |
8280 KB |
Output is correct |
2 |
Correct |
9 ms |
2396 KB |
Output is correct |
3 |
Correct |
18 ms |
8796 KB |
Output is correct |
4 |
Correct |
13 ms |
2776 KB |
Output is correct |
5 |
Correct |
21 ms |
7552 KB |
Output is correct |
6 |
Correct |
17 ms |
5980 KB |
Output is correct |
7 |
Correct |
4 ms |
2600 KB |
Output is correct |
8 |
Correct |
19 ms |
5980 KB |
Output is correct |
9 |
Correct |
41 ms |
6720 KB |
Output is correct |
10 |
Correct |
18 ms |
5720 KB |
Output is correct |
11 |
Incorrect |
0 ms |
348 KB |
Output isn't correct |
12 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
1 ms |
348 KB |
Output is correct |
2 |
Correct |
0 ms |
348 KB |
Output is correct |
3 |
Correct |
0 ms |
348 KB |
Output is correct |
4 |
Correct |
0 ms |
348 KB |
Output is correct |
5 |
Correct |
0 ms |
348 KB |
Output is correct |
6 |
Correct |
0 ms |
348 KB |
Output is correct |
7 |
Correct |
0 ms |
348 KB |
Output is correct |
8 |
Correct |
0 ms |
452 KB |
Output is correct |
9 |
Correct |
24 ms |
8280 KB |
Output is correct |
10 |
Correct |
9 ms |
2396 KB |
Output is correct |
11 |
Correct |
18 ms |
8796 KB |
Output is correct |
12 |
Correct |
13 ms |
2776 KB |
Output is correct |
13 |
Correct |
21 ms |
7552 KB |
Output is correct |
14 |
Correct |
17 ms |
5980 KB |
Output is correct |
15 |
Correct |
4 ms |
2600 KB |
Output is correct |
16 |
Correct |
19 ms |
5980 KB |
Output is correct |
17 |
Correct |
41 ms |
6720 KB |
Output is correct |
18 |
Correct |
18 ms |
5720 KB |
Output is correct |
19 |
Incorrect |
0 ms |
348 KB |
Output isn't correct |
20 |
Halted |
0 ms |
0 KB |
- |