8.13考试总结(NOIP模拟38)[a·b·c]
重要的不是你做了多少事，而是你放了多少心思进去。
T1 a
解题思路
总结一下，是双指针运用不够熟练（zxb笑了笑）。
其实这个题是可以用树状数组卡过的(众所周知我是一个正直的人)，但是一定是要打正解的。
树状数组比较好像，就和 入阵曲 一样只不过这个维护的是个范围。
因此需要树状数组维护前缀和，时间复杂度就多了一个 log
首先，这个数据范围明显是让我们 \(n^2m\) 跑过。
因此先枚举矩形的上边界，接着枚举列。
然后双指针扫比当前左边界为大矩形左边界，右边界为当前扫到列的矩形大的部分刚好在 \([l,r]\) 的列
注意卡一下边界，对于 \(l=0,r=n\times m\) 的需要特殊判断一下。
code
57pts 树状数组
#include<bits/stdc++.h>
#define int long long
#define ull unsigned long long
#define f() cout<<"Pass"<<endl
using namespace std;
inline int read()
{
int x=0,f=1;char ch=getchar();
while(ch>'9'||ch<'0'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
return x*f;
}
const int N=40,M=5e4+10,K=N*M;
int n,m,ans,li,ri,s[N][M],pre[N][M],tre[K];
char ch[M];
int lowbit(int x)
{
return x&(-x);
}
void insert(int x,int num)
{
for(int i=x+1;i<=pre[n][m]+1;i+=lowbit(i))
tre[i]+=num;
}
int query(int x)
{
if(x<0) return 0;
int sum=0;
for(int i=x+1;i;i-=lowbit(i))
sum+=tre[i];
return sum;
}
signed main()
{
scanf("%lld%lld",&n,&m);
for(int i=1;i<=n;i++)
{
scanf("%s",ch+1);
for(int j=1;j<=m;j++)
s[i][j]=ch[j]-'0';
}
scanf("%lld%lld",&li,&ri);
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
pre[i][j]=pre[i-1][j]+pre[i][j-1]-pre[i-1][j-1]+s[i][j];
ri=min(ri,pre[n][m]);
insert(0,1);
for(int i=0;i<n;i++)
for(int j=i+1;j<=n;j++)
{
for(int k=1;k<=m;k++)
{
int tmp=pre[j][k]-pre[i][k];
//cout<<i<<' '<<j<<' '<<k<<" "<<tmp<<' '<<query(tmp-li)<<' '<<query(tmp-ri-1)<<endl;
ans+=query(tmp-li)-query(tmp-ri-1);
//if(tmp)
insert(pre[j][k]-pre[i][k],1);
}
//memset(tre,0,sizeof(tre));
for(int k=1;k<=m;k++)
if(pre[j][k]-pre[i][k]>=0)
insert(pre[j][k]-pre[i][k],-1);
}
printf("%lld",ans);
return 0;
}
正解
#include<bits/stdc++.h>
#define int long long
#define ull unsigned long long
#define f() cout<<"Pass"<<endl
using namespace std;
inline int read()
{
int x=0,f=1;char ch=getchar();
while(ch>'9'||ch<'0'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
return x*f;
}
const int N=40,M=5e4+10,K=N*M;
int n,m,ans,li,ri,s[N][M],pre[N][M],l[M],r[M];
char ch[M];
signed main()
{
scanf("%lld%lld",&n,&m);
for(int i=1;i<=n;i++)
{
scanf("%s",ch+1);
for(int j=1;j<=m;j++)
s[i][j]=ch[j]-'0';
}
scanf("%lld%lld",&li,&ri);
if(!li&&ri==m*n)
{
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
ans+=(n-i+1)*(m-j+1);
printf("%lld",ans);
return 0;
}
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
pre[i][j]=pre[i-1][j]+pre[i][j-1]-pre[i-1][j-1]+s[i][j];
ri=min(ri,pre[n][m]);
for(int i=0;i<n;i++)
{
memset(l,0,sizeof(l));
memset(r,0,sizeof(r));
for(int j=0;j<m;j++)
{
for(int k=i+1;k<=n;k++)
{
int num=pre[k][j]-pre[i][j];
while(l[k]<m&&pre[k][l[k]]-pre[i][l[k]]-num<li) l[k]++;
while(r[k]<m&&pre[k][r[k]+1]-pre[i][r[k]+1]-num<=ri) r[k]++;
if(l[k]==m&&pre[k][l[k]]-pre[i][l[k]]-num<li) continue;
ans+=r[k]-l[k]+1;
}
}
}
printf("%lld",ans);
return 0;
}
T2 b
解题思路
显然是枚举 gcd 。
接下来对于重复的部分直接暴力容斥算回去。
问题转换为： 对于\(i \in [1,10^5]\)求多少种选择方案使得选的所有数均为 i 的倍数。
预处理出每一行中每一个数字的个数 记为 \(cnt_{i,j}\)
最后的答案就是 \(\prod\limits_{i=1}^{n}(cnt_{i,j}+1) -1\) 种方案。
毕竟要算上 0 的情况，去掉都是 0 的情况。
code
#include<bits/stdc++.h>
#define int long long
#define ull unsigned long long
#define f() cout<<"Pass"<<endl
using namespace std;
inline int read()
{
int x=0,f=1;char ch=getchar();
while(ch>'9'||ch<'0'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
return x*f;
}
const int N=30,M=1e5+10,mod=1e9+7;
int n,m,ans,mx,s[N][M],cnt[N][M],all[M];
signed main()
{
n=read(); m=read();
for(int i=1,x;i<=n;i++)
for(int j=1;j<=m;j++)
{
x=read();
s[i][x]++;
mx=max(mx,x);
}
for(int i=1;i<=n;i++)
for(int j=1;j<=mx;j++)
for(int k=1;k*j<=mx;k++)
cnt[i][j]+=s[i][j*k];
for(int i=mx;i>=1;i--)
{
bool flag=false;
for(int j=1;j<=n;j++)
flag|=(cnt[j][i]!=0);
if(!flag) goto V;
all[i]=1;
for(int j=1;j<=n;j++)
all[i]=all[i]*(cnt[j][i]+1)%mod;
all[i]=(all[i]-1+mod)%mod;
for(int j=2;i*j<=mx;j++)
all[i]=(all[i]-all[i*j]%mod+mod)%mod;
V:;
}
for(int i=1;i<=mx;i++)
ans=(ans+all[i]*i%mod)%mod;
printf("%lld",ans);
return 0;
}
T3 c
解题思路
点分治。
首先可以发现其实这就是一棵树，每一条边可能会有多种颜色。
然后每一条边最多留下三种颜色就不会影响我们的答案（贪心，抽屉原理）。
在对于颜色去重之后离线点分治。
每层跑一遍，\(dp_{i,j,k}\)：分治中心 focus 到 i 的路径中，第一条边的颜色为j，最后一条边颜色为 k 的最优解。
首先是点权下放，然后暴力枚举 focus，x，y 的四条边的眼色就好了。
code
#include<bits/stdc++.h>
#define int long long
#define ull unsigned long long
#define f() cout<<"Pass"<<endl
using namespace std;
inline int read()
{
int x=0,f=1;char ch=getchar();
while(ch>'9'||ch<'0'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
return x*f;
}
const int N=1e5+10,M=3e5+10,Q=1e5+10;
int n,m,qus,focus,root,tot,ans[Q],siz[N],rt[N],w[N][5],cnt[N],f[N][5][5];
bool vis[N];
vector<pair<int,int> > v[N],q[N];
vector<pair<int,vector<int> > > e[N];
void pre_dfs(int x,int fat)
{
sort(v[x].begin(),v[x].end());
v[x].erase(unique(v[x].begin(),v[x].end()),v[x].end());
for(int i=0,j;i<v[x].size();i=j)
{
int to=v[x][i].first;
for(j=i+1;j<v[x].size();j++)
if(v[x][j].first!=to)
break;
if(to==fat) continue;
e[x].push_back(make_pair(to,vector<int>()));
e[to].push_back(make_pair(x,vector<int>()));
for(int k=i;k<j&&k<i+3;k++)
{
e[x].back().second.push_back(v[x][k].second);
e[to].back().second.push_back(v[x][k].second);
}
pre_dfs(to,x);
}
}
void get_focus(int x,int all,int fat)
{
siz[x]=1;
int maxn=0;
for(int i=0;i<e[x].size();i++)
{
int to=e[x][i].first;
if(to==fat||vis[to]) continue;
get_focus(to,all,x);
maxn=max(maxn,siz[to]);
siz[x]+=siz[to];
}
maxn=max(all-siz[x],maxn);
if(maxn<tot)
tot=maxn,focus=x;
}
void work(int x,int fat)
{
siz[x]=1;
for(int i=0;i<e[x].size();i++)
{
int to=e[x][i].first;
if(to==fat||vis[to]) continue;
for(int j=0;j<e[x][i].second.size();j++)
w[to][j+1]=e[x][i].second[j];
cnt[to]=e[x][i].second.size();
for(int l=1;l<=cnt[root];l++)
for(int j=1;j<=cnt[to];j++)
for(int k=1;k<=cnt[x];k++)
if(w[x][k]==w[to][j]) f[to][l][j]=max(f[to][l][j],f[x][l][k]);
else f[to][l][j]=max(f[to][l][j],f[x][l][k]+1);
work(to,x);
siz[x]+=siz[to];
}
}
void query(int x,int fat)
{
for(int i=0;i<q[x].size();i++)
{
int num=q[x][i].second,id=q[x][i].first;
if(!rt[num]&&num!=focus) continue;
if(num==focus)
{
for(int j=1;j<=cnt[root];j++)
for(int k=1;k<=cnt[x];k++)
ans[id]=max(ans[id],f[x][j][k]);
continue;
}
for(int j=1;j<=cnt[root];j++)
for(int k=1;k<=cnt[rt[num]];k++)
for(int a=1;a<=cnt[x];a++)
for(int b=1;b<=cnt[num];b++)
if(w[rt[num]][k]==w[root][j]) ans[id]=max(ans[id],f[x][j][a]+f[num][k][b]-1);
else ans[id]=max(ans[id],f[x][j][a]+f[num][k][b]);
}
for(int i=0;i<e[x].size();i++)
{
int to=e[x][i].first;
if(to==fat||vis[to]) continue;
query(to,x);
}
}
void clear(int x,int fat)
{
rt[x]=0;
memset(f[x],0,sizeof(f[x]));
for(int i=0;i<e[x].size();i++)
{
int to=e[x][i].first;
if(to==fat||vis[to]) continue;
clear(to,x);
}
}
void mark(int x,int fat)
{
rt[x]=root;
for(int i=0;i<e[x].size();i++)
{
int to=e[x][i].first;
if(vis[to]||to==fat) continue;
mark(to,x);
}
}
void solve(int x,int all)
{
tot=all;
get_focus(x,all,0);
clear(focus,0);
vis[focus]=true;
int tmp=focus;
for(int i=0;i<e[tmp].size();i++)
{
int to=e[tmp][i].first;
if(vis[to]) continue;
cnt[to]=e[tmp][i].second.size();
for(int j=0;j<e[tmp][i].second.size();j++)
w[to][j+1]=e[tmp][i].second[j],f[to][j+1][j+1]=1;
root=to;
work(to,tmp);
query(to,tmp);
mark(to,tmp);
}
clear(focus,0);
for(int i=0;i<e[tmp].size();i++)
{
int to=e[tmp][i].first;
if(vis[to]) continue;
solve(to,siz[to]);
}
}
signed main()
{
n=read(); m=read();
for(int i=1,x,y,val;i<=m;i++)
{
x=read(); y=read(); val=read();
v[x].push_back(make_pair(y,val));
v[y].push_back(make_pair(x,val));
}
pre_dfs(1,0);
qus=read();
for(int i=1,x,y;i<=qus;i++)
{
x=read(); y=read();
q[x].push_back(make_pair(i,y));
q[y].push_back(make_pair(i,x));
}
solve(1,n);
for(int i=1;i<=qus;i++)
printf("%lld\n",ans[i]);
return 0;
}
8.13考试总结(NOIP模拟38)[a·b·c]的更多相关文章
- 2021.8.13考试总结[NOIP模拟38]
T1 a 入阵曲.枚举矩形上下界,之后从左到右扫一遍.用树状数组维护前缀和加特判可以$A$,更保险要脸的做法是双指针扫,因为前缀和单调不减. $code:$ 1 #include<bits/st ...
- 2021.9.13考试总结[NOIP模拟52]
T1 路径 考虑每一位的贡献,第$i$位每$2^i$个数会变一次,那么答案为$\sum_{i=1}^{log_2n} \frac{n}{2^i}$. $code:$ 1 #include<bit ...
- 5.23考试总结(NOIP模拟2)
5.23考试总结(NOIP模拟2) 洛谷题单 看第一题第一眼,不好打呀;看第一题样例又一眼,诶,我直接一手小阶乘走人 然后就急忙去干T2T3了 后来考完一看,只有\(T1\)骗到了\(15pts\)[ ...
- 6.17考试总结(NOIP模拟8)[星际旅行·砍树·超级树·求和]
6.17考试总结(NOIP模拟8) 背景 考得不咋样,有一个非常遗憾的地方:最后一题少取膜了,\(100pts->40pts\),改了这么多年的错还是头一回看见以下的情景... T1星际旅行 前 ...
- 5.22考试总结(NOIP模拟1)
5.22考试总结(NOIP模拟1) 改题记录 T1 序列 题解 暴力思路很好想,分数也很好想\(QAQ\) (反正我只拿了5pts) 正解的话: 先用欧拉筛把1-n的素数筛出来 void get_Pr ...
- 2021.9.17考试总结[NOIP模拟55]
有的考试表面上自称NOIP模拟,背地里却是绍兴一中NOI模拟 吓得我直接文件打错 T1 Skip 设状态$f_i$为最后一次选$i$在$i$时的最优解.有$f_i=max_{j<i}[f_j+a ...
- [考试总结]noip模拟23
因为考试过多,所以学校的博客就暂时咕掉了,放到家里来写 不过话说,vscode的markdown编辑器还是真的很好用 先把 \(noip\) 模拟 \(23\) 的总结写了吧.. 俗话说:" ...
- 「考试」noip模拟9,11,13
9.1 辣鸡 可以把答案分成 每个矩形内部连线 和 矩形之间的连线 两部分 前半部分即为\(2(w-1)(h-1)\),后半部分可以模拟求(就是讨论四种相邻的情况) 如果\(n^2\)选择暴力模拟是有 ...
- Noip模拟38 2021.8.13
T1 a 跟入阵曲很像,但是忘记入阵曲这题的思路是什么了 这里再提一下,入阵曲是子矩阵和是$k$的倍数,这道题目是子矩阵和是在一段区间内$[L,R]$ 因为这道题$n$特别小,$m$较大,考虑复杂度为 ...
- [考试总结]noip模拟13
因为最近考试频繁,所以咕掉了好长时间... 淦,刚说完又来一场... 先咕了,等以后有时间再写.... 回来了... 首先看到这个题目们,感觉就不存好意... 然后开始开 \(T1\). 只能蒻蒻地按 ...
随机推荐
- xml转voc，voc转coco，coco转yolo，coco划分，coco检查，yolo检查，coco可视化
平常用coco格式的数据集比较多,所有这里整合一下数据集相关的常用的脚本. pycocotools安装 这个非常重要,因为处理coco数据集时,用pycocotools包非常方便. 自行搜索一下怎么安 ...
- vue-cli4.0 （vue3.0 的脚手架）
前言: 这个搭建脚手架的话实际是我们创建一个新项目的第一步,当然,现在脚手架4.0都出来了,经过使用后发现跟我们之前的3.0使用方法是答题一样的,其中用vue-cli3.0来搭建我们的项目的话又分为两 ...
- 修改中文、英文参考文献在文末列表中的顺序：EndNote
本文介绍在EndNote软件中,使得参考文献按照语种排列,中文在前.英文在后的方法. 前期我们在EndNote参考文献格式Output Styles界面介绍一文中,详细介绍了文献管理软件End ...
- 力扣142(Java)-环形链表Ⅱ（中等）
题目: 给定一个链表的头节点 head ,返回链表开始入环的第一个节点. 如果链表无环,则返回 null. 如果链表中有某个节点,可以通过连续跟踪 next 指针再次到达,则链表中存在环. 为了表示 ...
- 力扣614(MySQL)-二级关注者（中等）
题目: 在 facebook 中,表 follow 会有 2 个字段: followee, follower ,分别表示被关注者和关注者. 请写一个 sql 查询语句,对每一个关注者,查询关注他的关注 ...
- CF1913C Game with Multiset 题解
[题目描述] 你有一个空的多重集,你需要处理若干下列询问: ADD $ x $:加入一个数值为 $ 2^x $ 的元素到该多重集. GET $ w $:判断是否存在一个该多重集的子集,使得这个子集的所 ...
- KubeVela 正式开源：一个高可扩展的云原生应用平台与核心引擎
美国西部时间 2020 年 11 月 18 日,在云原生技术"最高盛宴"的 KubeCon 北美峰会 2020 上,CNCF 应用交付领域小组(CNCF SIG App Deliv ...
- 网不好怎么办？TLS握手带宽直降80%，BabaSSL是怎么做到的？| 龙蜥技术
简介:为了保障数据的安全性,客户端会先和服务器进行 TLS 握手,有什么办法可以减少 TLS 握手的带宽消耗呢? 编者按:BabaSSL 是一款开源的密码库产品,在 GitHub 和龙蜥社区开源,并 ...
- 记 Win8.1 某应用渲染抛出 OutOfMemoryException 异常及修复方法
本文记录某个应用在某台 Windows 8.1 x86 系统上,运行时抛出 OutOfMemoryException 异常,启动失败.应用程序能启动,但是在第一次碰到渲染时,就发现渲染初始化失败,从而 ...
- 2019-4-29-Roslyn-将这个文件放在你的项目文件夹，无论哪个控制台项目都会输出林德熙是逗比...
title author date CreateTime categories Roslyn 将这个文件放在你的项目文件夹,无论哪个控制台项目都会输出林德熙是逗比 lindexi 2019-4-29 ...