前言

240分
100+80+20+40
T3少取了一个模结果全挂掉了（好不容易推出来了…）
T2也因为各种奇怪的错误挂了分
qwq
吸取教训吧

考场

今天先看题
T1第一眼看错了题意觉得水的不行
T3YBT原题且很水
T2T4乍看不太可做

因为错误审题，决定先写自以为很水的T1
写完样例不对，才发现自己审题完全审错了
那个翻折的关键性质又没有看出来（都看出来了吗，我感觉并不显然啊qwq）
分析了一会性质感觉打表或许可行
但可能是道搬砖题
此时大概8:40
于是转T3稳一稳

此时心态有些不稳了
（我甚至觉得T3的3e7会炸）
把T3超级暴力的《正解》floyd敲完还不太放心
但是也没有什么太好的方法
转回T1死磕
9:30

T1走上了慢慢打表之路
打了30min终于把恶心的转移数组敲完了
此时已经感觉药丸
但是编译结果有点惊喜
一遍过了样例？？？
又测了几个位置也对了
说实话那个搬砖代码一个bug没有真的挺神奇的
我本来已经做好了和它死磕的准备
结果就这么过了？
重新燃起希望
10:20

来到T2
很不错的时我很快发现了它的关键性质
利用一半为分界点递归转移
然后后面再跟个二分
check的地方递归求个数
很愉快的切掉了这道挺难的题
有一说一这道应该是我这几次模拟做的最好的一道题
确实有逻辑分析在里面
遗憾的是我切掉的时候也有点激动
手一抖最后的一步乘法没有取模
…
11:00

此时我以为我已经切了3题
心态就平和了起来
推了一会T4没有任何思路
就决定写个好一点的部分分下班
20分还是很好写的
40分也不太难
（又开始暴力搬砖）
11:30

这次时间有些紧凑
T1和T3浪费了太多时间
我没有很多的时间检查了
（说不定再看看T2的沙雕错误就瞅出来了呢qwq）
最后就带着那个bug走出了考场
qwq

复盘

T1 road

正解简单的离谱
关键是回溯时对原图进行一个神奇的反转操作
就很easy了
但为了纪念还是附上打表代码

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int N=1e6+100;
int n;
int a,b;
int mi[32];
int shunpl[5][5]={{},{0,1,1,1,3},{0,4,2,2,2},{0,3,1,3,3},{0,4,4,2,4}};
int shundir[5][5]={{},{0,2,1,1,2},{0,2,2,1,1},{0,1,2,2,1},{0,1,1,2,2}};
int shunid[5][5]={{},{0,1,2,3,4},{0,4,1,2,3},{0,3,4,1,2},{0,2,3,4,1}};
int nipl[5][5]={{},{0,1,3,1,1},{0,2,2,4,2},{0,3,3,3,1},{0,2,4,4,4}};
int nidir[5][5]={{},{0,1,1,2,2},{0,2,1,1,2},{0,2,2,1,1},{0,1,2,2,1}};
int niid[5][5]={{},{0,1,4,3,2},{0,2,1,4,3},{0,3,2,1,4},{0,4,3,2,1}};
int dx[5]={0,0,0,1,1},dy[5]={0,0,1,1,0};
/*
1:shun
2:ni
*/
int x2,y2,x3,y3;
void find(int k,int pl,int dir,int id,int x,int y){int ed=id+mi[2*k]-1;if((a<id||a>ed)&&(b<id||b>ed)) return;if(k==0){if(id==a){x2=x;y2=y;}if(id==b){x3=x;y3=y;}return;}if(dir==1){//shunint len=mi[k-1],tot=mi[(k-1)*2];for(int i=1;i<=4;i++){find(k-1,shunpl[pl][i],shundir[pl][i],id+(shunid[pl][i]-1)*tot,x+dx[i]*len,y+dy[i]*len);}return;}else{int len=mi[k-1],tot=mi[(k-1)*2];for(int i=1;i<=4;i++){find(k-1,nipl[pl][i],nidir[pl][i],id+(niid[pl][i]-1)*tot,x+dx[i]*len,y+dy[i]*len);}return;}
}
int main(){//  freopen("road.in","r",stdin);
//  freopen("road.out","w",stdout);int t;scanf("%d",&t);mi[0]=1;for(int i=1;i<=30;i++) mi[i]=mi[i-1]<<1;while(t--){scanf("%d%d%d",&n,&a,&b);find(n,1,1,1,1,1);//printf("a=(%d,%d) b=(%d,%d)\n",x2,y2,x3,y3);printf("%.0lf\n",10.0*sqrt((x3-x2)*(x3-x2)+(y3-y2)*(y3-y2)));}return 0;
}
/*
3
1 1 2
2 16 1
3 4 33
*/

T2 shop

本题std的二分很神奇
类似于一个二分相对大小的最大值
然后硬限制一个范围
就过了…
觉得还是我的代码更好理解一写qwq（所有人对自己的代码都是这么想的）
不取模，见祖宗

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int N=105;
const int mod=1e9+7;
int n,k;
int solve(int x){if(x==0) return 0;return x+solve(x/2);
}
int find(int l,int r){if(l>r) return 0;if(r==0) return 0;return r-l+1 + find((l+1)/2,r/2);
}
ll ksm(ll x,int k){ll ans=1,res=x;while(k){if(k&1) ans=ans*res%mod;res=res*res%mod;k>>=1;}return ans;
}
void work(int k,int n){int tot=solve(n/2);if(k<=tot){work(k,n/2);return;}k-=tot;
//  printf("n=%d tot=%d k=%d\n",n,tot,k);int o=k/n,p=k%n;if(p==0){printf("%lld\n",n*ksm(2,o-1)%mod);return;}int st=n/2+1,ed=n;while(st<ed){int mid=st+ed>>1;
//      printf("  l=%d r=%d find=%d\n",n/2+1,mid,find(n/2+1,mid));if(find(n/2+1,mid)>=p) ed=mid;else st=mid+1;}printf("%lld\n",st*ksm(2,o)%mod);
}
int main(){//  freopen("shop.in","r",stdin);
//  freopen("shop.out","w",stdout);int t;scanf("%d",&t);while(t--){scanf("%d%d",&k,&n);work(k,n);}return 0;
}
/*
5
3 1
3 2
5 3
1000000000 1
987654321 8765432101
32 10
*/

T3 run

YBT原题且很水
这也挂分就离谱
有一个是没有特判m=1的情况
还有一个是因为int的范围k只枚举到30而不是31！
头疼

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int N=105;
int n,m;
bool f[33][N][N];
int dis[N][N];
int main(){//  freopen("run.in","r",stdin);
//  freopen("run.out","w",stdout);int t;scanf("%d",&t);while(t--){memset(f,0,sizeof(f));memset(dis,0x3f,sizeof(dis));
//      if(n==1){//          printf("0\n");
//          continue;
//      }scanf("%d%d",&n,&m);for(int i=1;i<=m;i++){int x,y;scanf("%d%d",&x,&y);f[0][x][y]=1;}for(int k=1;k<=30;k++){for(int i=1;i<=n;i++){for(int j=1;j<=n;j++){//if(i==j) continue;for(int p=1;p<=n;p++){//if(p==i||p==j) continue;if(f[k-1][i][p]&&f[k-1][p][j]){f[k][i][j]=1;
//                          printf("  p=%d\n",p);break;}}
//                  printf("k=%d i=%d j=%d f=%d\n",k,i,j,f[k][i][j]);}}}for(int i=1;i<=n;i++){for(int j=1;j<=n;j++){//if(i==j) continue;for(int k=0;k<=30;k++){if(f[k][i][j]){//                      printf("ok:k=%d i=%d j=%d\n",k,i,j);dis[i][j]=1;break;}}}dis[i][i]=0;}for(int k=1;k<=n;k++){for(int i=1;i<=n;i++){for(int j=1;j<=n;j++){dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);}}}printf("%d\n",dis[1][n]);}return 0;
}
/*
1
5 5
1 4
4 1
1 2
2 3
3 5
*/

T4 stairs

这题考场确实不太可做
考后写起来发现反而不难了
状压求出转移矩阵再快速幂加速即可
一个重要的技巧是把0规定为有挡板
这样不同高度之间可以无缝连接

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int N=105;
const int mod=1e9+7;
int n,k;
int w[8];
struct matrix{int x,y;ll a[130][130];void clear(){for(int i=0;i<=x;i++)for(int j=0;j<=y;j++) a[i][j]=0;}
};
int mi[8];
void print(matrix o){for(int i=0;i<=min(3,o.x);i++){for(int j=0;j<=3;j++) printf("%d ",o.a[i][j]);printf("\n");}return;
}
matrix cheng(matrix a,matrix b){matrix o;//printf("Mul:\n") ;//print(a);print(b);o.x=a.x;o.y=b.y;o.clear();for(int i=0;i<=o.x;i++){for(int j=0;j<=o.y;j++){for(int k=0;k<=a.y;k++){o.a[i][j]=(o.a[i][j]+a.a[i][k]*b.a[k][j])%mod;//if(a.a[i][k]*b.a[k][j])   printf("  i=%d j=%d k=%d %d*%d a=%d\n",i,j,k,a.a[i][k],b.a[k][j],o.a[i][j]);}}}return o;
}
int dp[130][2];//dp[i][op]op:下面的边
matrix trans[8];
void init(){mi[0]=1;for(int i=1;i<=7;i++) mi[i]=mi[i-1]<<1;for(int k=1;k<=7;k++){trans[k].x=trans[k].y=mi[7]-1;for(int i=0;i<mi[k];i++){for(int j=0;j<mi[k];j++){dp[0][1]=1;dp[0][0]=0;for(int p=1;p<=k-1;p++){if((i|j)&mi[p-1]){dp[p][1]=dp[p-1][0]+dp[p-1][1];dp[p][0]=dp[p-1][1]+dp[p-1][0];}else{dp[p][1]=dp[p-1][0];dp[p][0]=dp[p-1][1]+dp[p-1][0];}dp[p][0]%=mod;dp[p][1]%=mod;}if((i|j)&mi[k-1]) trans[k].a[i][j]=dp[k-1][0]+dp[k-1][1];else trans[k].a[i][j]=dp[k-1][0];//printf("k=%d i=%d j=%d trans=%d\n",k,i,j,trans[k].a[i][j]);}}}
}
matrix ans;
int main(){for(int i=1;i<=7;i++) scanf("%d",&w[i]);init();ans.x=0;ans.y=mi[7]-1;ans.a[0][0]=1;for(int i=1;i<=ans.y;i++) ans.a[0][i]=0;for(int k=1;k<=7;k++){//printf("k=%d trans:\n",k);//print(trans[k]);while(w[k]){//printf("k=%d\n",w[k]);if(w[k]&1){ans=cheng(ans,trans[k]);//printf("ans:\n");//for(int i=0;i<=4;i++) printf("%d ",ans.a[0][i]);//printf("\n");}trans[k]=cheng(trans[k],trans[k]);//printf("trans:\n");//print(trans[k]);w[k]>>=1;}}printf("%lld\n",ans.a[0][0]);return 0;
}
/*
5
3 1
3 2
5 3
1000000000 1
987654321 8765432101
32 10
*/

总结

争取吃一堑长一智吧
以后检查代码注意的东西又多了一个：取模是否彻底！
明天：字符串加油！awa

8.13模拟：分治二分倍增快速幂相关推荐

BZOJ.4180.字符串计数(后缀自动机二分矩阵快速幂/倍增Floyd)
题目链接先考虑假设S确定,使构造S操作次数最小的方案应是:对T建SAM,S在SAM上匹配,如果有S的转移就转移,否则操作数++,回到根节点继续匹配S.即每次操作一定是一次极大匹配. 简单证明:假设 ...
【PTA】7-2 国王的奖励——分数取模、分治思想、快速幂、int64的乘法模运算【C/C++】
文章目录 1 问题 1.1 题面描述 1.2 输入描述 1.3 输出描述 1.4 样例描述 1.5 样例解释 2 分析 2.1 数学抽象 2.2 解决方法 2.2.1 等比求和 2.2.1.1 求和公 ...
LeetCode 1498. 满足条件的子序列数目（排序+二分查找+快速幂）
1. 题目给你一个整数数组 nums 和一个整数 target . 请你统计并返回 nums 中能满足其最小元素与最大元素的和小于或等于 target 的非空子序列的数目. 由于答案可能很大 ...
[二分] WYF快速幂魔改版
题目出处 qbxt夏北京二期,提高组精英版,WYF出题. 题目给定一个1~n的排列a,再给出m组询问.每次询问,你需要计算a[a[a[a-[x]]-],总共有t组括号. n,m<=100000 ...
[NOIP2003普及组]麦森数（快速幂+高精度）
[NOIP2003普及组]麦森数(快速幂+高精度) Description 形如2^P-1的素数称为麦森数,这时P一定也是个素数.但反过来不一定,即如果P是个素数,2^P-1不一定也是素数.到1998 ...
NBUT 1225 NEW RDSP MODE I（找规律）（快速幂）
NEW RDSP MODE I 问题描述 Little A has became fascinated with the game Dota recently, but he is not a goo ...
快速幂、快速乘、矩阵快速幂
文章目录一快速幂模板模板例题二快速乘模板三矩阵快速幂例题模板题应用题1 应用题2 一快速幂目的:当我们在求f(x) = a ^ x % mod时,f(x)的结果会呈指数型增 ...
2018.11.08 NOIP模拟景点（倍增+矩阵快速幂优化dp）
传送门首先按照题意构造出转移矩阵. 然后可以矩阵快速幂求出答案. 但是直接做是O(n3qlogm)O(n^3qlogm)O(n3qlogm)的会TTT掉. 观察要求的东西发现我们只关系一行的答案. ...
分治算法——快速幂（平方求幂）
分治算法--快速幂(平方求幂) 什么是快速幂呢,我们先来看维基百科对快速幂的解释: 在数学和程序设计中,平方求幂(英语:exponentiating by squaring)或快速幂是快速计算一个数( ...

8.13模拟：分治二分倍增快速幂

文章目录

前言

考场

复盘

T1 road

T2 shop

T3 run

T4 stairs

总结

8.13模拟：分治二分倍增快速幂相关推荐

最新文章

热门文章