NOIP模拟测试23「mine·water·gcd」
mine
题解
一道比较水的dp
考试因为初始化挂掉了只有$80$分
代码有注释
#include<bits/stdc++.h>
using namespace std;
//无脑dp
#define ll long long
#define A 2222222
const ll mod=1e9+7;
//设f[x][5][2]为当前是第x个格子,当前填什么 0前面填1后面填
ll f[A][5][2];
//0表示0雷,1表示一个雷 2表示两个雷 3*
ll n,ans=0;
char ch[A];
void debuger(){
for(ll i=1;i<=n;i++,puts(""))
for(ll j=0;j<=3;j++){
if(j==1){
printf("f[%lld][1][0]=%lld f[%lld][1][1]=%lld ",
i,f[i][1][0],i,f[i][1][1]);
}
else printf("f[%lld][%lld]=%lld ",i,j,f[i][j][0]);
}
}
int main(){
scanf("%s",ch+1);
n=strlen(ch+1);
if(ch[1]=='*'){
f[0][1][1]=1;
}
else if(ch[1]=='?'){
f[0][1][1]=1;
f[0][0][0]=1;
}
else if(ch[1]=='1'){
f[0][0][0]=1;
}
else if(ch[1]=='0'){
f[0][0][0]=1;
}
for(ll x=1;x<=n;x++){
if(ch[x]=='?'){
f[x][3][0]=(f[x-1][1][1]+f[x-1][2][0]+f[x-1][3][0])%mod;
f[x][0][0]=(f[x-1][1][0]+f[x-1][0][0])%mod;
f[x][2][0]=(f[x-1][3][0])%mod;
f[x][1][0]=(f[x-1][3][0])%mod;
f[x][1][1]=(f[x-1][1][0]+f[x-1][0][0])%mod;
}
else if(ch[x]=='*'){
f[x][3][0]=(f[x-1][1][1]+f[x-1][2][0]+f[x-1][3][0])%mod;
//若当前格子是雷,前面格子可以是1且雷在后,两个雷标记,雷
}
else if(ch[x]=='0'){
f[x][0][0]=(f[x-1][1][0]+f[x-1][0][0])%mod;
//若当前格子是0雷标记,前面格子可以是0雷标记或1雷但雷在前
}
else if(ch[x]=='2'){
f[x][2][0]=(f[x-1][3][0])%mod;
//若当前格子是2雷标记,前面格子必须是雷
}
else if(ch[x]=='1'){
f[x][1][0]=f[x-1][3][0];
//若雷在前,前面格子必须是雷
f[x][1][1]=(f[x-1][1][0]+f[x-1][0][0])%mod;
//若雷在后,前面格子可以是1雷标记但雷在前,或是无雷标记
}
}
// debuger();
ans=(ans+f[n][1][0]+f[n][0][0])%mod;
ans=(ans+f[n][3][0])%mod;
printf("%lld\n",ans);
}
water
题解
一道比较水的贪心
或是一道比较难的图论
看你怎么做了
下面是贪心做法
不断地去找最低的那个小块进行灌水,这是肯定的。
由于需要一直找最小值,所以用到最小堆来进行优化。
初始化,边界是肯定不能积水的,因为肯定会流出来。所以在一开始所有的边界都设置为已访问标记。将边界上所有的点都加进堆里,并不断地调整堆。
在边界中取最低的那一个小块进行扩展,记录这个最小的节点高度为t ,扩展到的内部节点如果比t小,那么就灌水,也就是将扩展到的节点高度上升到t,不要忘了统计ans。并同时将上升后的点加进堆。如果扩展到的点比t高,也就相当于一个边界,将这个点直接加进堆。
注意如果为$-$那么它一定是要积水的而且至少积水到0
代码
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define A 510
ll h[A][A],a[A][A];
bool vis[A][A];
struct node{
ll x,y,high;
friend bool operator < (const node &a,const node &b){
return a.high>b.high;
}
}; priority_queue<node> q;
ll n,m;
const ll nowx[5]={0,1,-1,0,0};
const ll nowy[5]={0,0,0,1,-1};
bool pan(ll x,ll y){
if(x<1||x>n||y<1||y>m||vis[x][y]) return 0;
return 1;
}
void dfs(ll x,ll y,ll v){
vis[x][y]=1;
if(a[x][y]>v){
node o;
o.x=x,o.y=y,o.high=a[x][y];
q.push(o);
return ;
}
h[x][y]+=v-a[x][y];
for(ll i=1;i<=4;i++){
ll xnow=x+nowx[i],ynow=y+nowy[i];
if(pan(xnow,ynow))
dfs(xnow,ynow,v);
}
}
int main(){
scanf("%lld%lld",&n,&m);
for(ll i=1;i<=n;i++)
for(ll j=1;j<=m;j++){
scanf("%lld",&a[i][j]);
if(a[i][j]<0){
h[i][j]=0-a[i][j];
a[i][j]=0;
}
if(i==1||j==1||i==n||j==m){
node b;
b.x=i,b.y=j,b.high=a[i][j];
q.push(b);
vis[i][j]=1;
}
}
while(!q.empty()){
node o=q.top();
q.pop();
ll x=o.x,y=o.y,high=o.high;
for(ll i=1;i<=4;i++){
ll xnow=x+nowx[i],ynow=y+nowy[i];
if(pan(xnow,ynow))
dfs(xnow,ynow,high);
}
}
for(ll i=1;i<=n;i++,puts(""))
for(ll j=1;j<=m;j++){
printf("%lld ",h[i][j]);
}
}
gcd
题解
奇妙的莫比乌斯,等我强了再补
代码
#include<bits/stdc++.h>
using namespace std;
#define ll int
#define A 1001010
ll vis[A],mu[A],prime[A];
long long g[A],s[A],su[A];
ll get_[A];
long long ans,cnt=0;
void get__mu(ll n){
mu[1]=1;
for(ll i=2;i<=500000;i++){
if(!vis[i]){
prime[++cnt]=i,mu[i]=-1;
}
for(ll j=1;j<=cnt&&prime[j]*1ll*i<=500000;j++){
vis[prime[j]*1ll*i]=1;
if(i%prime[j]==0) break;
else mu[i*1ll*prime[j]]=-mu[i];
}
}
}
ll n,m;
void work(ll x,ll d){
// printf("x=%lld\n",x);
for(ll i=1;i<=sqrt(x);i++)
if(!(x%i)){
s[i]+=d;
// printf("x=%lld i=%lld \n",x,i);
if(i*1ll*i!=x)s[x/i]+=d;
ans-=mu[i]*1ll*g[i];
if(i*1ll*i!=x)ans-=mu[x/i]*1ll*g[x/i];
g[i]=s[i]*1ll*(s[i]-1)/2;
if(i*1ll*i!=x)g[x/i]=s[x/i]*1ll*(s[x/i]-1)/2;
ans+=mu[i]*1ll*g[i];
if(i*1ll*i!=x)ans+=mu[x/i]*1ll*g[x/i];
}
// for(ll i=2;i<=x;i++){
// printf("s[%lld]=%lld\n",i,s[i]);
// }
printf("%lld\n",ans);
}
int main(){
scanf("%d%d",&n,&m);
get__mu(510000);
// printf("mu=%lld\n",prime[3]);
for(ll i=1;i<=n;i++)
scanf("%d",&su[i]);
for(ll i=1,a;i<=m;i++){
scanf("%d",&a);
ll x=su[a];
if(get_[a]==0)
work(x,1),get_[a]=1;
else work(x,-1),get_[a]=0;
}
}
NOIP模拟测试23「mine·water·gcd」的更多相关文章
- NOIP模拟测试19「count·dinner·chess」
反思: 我考得最炸的一次 怎么说呢?简单的两个题0分,稍难(我还不敢说难,肯定又有人喷我)42分 前10分钟看T1,不会,觉得不可做,完全不可做,把它跳了 最后10分钟看T1,发现一个有点用的性质,仍 ...
- NOIP模拟测试30「return·one·magic」
magic 题解 首先原式指数肯定会爆$long$ $long$ 首先根据欧拉定理我们可以将原式换成$N^{\sum\limits_{i=1}^{i<=N} [gcd(i,N)==1] C_{G ...
- NOIP模拟测试38「金·斯诺·赤」
金 辗转相减见祖宗 高精 #include<bits/stdc++.h> using namespace std; #define A 2000 #define P 1 #define N ...
- NOIP模拟测试28「阴阳·虎·山洞」
写这几个题解我觉得我就像在按照官方题解抄一样 阴阳 题解 将题目中给的阴阳看作黑色和白色 首先我们观察到最后生成图中某种颜色必须是竖着单调递增或竖着单调递减 类似这样 否则不满足这个条件 但合法染色方 ...
- NOIP模拟测试20「周·任·飞」
liu_runda出的题再次$\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%\%$ 任 题解 题目中为什么反复强调简单路径,没有环 没有环的图中点数-边数=联通块数 前缀和维护边 ...
- NOIP模拟测试16「Drink·blue·weed」
话说这次考试 Drink 非常棒的一道卡常练习题,适合练习卡常 真的很棒 前置卡常知识 1.char要比int快 char是最快的 输出putchar,输入getchar 在这个题快了7000豪 2. ...
- NOIP模拟测试9「随·单·题」
liu_runda出的题,先$\%\%\%\%\%\%\%\%\%\%\%$为敬 随 考试时没有Qj 然后甚至没做,甚至没交 我不知道我怎么想的 这个题挺难改 你需要用到 循环矩阵快速幂,矩阵快速幂优 ...
- NOIP模拟测试4「礼物·通讯·奇袭」
礼物. 首先见到期望一定要想dp,看到n的范围无脑想状压, 然后我就只想到这了. dp方程式还是比较好想的,但是我依然想不出来 略经思考 颓题解 依然不会,随便写了个式子 i状态中不含j $f[i ...
- NOIP模拟测试26「嚎叫响彻在贪婪的机房·主仆见证了 Hobo 的离别·征途堆积出友情的永恒」
题目比较神仙,注意是题目神仙 贪婪暗示贪心,堆积暗示堆优化$\%\%\%\%\%\%\%$ 两个乱搞$+$一个堆优化$dp$ 嚎叫响彻在贪婪的机房 题解 对于一个序列来说只要他们差的$gcd$不为$1 ...
随机推荐
- C#中的元组(Tuple)和结构体(struct)
在正常的函数调用中,一个函数只能返回一个类型的值,但在某些特殊情况下,我们可能需要一个方法返回多个类型的值,除了通过ref,out或者泛型集合可以实现这种需求外,今天,讲一下元组和结构体在这一方面的应 ...
- Python数模笔记-Sklearn(5)支持向量机
支持向量机(Support vector machine, SVM)是一种二分类模型,是按有监督学习方式对数据进行二元分类的广义线性分类器. 支持向量机经常应用于模式识别问题,如人像识别.文本分类.手 ...
- 关于MDI多文档程序的一些思考
MDI程序的框架 客户窗口是一个预定义的窗口类(MDICLIENT),它是框架窗口的子窗口同时也是各个子文档窗口的父窗口.框架窗口和各个子文档窗口都是自定义的窗口类. MDI程序中的一些要点 窗口中的 ...
- PE文件中的输入表
前言 PE文件中的输入表含有三个重要结构IID,IDT,IAT.PE文件为需要加载的DLL文件创建一个IID结构,一个DLL与一个IID对应.IDT是输入名称表,IAT输入地址表,在没有绑定输入的情况 ...
- redis分布式锁-可重入锁
redis分布式锁-可重入锁 上篇redis实现的分布式锁,有一个问题,它不可重入. 所谓不可重入锁,即若当前线程执行某个方法已经获取了该锁,那么在方法中尝试再次获取锁时,就会获取不到被阻塞. 同一个 ...
- ES6中的箭头函数的语法、指向、不定参数
箭头函数的语法 function fn1() { console.log(1); } let fn2 = () => { console.log(2); } fn1()//1 fn2()//2 ...
- ruby基础(二)
ruby语法基础 1.方法 方法时对象定义的与该对象相关的操作.在Ruby中,对象的所有的操作都被封装成 方法. 语法糖:语法糖是一种为了照顾一般人的习惯而产生的特殊语法. ruby中一切数据都是对象 ...
- [DB] 大数据集群安装
学习要点 体系架构.原理 多做练习.试验 装虚拟机 网络模式:仅主机模式 software selection:development tools, GUI network & host na ...
- Win10屏幕亮度不能调节,调节无效怎么办?
Win10屏幕亮度不能调节,调节无效怎么办? 听语音 浏览:1027 | 更新:2019-11-22 11:43 1 2 3 4 5 6 7 分步阅读 一些用户在使用win10系统之后,出现了电脑屏幕 ...
- Centos7.4 docker安装包下载以及离线安装
docker安装包下载地址:https://download.docker.com/linux/centos/7/x86_64/stable/Packages/ 需要下载一个selinux包:dock ...