2024省选联测10
A. 小幸运
给出平面上 \(n\) 个点的坐标，以及整数 \(W，H\)。以每个点为底边中点构造底边长度相等且底边与一坐标轴平行的等腰直角三角形，满足三角形在 \((0,0),(W,0),(W,H),(0,H)\) 四点构成的矩形内部且三角形内部区域互不重叠。求每个三角形底边长度的最大值。
把所有坐标扩大两倍后答案只可能是整数。二分答案。每个点拆成 \(4\) 个小三角形，需要选择两个相邻的三角形。判断这 \(4n\) 个三角形是否相交，做 2-SAT。
如果两个三角形有边相交，那么这两个三角形相交。
点击查看代码
#include <cstdio>
#include <iostream>
#include <algorithm>
#define int long long
using namespace std;
const int N = 1000;
int n, w, h, tot, top, cnt, pos, x[N], y[N], s[N], sd[N], dfn[N], low[N], vis[N], head[N];
struct Edge {int to, nex;} e[1000005];
void add(int u, int v) {e[++tot] = {v, head[u]}; head[u] = tot;}
void tarjan(int u)
{
dfn[u] = low[u] = ++cnt; vis[s[++top] = u] = 1;
for(int i = head[u], v; i; i = e[i].nex)
{
if(!dfn[v = e[i].to]) tarjan(v), low[u] = min(low[u], low[v]);
else if(vis[v]) low[u] = min(low[u], low[v]);
}
if(dfn[u] != low[u]) return; ++ pos; while(top)
{
sd[s[top]] = pos; vis[s[top]] = 0;
if(s[top --] == u) break;
}
}
#define id(x, y) x + y * n
#define fid(x, y) x + y * n + 4 * n
bool t(int i, int j, int d)
{
if(j <= 1 && y[i] + d > h) return 0;
if(j >= 2 && y[i] - d < 0) return 0;
if(j % 2 == 0 && x[i] - d < 0) return 0;
if(j % 2 == 1 && x[i] + d > w) return 0;
return 1;
}
bool g(int x01, int y01, int x02, int y02, int x11, int y11, int x12, int y12)
{
if(x01 == x02 && x11 == x12) return 0;
if(x01 > x02) swap(x01, x02), swap(y01, y02);
if(x11 > x12) swap(x11, x12), swap(y11, y12);
if(x01 == x02 || x11 == x12)
{
if(x11 == x12) swap(x01, x11), swap(y01, y11), swap(x02, x12), swap(y02, y12);
double k = 1.0 * (y11 - y12) / (x11 - x12), b = y11 - k * x11, y = k * x01 + b;
if(y01 > y02) swap(y01, y02); if(y > y01 && y < y02 && x01 > x11 && x01 < x12) return 1; return 0;
}
double k1 = 1.0 * (y01 - y02) / (x01 - x02), b1 = y01 - k1 * x01;
double k2 = 1.0 * (y11 - y12) / (x11 - x12), b2 = y11 - k2 * x11;
if(k1 == k2) return 0; double x = - (b1 - b2) / (k1 - k2);
return x > x01 && x < x02 && x > x11 && x < x12;
}
int X[2][3], Y[2][3];
bool f(int p, int a, int q, int b, int d)
{
X[0][0] = x[p]; X[0][1] = x[p]; X[0][2] = x[p] + (a % 2 ? d : -d);
Y[0][0] = y[p]; Y[0][1] = y[p] + (a < 2 ? d : -d); Y[0][2] = y[p];
X[1][0] = x[q]; X[1][1] = x[q]; X[1][2] = x[q] + (b % 2 ? d : -d);
Y[1][0] = y[q]; Y[1][1] = y[q] + (b < 2 ? d : -d); Y[1][2] = y[q];
for(int i = 0; i <= 2; ++ i) for(int j = 0; j <= 2; ++ j)
if(g(X[0][i], Y[0][i], X[0][(i + 1) % 3], Y[0][(i + 1) % 3], X[1][j], Y[1][j], X[1][(j + 1) % 3], Y[1][(j + 1) % 3])) return 1;
return 0;
}
bool check(int d)
{
tot = top = cnt = pos = 0;
for(int i = 1; i <= 8 * n; ++ i) dfn[i] = low[i] = head[i] = 0;
for(int i = 1; i <= n; ++ i)
{
add(id(i, 0), fid(i, 3)); add(fid(i, 0), id(i, 3));
add(id(i, 3), fid(i, 0)); add(fid(i, 3), id(i, 0));
add(id(i, 1), fid(i, 2)); add(fid(i, 1), id(i, 2));
add(id(i, 2), fid(i, 1)); add(fid(i, 2), id(i, 1));
for(int j = 0; j <= 3; ++ j) if(!t(i, j, d)) add(id(i, j), fid(i, j));
}
for(int i = 1; i < n; ++ i)
for(int a = 0; a <= 3; ++ a)
for(int j = i + 1; j <= n; ++ j)
for(int b = 0; b <= 3; ++ b)
if(f(i, a, j, b, d)) add(id(i, a), fid(j, b)), add(id(j, b), fid(i, a));
for(int i = 1; i <= n * 8; ++ i) if(!dfn[i]) tarjan(i);
for(int i = 1; i <= n * 4; ++ i) if(sd[i] == sd[i + n * 4]) return 0;
return 1;
}
signed main()
{
freopen("lucky.in", "r", stdin);
freopen("lucky.out", "w", stdout);
scanf("%lld%lld%lld", &n, &w, &h); w *= 2; h *= 2;
for(int i = 1; i <= n; ++ i) scanf("%lld", &x[i]), x[i] *= 2;
for(int i = 1; i <= n; ++ i) scanf("%lld", &y[i]), y[i] *= 2;
int l = 0, r = min(w, h), ans = 0;
while(l <= r)
{
int mid = l + r >> 1;
check(mid) ? l = mid + 1, ans = mid : r = mid - 1;
}
printf("%.6lf", 1.0 * ans);
fclose(stdin);
fclose(stdout);
return 0;
}
B. 山河入梦来
给定一个 \(n\times n\) 的矩阵，第 \(i\) 行的第 \(L_i\) 到 \(R_i\) 列为 \(1\)，其余为 \(0\)。求行列式。\(n\le 10^5\)
求行列式可以转换成上三角矩阵。
对于 \(L_i\) 相同的一些值，我们会保留 \(R_i\) 最小的行。使用可并堆维护。
点击查看代码
#include <cstdio>
#include <iostream>
using namespace std;
const int N = 1000005;
int T, n, ans, p[N], rt[N], r[N];
struct Node {int x, dep, ls, rs;} t[N];
int merge(int u, int v)
{
if(!u || !v) return u + v;
if(t[u].x > t[v].x) swap(u, v);
t[u].rs = merge(t[u].rs, v);
if(t[t[u].ls].dep < t[t[u].rs].dep) swap(t[u].ls, t[u].rs);
t[u].dep = t[t[u].rs].dep + 1; return u;
}
void work()
{
scanf("%d", &n); ans = 1;
for(int i = 1; i <= n; ++ i) rt[i] = 0, p[i] = r[i] = i;
for(int i = 1, l, r; i <= n; ++ i)
scanf("%d%d", &l, &r), t[i] = {r, 1, 0, 0}, rt[l] = merge(rt[l], i);
// p[i] 左端点为 i 的行数
// r[i] 第 i 行的左端点
for(int i = 1; i <= n; ++ i)
{
if(!rt[i]) {ans = 0; break;}
int x = rt[i], w = t[rt[i]].x, y;
if(p[x] != i)
p[y = r[i]] = p[x], r[p[x] = i] = x, r[p[y]] = y, ans *= -1;
rt[i] = merge(t[rt[i]].ls, t[rt[i]].rs);
if(rt[i] && t[rt[i]].x == w) {ans = 0; break;}
rt[w + 1] = merge(rt[w + 1], rt[i]);
}
if(ans == 0) printf("tie\n");
if(ans > 0) printf("pp\n");
if(ans < 0) printf("xx\n");
}
int main()
{
freopen("river.in", "r", stdin);
freopen("river.out", "w", stdout);
scanf("%d", &T); while(T --) work();
return 0;
}
C. 遇见
有一个长度为 \(n\) 的序列，每个元素是 \(1\) 到 \(m\) 的整数。有多少种序列满足在任何一个长度为 \(k\) 的子区间线性相关。
当 \(k > \log_2m\) 时，答案为 \(m^n\)。
当 \(k \le \log_2m\) 时，条件可以转换为任何一个线性无关的子区间的长度都 \(<k\)。
设 \(f_{i,j}\) 表示前 \(i\) 个数，后 \(j\) 个数为极长线性无关的序列个数。
如果新加入一个数没有缩短之前的线性无关区间长度，即与后 \(j\) 个数线性无关。
\(f_{i+1,j+1}+=f_{i,j}\times(m+1-2^j)\)
如果新加入的数与后 \(j\) 个数线性相关，与后 \(t\) 个数线性无关。
\(f_{i+1,t+1}+=f_{i,j}\times2^t\)
答案为 \(\sum_{i<k}f_{n,i}\)
使用矩阵快速幂优化。
点击查看代码
#include <cstdio>
#include <cstring>
#include <algorithm>
#define int long long
using namespace std;
const int mod = 1000000007;
int n, m, k, ans;
struct Matrix {int n, m, c[35][35];} A, B;
int Pow(int x, int y)
{
int r = 1; for(; y; y >>= 1, x = x * x % mod)
if(y & 1) r = r * x % mod; return r;
}
Matrix operator * (Matrix x, Matrix y)
{
Matrix a; memset(a.c, 0, sizeof(a.c));
for(int i = 1; i <= x.n; ++ i)
for(int j = 1; j <= x.m; ++ j)
for(int k = 1; k <= y.m; ++ k)
a.c[i][k] = (a.c[i][k] + x.c[i][j] * y.c[j][k] % mod) % mod;
a.n = x.n; a.m = y.m; return a;
}
Matrix operator ^ (Matrix x, int y)
{
Matrix r; r.n = r.m = x.n; memset(r.c, 0, sizeof(r.c));
for(int i = 1; i <= r.n; ++ i) r.c[i][i] = 1;
for(; y; y >>= 1, x = x * x) if(y & 1) r = r * x; return r;
}
signed main()
{
freopen("meet.in", "r", stdin);
freopen("meet.out", "w", stdout);
scanf("%lld%lld%lld", &n, &m, &k);
if(k > __lg(m + 1)) {printf("%lld", Pow(m, n)); return 0;}
A.n = 1; A.m = B.n = B.m = k - 1; A.c[1][1] = m;
for(int i = 1; i < k; ++ i)
{
if(i + 1 < k)
B.c[i][i + 1] = (m + 1 - Pow(2, i) + mod) % mod;
for(int j = 1; j <= i; ++ j)
B.c[i][j] = Pow(2, j - 1);
}
A = A * (B ^ n - 1);
for(int i = 1; i < k; ++ i) ans = (ans + A.c[1][i]) % mod;
printf("%lld", ans);
return 0;
}
浙公网安备 33010602011771号