比赛过程
17:50
cccc来晚了来晚了gggggg….. 看了一遍题, woc今年的省选??? 感觉要gg. 怎么一道题都没有见过啊, 今天又要爆零啦…. 看了一遍t1, 好长啊? 所以看完了还要数据分治? 看第二题… 期望? 省选级别的期望都不简单啊, 跪了跪了跪了…. 看第三题? dp? 100数据范围感觉可以乱搞啊?
20:00
想了三道题无果, 开始t1数据分治, 感觉能拿84pts? 这也不差了吧… 但是实在写的太烦躁了决定冷静一下, 写完懒得编译, 先去做其他题(千万不要学我). 第三题看了样例还有负数? 再看了一下数据范围100… 这不最大权闭合子图吗? 这玩意没敲过啊, 但是大概建图我知道. 于是开始码码码码… 还是挺好写的? 感觉建边方式有点多啊, 有点虚? 调了一遍样例, 死了? 检查一下发现bfs没pop, 老毛病了? 还死? 检查了一下读优发现没判负数. 怎么还死? 突然发现没死原来是还可以输入… 检查了一下读入发现有个地方写错了, 改了一下所有样例瞬间全部过完了, 心中舒畅…dinic就是好写啊.
21:30
重新看t2, 感觉还是不会做, 也就没有看部分分的数据保证(事实证明这是今日最后悔的事情). 看t1去了. 继续修改了t1的暴力, 发现到了22:00没调出来? 十分钟狂飙手速一波改成低级暴力… 发现还是过不了? 心中慌.
22:10
t1弃疗, 估分100+, 考试结束.
赛后
t3 a了, 爽歪歪. tm t1真没分??? 我写了那么久的数据分治啊… 太贪心了想去拿84pts. 早知道写个最朴素的暴力把他写正确也很不错啊. t2居然有n==k的时候???? 那不是倒着for一遍就有分? 瞬间敲完暴力上去发现有80pts??? woc怎么没有去看t2的部分数据啊我好蠢啊…
反思
1.无论多虚都要看一下部分数据保证, 说不定能骗很多分.
2. 写哪道题一定要想清楚了写, 写就一定要写的能得分, 要不然等于白费时间, 即使是暴力也要想清楚了再写.
题解
分手是祝愿
#include<bits/stdc++.h>
using namespace std;
typedef long long lnt;
const int mod = 1e5 + 3;
const int maxn = 1e5 + 5;
int n, k, num, a[maxn];
lnt inv[maxn], b[maxn], ans;
int main() {
inv[1] = 1;
scanf("%d%d", &n, &k);
for (int i = 1; i <= n; ++ i) scanf("%d", &a[i]);
for (int i = 2; i <= n; ++ i)
inv[i] = (-(mod / i) * inv[mod % i] % mod + mod) % mod;
for (int i = n; i; -- i)
if (a[i]) {
++ num;
for (int j = 1; j * j <= i; ++ j)
if (i % j == 0) {
a[j] ^= 1;
if (j * j != i) a[i / j] ^= 1;
}
}
for (int i = n; i; -- i) b[i] = ((n - i) * b[i + 1] % mod + n) * inv[i] % mod;
if (k >= num) ans = num;
else {
b[k] = k;
for (int i = k; i <= num; ++ i)
ans = (ans + b[i]) % mod;
}
for (int i = 2; i <= n; ++ i)
ans = ans * i % mod;
printf("%lld\n", ans);
return 0;
}
寿司
#include<bits/stdc++.h>
#define deeper(a) memset(a, -1, sizeof(a))
using namespace std;
const int inf = 1e9;
const int maxm = 2e6;
const int maxn = 4e4;
queue<int> q;
long long ans;
bool vis[1005];
int n, num, S, T, cnt, m;
int h[maxn], dis[maxn], a[maxn], id[105][105], d[105][105];
struct edge {
int nxt, v, c;
}e[maxm];
inline void add(int u, int v, int w) {
e[++ num].v = v, e[num].nxt = h[u], e[num].c = w, h[u] = num;
e[++ num].v = u, e[num].nxt = h[v], e[num].c = 0, h[v] = num;
}
inline const int read() {
register int x = 0, f = 1;
register char ch = getchar();
while (ch < '0' || ch > '9') {if (ch == '-') f = -1; ch = getchar();}
while (ch >= '0' && ch <= '9') x = (x << 3) + (x << 1) + ch - '0', ch = getchar();
return f * x;
}
inline bool bfs() {
deeper(dis);
dis[S] = 0, q.push(S);
while (!q.empty()) {
int u = q.front(); q.pop();
for (int i = h[u]; i; i = e[i].nxt) {
int v = e[i].v;
if (dis[v] == -1 && e[i].c)
dis[v] = dis[u] + 1, q.push(v);
}
}
return dis[T] != -1;
}
inline int dfs(int u, int a) {
if (u == T) return a;
int f, flow = 0;
for (int i = h[u]; i; i = e[i].nxt) {
int v = e[i].v;
if (dis[v] == dis[u] + 1 && e[i].c) {
f = dfs(v, min(e[i].c, a));
e[i].c -= f, e[i ^ 1].c += f;
a -= f, flow += f;
if (!a) return flow;
}
}
if (!flow) dis[u] = -1;
return flow;
}
int main() {
freopen("sushi.in", "r", stdin);
freopen("sushi.out", "w", stdout);
n = read(), m = read();
for (int i = 1; i <= n; ++ i)
a[i] = read(), vis[a[i]] = true;
T = cnt = 1001, num = 1;
for (int i = 1; i <= 1000; ++ i)
if (vis[i]) add(i, T, m * i * i);
for (int i = 1; i <= n; ++ i)
for (int j = i; j <= n; ++ j) d[i][j] = read(), id[i][j] = ++ cnt;
for (int i = 1; i <= n; ++ i)
for (int j = 1; j <= n; ++ j) {
if (i == j)
d[i][j] -= a[i], add(id[i][j], a[i], inf);
if (d[i][j] > 0) {
ans += d[i][j];
add(S, id[i][j], d[i][j]);
}
if (d[i][j] < 0) add(id[i][j], T, -d[i][j]);
if (i != j) {
add(id[i][j], id[i + 1][j], inf);
add(id[i][j], id[i][j - 1], inf);
}
}
while (bfs()) ans -= dfs(S, inf);
printf("%lld\n", ans);
return 0;
}