A. 第一题

儿子遍历顺序按深度由小到大即可


B. 第二题

二分最小值,以点权作为初始距离跑最长路即可

直接用大根堆跑 \(dij\) 会 \(T\),考虑初始权值可以处理,且边权一定,用类似蚯蚓的方法开两个队列,一开始都在第一个队列,松弛后的点放第二个,容易发现两个队列都是单调递减的,比较队头即可

也可以直接 \(spfa\) 到不能松弛为止,由于图的特殊性不会更新很多次

代码实现
#include<bits/stdc++.h>
using namespace std;
#define int long long
const int maxn=1e7+5;
const int maxm=5e7+5;
int n,m,k,val[maxn],dis[maxn],hd[maxn],cnt;
bool vis[maxn];
int movx[10]={0,0,-1,-1,1,1};
int movy[10]={-1,1,0,1,-1,0};
int read(){
int x=0,f=1;
char ch=getchar();
while(!isdigit(ch)){
if(ch=='-')f=-1;
ch=getchar();
}
while(isdigit(ch)){
x=x*10+ch-48;
ch=getchar();
}
return x*f;
}
struct Edge{
int nxt,to;
}edge[maxm];
void add(int u,int v){
edge[++cnt].nxt=hd[u];
edge[cnt].to=v;
hd[u]=cnt;
return ;
}
struct Node{
int dis,id;
Node(){}
Node(int x,int y):dis(x),id(y){}
bool operator < (const Node &x)const{
return dis>x.dis;
}
}p[maxn];
queue<int>q1,q2;
bool check(int limit){
for(int i=1;i<=n*m;i++){
dis[p[i].id]=p[i].dis;
q1.push(p[i].id);
}
memset(vis,0,sizeof vis);
while((!q1.empty())||(!q2.empty())){
int x=0,y=0,u=0;
if(!q1.empty())x=q1.front();
if(!q2.empty())y=q2.front();
if(dis[x]>dis[y]){
q1.pop();
u=x;
}
else q2.pop(),u=y;
//cout<<"hhh "<<u<<" "<<dis[u]<<endl;
if(vis[u])continue; vis[u]=true;
for(int i=hd[u];i;i=edge[i].nxt){
int v=edge[i].to;
if(dis[v]<dis[u]-limit){
dis[v]=dis[u]-limit;
//if(!vis[v])
q2.push(v);
//,vis[v]=true;
}
}
}
int sum=0;
for(int i=1;i<=n*m;i++){
sum+=dis[p[i].id]-p[i].dis;
if(sum>k)return false;
}
return true;
}
bool pd(int x,int y){
return x>=1&&x<=n&&y>=1&&y<=m;
}
int id(int x,int y){
return (x-1)*m+y;
}
signed main(){
// freopen("s.out","r",stdin);
// freopen("my.out","w",stdout);
n=read();
m=read();
k=read();
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++){
val[id(i,j)]=read();
for(int l=0;l<6;l++){
int x=i+movx[l];
int y=j+movy[l];
if(pd(x,y))add(id(i,j),id(x,y));
}
}
}
for(int i=1;i<=n*m;i++){
p[i].dis=val[i];
p[i].id=i;
}
sort(p+1,p+n*m+1);
dis[0]=-1;
// for(int i=1;i<=n*m;i++)cout<<p[i].dis<<" "<<p[i].id<<endl;;
int l=0,r=1000000000000;
while(l<r){
int mid=(l+r)/2;
if(check(mid))r=mid;
else l=mid+1;
}
cout<<l;
return 0;
}

C. 第三题

\(1\) 的个数从大到小依次考虑即可,用数位 \(dp\) 计算个数和总和

这道题是同时卡上下边界的数位 \(dp\)


D. 第四题

设 \(f[i][j]\) 表示前 \(i\) 个数,最大值为 \(j\) 的方案数,转移 \(f[i][j]=f[i-1][j]*j+f[i-1][j-1]\)

设 \(g[i][j]\) 表示从第 \(i\) 个数开始到结尾,开头初始最大值为 \(j\) 的方案数,转移 \(g[i][j]=g[i+1][j]*j+g[i+1][j+1]\)

观察答案形式,由于 \(\displaystyle x^2=\binom{x}{2}*2+1\)

那么强制一个数在数列里出现两次

数 \(j\) 在 \(i\) 位置第一次出现的贡献为:

\[f[i-1][j-1]*(g[i+1][j]+2*(n-i)*g[i+2][j])+\sum_{k>=j}f[i-1][k]*(g[i+1][k]+2*(n-i)*g[i+2][k])
\]

即讨论上一个数的大小即可

然后用前缀和优化为 \(n^2\)

代码实现
#include<bits/stdc++.h>
using namespace std;
#define int long long
const int maxn=3005;
int n,mod,f[maxn][maxn],g[maxn][maxn],ans,sum2[maxn],sum1[maxn][maxn],sum[maxn][maxn];
signed main(){
cin>>n>>mod;
f[0][0]=1;
// for(int i=0;i<=n;i++)f[0][i]=1;
for(int i=1;i<=n;i++){
for(int j=1;j<=i;j++){
f[i][j]=(f[i-1][j-1]+f[i-1][j]*j%mod)%mod;
// cout<<f[i][j]<<" ";
}
}
// g[n+1][n]=g[n+1][n]=1;
for(int i=1;i<=n;i++)g[n+1][i]=1;
for(int i=n;i>=1;i--){
for(int j=n;j>=1;j--){
g[i][j]=(g[i+1][min(j+1,n)]+g[i+1][j]*j%mod)%mod;
}
}
// cout<<g[5][5]<<endl;
for(int i=1;i<=n;i++){
for(int j=n;j>=1;j--){
sum[i][j]=(sum[i][j+1]+f[i-1][j]*g[i+1][j]%mod)%mod;
sum1[i][j]=(sum1[i][j+1]+g[i+2][j]*f[i-1][j]%mod)%mod;
}
}
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++){
ans=0;
ans=(sum[i][j]+2*(n-i)%mod*sum1[i][j]%mod)%mod;
// for(int k=j;k<=n;k++){
// ans=(ans+f[i-1][k]*(g[i+1][k]+2*(n-i)%mod*g[i+2][k]%mod)%mod)%mod;
// }
ans=(ans+f[i-1][j-1]*(g[i+1][j]+2*(n-i)%mod*g[i+2][j]%mod)%mod)%mod;
sum2[j]=(sum2[j]+ans)%mod;
}
}
for(int i=1;i<=n;i++)cout<<sum2[i]<<" ";
return 0;
}

noip模拟43的更多相关文章

  1. 8.18考试总结[NOIP模拟43]

    又挂了$80$ 好气哦,但要保持优雅.(草 T1 地衣体 小小的贪心:每次肯定从深度较小的点向深度较大的点转移更优. 模拟一下,把边按链接点的子树最大深度排序,发现实际上只有上一个遍历到的点是对当前考 ...

  2. Noip模拟43 2021.8.18

    T1 地一体 可以树形$dp$,但考场没写出来,只打了没正确性的贪心水了$30$ 然后讲题的时候B哥讲了如何正确的贪心,喜出望外的学习了一下 不难发现 每次士兵都会直接冲到叶子节点 从深的点再返回到另 ...

  3. [考试总结]noip模拟43

    这个题目出的还是很偷懒.... 第一题...第二题...第三题...四.... 好吧... 这几次考得都有些问题,似乎可能是有些疲惫,脑袋也是转不太动,考完总觉得自己是能力的问题,但是改一分钟之后会发 ...

  4. NOIP模拟 1

    NOIP模拟1,到现在时间已经比较长了.. 那天是6.14,今天7.18了 //然鹅我看着最前边缺失的模拟1,还是终于忍不住把它补上,为了保持顺序2345重新发布了一遍.. #   用  户  名   ...

  5. NOIP模拟17.9.22

    NOIP模拟17.9.22 前进![问题描述]数轴的原点上有一只青蛙.青蛙要跳到数轴上≥

  6. NOIP 模拟4 T2

    本题属于二和一问题 子问题相互对称 考虑对于问题一:知a求b 那么根据b数组定义式 显然能发现问题在于如何求dis(最短路) 有很多算法可供选择 dijsktra,floyed,bfs/dfs,spf ...

  7. noip模拟32[好数学啊]

    noip模拟32 solutions 真是无语子,又没上100,无奈死了 虽然我每次都觉得题很难,但是还是有好多上100的 战神都200多了,好生气啊啊啊 从题开始变难之后,我的时间分配越来越不均匀, ...

  8. 2021.9.17考试总结[NOIP模拟55]

    有的考试表面上自称NOIP模拟,背地里却是绍兴一中NOI模拟 吓得我直接文件打错 T1 Skip 设状态$f_i$为最后一次选$i$在$i$时的最优解.有$f_i=max_{j<i}[f_j+a ...

  9. NOIP模拟赛20161022

    NOIP模拟赛2016-10-22 题目名 东风谷早苗 西行寺幽幽子 琪露诺 上白泽慧音 源文件 robot.cpp/c/pas spring.cpp/c/pas iceroad.cpp/c/pas ...

随机推荐

  1. Spring Data Commons 远程命令执行漏洞(CVE-2018-1273)

    影响版本 Spring Framework 5.0 to 5.0.4 Spring Framework 4.3 to 4.3.14 poc https://github.com/zhzyker/exp ...

  2. 官宣|VS Code 1.59 发布!Terminal 可以跨窗口移动了!

    欢迎使用 2021 年 7 月版的 Visual Studio Code.我们希望您会喜欢此版本中的许多更新与改进,以下是其中的一些高亮: 扩展视图的改进 - 丰富的扩展详细信息悬停,新的运行状态标签 ...

  3. 文本生成器(AC自动机+dp)

    F. 文本生成器 内存限制:512 MiB 时间限制:1000 ms 标准输入输出 题目类型:传统 评测方式:文本比较 题目描述 原题来自:JSOI 2007 JSOI 交给队员 ZYX 一个任务,编 ...

  4. 我为什么选Markdown

    前沿说明:Yaml Front Matter MarkDown 目录 前沿说明:Yaml Front Matter 什么是MarkDown Markdown是一种轻量级标记语言, 它允许人们使用易读易 ...

  5. Android无障碍宝典-talkback

    http://geek.csdn.net/news/detail/93269 http://geek.csdn.net/news/detail/135867

  6. 20初识前端HTML(1)

    1 .HTML 1.1 网页的组成 文字 图片 链接 等元素构成.除了这些元素之外 网页中还可以包含音频 视频 等 1.2 WEB前端开发的流程 现在主流的开发流程: 前后端分离的开发模式. 美工:p ...

  7. C++//菱形继承 //俩个派生类继承同一个基类 //又有某个类同时继承俩个派生类 //成为 菱形继承 或者 钻石 继承//+解决

    1 //菱形继承 2 //俩个派生类继承同一个基类 3 //又有某个类同时继承俩个派生类 4 //成为 菱形继承 或者 钻石 继承 5 6 #include <iostream> 7 #i ...

  8. 使用VNC远程安装CentOS 7操作系统

    使用VNC远程安装CentOS 7操作系统 by 无若 数据中心一般都不在本地,如果希望重新安装系统,难道还要跑到数据中心...所以必须要有一种方式来远程解决这个问题. 目前CentOS 7主要使用的 ...

  9. [TensorFow2.0]-MNIST手写数字识别

    本人人工智能初学者,现在在学习TensorFlow2.0,对一些学习内容做一下笔记.笔记中,有些内容理解可能较为肤浅.有偏差等,各位在阅读时如有发现问题,请评论或者邮箱(右侧边栏有邮箱地址)提醒. 若 ...

  10. WPF(MVVM) 利用资源字典实现中英文动态切换

    1.首先新建两个字典文件en-us.xaml.zh-cn.xaml.定义中英文的字符串在这里面. 2.将两个资源字典添加到App.xaml中,这里注意下,因为两个字典中有同样字符,如果没有动态更改,默 ...