Skip to content

泉州一中信息学Blog

信息学奥赛技术分享博客

  • 登录
  • 小学oj
  • 中学oj
  • 测试页面1
  • Toggle search form

NOIP基地校Day1

Posted on 2023年11月10日2023年11月10日 By 陈骏翔 NOIP基地校Day1无评论

T1 朴素dp都会吧,上矩阵秒了,但其实复杂度O(极大常数*lognT)是错的,大样例四秒钟,没被卡掉只能说有点遗憾吧

#include <bits/stdc++.h>
#define rep(i, j, k) for(ll i = (j), _i = (k); i <= _i; i++)
#define drp(i, j, k) for(ll i = (j), _i = (k); i >= _i; i--)
#define ll long long
using namespace std;
#define fo(s) freopen(s ".in","r",stdin),freopen(s ".out","w",stdout)
//-------------------------------------//
ll read() {
    ll x = 0, f = 0; char ch;
    while(!isdigit(ch = getchar())) f |= ch == &#039;-&#039;;
    while(isdigit(ch)) x = (x<<3) + (x<<1) + (ch^48), ch = getchar();
    return f ? -x : x;
}
//-------------------------------------//
static const ll mod = 998244353;
struct Matrix {
    ll a[3][3], n, m;
    Matrix(ll N): n(N), m(N) { 
        rep(i, 1, N) rep(j, 1, N) a[i][j] = 0; rep(i, 1, N) a[i][i] = 1;
    }
    Matrix(ll N, ll M): n(N), m(M) { rep(i, 1, N) rep(j, 1, M) a[i][j] = 0; }
    Matrix() {}
    friend Matrix operator*(Matrix a, Matrix b) {
        Matrix c;
        c = Matrix(a.n, b.m);
        rep(i, 1, a.n) rep(j, 1, b.m) rep(k, 1, a.m) c[i][j] = (c[i][j] + a[i][k] * b[k][j] % mod) % mod; 
        return c;
    }
    ll* operator[](ll x) { return a[x]; }
};
Matrix I, A, B;
ll n, m, ans, mk;

ll qpow(ll x, ll k){ return k ? (ll)pow(x, k&1) * qpow(x * x % mod, k>>1) % mod : 1ll; }

void initMat() {
    I = Matrix(2), A = Matrix(2, 2), B = Matrix(1, 2);
    ll N = (n - 1 + mod) % mod;
    ll k = ((N * N % mod - 1) + mod) % mod;
    A[1][1] = k, A[1][2] = 0, A[2][1] = 2, A[2][2] = n * n % mod, B[1][1] = 1, B[1][2] = 1;
}
ll solve() {
    n = read() % mod, mk = m = read();
    initMat();
    ll inv = qpow(qpow(n * n % mod, m), mod - 2);
    while(mk) { 
        if(mk&1) I = I * A; 
        A = A * A, mk >>= 1; 
    }
    return (qpow(n * n % mod, m) - (B * I)[1][1] + mod) % mod * n % mod * inv % mod;
}
int main() 
{
    fo("poster");
    rep(T, 1, read()) printf("%lld\n",solve());
    return 0;
}

T2 dp[u][0/1/2/3]表示u这个点自己放路由器/只有儿子放路由器/只有父亲放路由器/儿子和父亲都有都有路由器。转移显然

#include <bits/stdc++.h>
#define rep(i, j, k) for(ll i = (j), _i = (k); i <= _i; i++)
#define drp(i, j, k) for(ll i = (j), _i = (k); i >= _i; i--)
#define ll long long
#define Nx 1000050
using namespace std;
#define fo(s) freopen(s ".in","r",stdin),freopen(s ".out","w",stdout)
//-------------------------------------//
ll read() {
    ll x = 0, f = 0; char ch;
    while(!isdigit(ch = getchar())) f |= ch == &#039;-&#039;;
    while(isdigit(ch)) x = (x<<3) + (x<<1) + (ch^48), ch = getchar();
    return f ? -x : x;
}
//-------------------------------------//
static const ll inf = 1e17;
ll n, a[Nx], f[Nx][4];// 0 self 1 son 2 father 3 both
vector<ll> G[Nx];
void dfs(ll u, ll fa) {
    f[u][0] = a[u];
    if(!G[u].size()) {
        f[u][1] = f[u][3] = inf; return ;
    }
    ll sum = 0, mn = inf, flag = 0;
    for(auto v : G[u]) {
        if(v == fa) continue;
        dfs(v, u);
        f[u][0] += min(f[v][0], min(f[v][2], f[v][3]));
        if(f[v][1] != inf) f[u][2] += f[v][1];
        else f[u][2] = inf;
        if(f[v][0] < f[v][1]) flag = 1, sum += f[v][0];
        else sum += f[v][1], mn = min(mn, f[v][0] - f[v][1]);
    }
    f[u][1] = f[u][3] = flag ? sum : sum + mn; 
}
int main() 
{
    fo("router");
    rep(i, 1, n = read()) a[i] = read();
    rep(i, 1, n - 1) {
        ll u = read(), v = read();
        G[u].push_back(v), G[v].push_back(u);
    }
    dfs(1, 0);
    printf("%lld", min(f[1][0], f[1][1]));
    return 0;
}

T3 & T4 不熟

训练日志

文章导航

Previous Post: 只因地校集训DAY1
Next Post: Nov. 10th 博客

发表回复 取消回复

要发表评论,您必须先登录。

2025年 6月
一 二 三 四 五 六 日
 1
2345678
9101112131415
16171819202122
23242526272829
30  
« 2月    

2024常州 Class Classic OI Problems Contest cqr的长乐集训2023 CZYZ LOC New Game NOI NOIP Password Protected PM_PK Preview Problems Retrospect Selfmade Qusetion STL The end Training Uneasy Problem 蒟蒻 通报

  • 训练日志
  • 链表
  • 入门
  • 模拟
  • dfs序
  • 并查集
  • spfa
  • 最小割
  • 矩阵树定理
  • 仙人掌
  • BSGS
  • 凸包
  • 回文自动机
  • 递推与动归
  • 堆
  • 莫队算法
  • ST表
  • Treap
  • 树套树
  • 可持久化线段树
  • 初赛
  • 搜索
  • 贪心
  • 深度优先搜索
  • 欧拉图
  • dijkstra
  • 费用流
  • 哈夫曼树
  • kruskual
  • 置换
  • 旋转卡壳
  • KMP
  • 区间动归
  • STL
  • 链表
  • 可并堆
  • sply
  • 主席树
  • 可持久化字典树
  • 算法
  • 动态规划
  • 构造
  • 广度优先搜索
  • 最短路
  • floyd
  • 最大流
  • 虚树
  • prim
  • 筛法
  • 半平面交
  • 字典树
  • 背包动归
  • 基础数据结构
  • 分块
  • 线段树
  • 替罪羊树
  • K-DTree
  • 图论
  • 二分法
  • 迭代搜索
  • 拓扑排序
  • 有上下界网络流
  • 生成树
  • 快速幂
  • 后缀数组
  • 树形动归
  • 哈希表
  • 中级数据结构
  • 平衡树
  • 可持久化数据结构
  • 数据结构
  • 三分法
  • 启发式搜索
  • 图的连通
  • 点分治
  • 博弈论
  • AC自动机
  • 状压动归
  • 单调栈
  • 树状数组
  • 高级数据结构
  • OI资料
  • 数学
  • 高精度
  • 差分约束
  • 树上倍增
  • 素数测试
  • 后缀自动机
  • 数位动归
  • 单调队列
  • 新闻
  • 几何
  • 随机化
  • 二分图染色
  • 树链剖分
  • 欧拉函数
  • manacher
  • 斜率优化
  • 离线处理
  • 信息学奥赛学长风采
  • 字符串
  • 二分图匹配
  • prufer编码
  • 卡特兰数
  • 密码学
  • 决策单调
  • 赛后总结
  • 其他
  • 2-SAT
  • 最近公共祖先
  • 矩阵乘法
  • 记忆化搜索
  • 网络流
  • Link cut tree
  • 排列组合
  • 树
  • 高斯消元
  • 乘法逆元
  • 容斥原理
  • 调和级数
  • 概率与期望
  • 模线性方程组
  • 莫比乌斯反演
  • 快速傅里叶变换
  • 扩展欧几里德
  • 最大公约数与最小公倍数

近期文章

  • DP杂题
  • 2025年2月13日模拟赛
  • HLOJ-TEST ROUND 4-T1/T2(构造)- 3
  • HLOJ-TEST ROUND 4-T1/T2(构造)- 2
  • HLOJ-TEST ROUND 4-T1/T2(构造)- 1

近期评论

归档

  • 2025年2月
  • 2025年1月
  • 2024年11月
  • 2024年10月
  • 2024年9月
  • 2024年8月
  • 2024年7月
  • 2024年3月
  • 2024年2月
  • 2024年1月
  • 2023年12月
  • 2023年11月
  • 2023年10月
  • 2023年9月
  • 2023年8月
  • 2023年7月
  • 2023年3月
  • 2023年2月
  • 2023年1月
  • 2022年12月

Copyright © 2025 泉州一中信息学Blog.

Powered by PressBook WordPress theme