形象理解这一场的 C
A.Seats \(\text{diff }20\)
对给定序列 \(S\) 找出 \(i\) 的个数,使得 \(S_{i}=0,S_{i+1}=1,S_{i+2}=0\)
#define int long long
string x;signed main(){int n;cin>>n;cin>>x;int ans=0;for(int i=0;i<=(int)x.length()-3;++i){if(x[i]=='#' and x[i+2]=='#' and x[i+1]=='.'){ans++;}}cout<<ans;
}
B.Traveling Takahashi Problem \(\text{diff }65\)
给定平面上 \(N\) 个点,求从原点出发按编号顺序走完 \(N\) 个点并走回原点的直线距离之和
#include<bits/stdc++.h>
using namespace std;
#define int long long
long double dist(int x1,int y1,int x2,int y2){return sqrtl((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
int n;
int nx=0,ny=0;
long double ans=0;
signed main(){cin>>n;for(int i=1;i<=n;++i){int x,y;cin>>x>>y;ans+=dist(nx,ny,x,y);nx=x;ny=y;}ans+=dist(nx,ny,0,0);printf("%.20Lf",ans);
}
C.Spiral Rotation \(\text{diff }972\)
给定一个 \(N\times N\) 的 \(01\) 矩阵,\(N\) 为偶数,对每个 \(i\in[1,\frac{N}{2}]\) 执行如下操作
- 选择所有 \(i\le x,y\le N-i+1\),将 \((y,N-x+1)\) 上的值替换为 \((x,y)\) 上的值
对于每一个 \(i\),对所有 \((x,y)\) 的替换同时发生,也即在每一个 \(i\) 内不存在覆盖的情况
输出最终的矩阵\(N\le 10^3\)
通过手摸可以发现,对单独的一个 \(i\) 来说,相当于四个一组发生旋转
点击查看
但是这里有一个 \(i\in[1,\frac{N}{2}]\),这其实相当于是分层的,先整体转,再转第二层和以内的,再转第三层和以内的,相当于第一层转一圈,第二层转两圈,以此类推
可以发现转四圈不会影响结果,因此第 \(k\) 层转 \(k\mod 4\) 圈就行了
#include<bits/stdc++.h>
using namespace std;
#define int long long
int n;
char mp[3001][3001];
bool vis[3001][3001];
signed main(){cin>>n;for(int i=1;i<=n;++i){for(int j=1;j<=n;++j){cin>>mp[i][j];}}for(int k=1;k<=n/2;++k){for(int i=k;i<=n-k+1;++i){if(!vis[i][k]){if(k%4==1){int x1=i,y1=k;int x2=y1,y2=n-x1+1;int x3=y2,y3=n-x2+1;int x4=y3,y4=n-x3+1;char tmp=mp[x4][y4];mp[x4][y4]=mp[x3][y3];mp[x3][y3]=mp[x2][y2];mp[x2][y2]=mp[x1][y1];mp[x1][y1]=tmp;vis[x1][y1]=true;vis[x2][y2]=true;vis[x3][y3]=true;vis[x4][y4]=true;}else if(k%4==2){int x1=i,y1=k;int x2=y1,y2=n-x1+1;int x3=y2,y3=n-x2+1;int x4=y3,y4=n-x3+1;char tmp=mp[x4][y4];mp[x4][y4]=mp[x3][y3];mp[x3][y3]=mp[x2][y2];mp[x2][y2]=mp[x1][y1];mp[x1][y1]=tmp;tmp=mp[x4][y4];mp[x4][y4]=mp[x3][y3];mp[x3][y3]=mp[x2][y2];mp[x2][y2]=mp[x1][y1];mp[x1][y1]=tmp;vis[x1][y1]=true;vis[x2][y2]=true;vis[x3][y3]=true;vis[x4][y4]=true;}else if(k%4==3){int x1=i,y1=k;int x2=y1,y2=n-x1+1;int x3=y2,y3=n-x2+1;int x4=y3,y4=n-x3+1;char tmp=mp[x4][y4];mp[x4][y4]=mp[x3][y3];mp[x3][y3]=mp[x2][y2];mp[x2][y2]=mp[x1][y1];mp[x1][y1]=tmp;tmp=mp[x4][y4];mp[x4][y4]=mp[x3][y3];mp[x3][y3]=mp[x2][y2];mp[x2][y2]=mp[x1][y1];mp[x1][y1]=tmp;tmp=mp[x4][y4];mp[x4][y4]=mp[x3][y3];mp[x3][y3]=mp[x2][y2];mp[x2][y2]=mp[x1][y1];mp[x1][y1]=tmp;vis[x1][y1]=true;vis[x2][y2]=true;vis[x3][y3]=true;vis[x4][y4]=true;}}}for(int j=k+1;j<=n-k;++j){if(!vis[k][j]){if(k%4==1){int x1=k,y1=j;int x2=y1,y2=n-x1+1;int x3=y2,y3=n-x2+1;int x4=y3,y4=n-x3+1;char tmp=mp[x4][y4];mp[x4][y4]=mp[x3][y3];mp[x3][y3]=mp[x2][y2];mp[x2][y2]=mp[x1][y1];mp[x1][y1]=tmp;vis[x1][y1]=true;vis[x2][y2]=true;vis[x3][y3]=true;vis[x4][y4]=true;}else if(k%4==2){int x1=k,y1=j;int x2=y1,y2=n-x1+1;int x3=y2,y3=n-x2+1;int x4=y3,y4=n-x3+1;char tmp=mp[x4][y4];mp[x4][y4]=mp[x3][y3];mp[x3][y3]=mp[x2][y2];mp[x2][y2]=mp[x1][y1];mp[x1][y1]=tmp;tmp=mp[x4][y4];mp[x4][y4]=mp[x3][y3];mp[x3][y3]=mp[x2][y2];mp[x2][y2]=mp[x1][y1];mp[x1][y1]=tmp;vis[x1][y1]=true;vis[x2][y2]=true;vis[x3][y3]=true;vis[x4][y4]=true;}else if(k%4==3){int x1=k,y1=j;int x2=y1,y2=n-x1+1;int x3=y2,y3=n-x2+1;int x4=y3,y4=n-x3+1;char tmp=mp[x4][y4];mp[x4][y4]=mp[x3][y3];mp[x3][y3]=mp[x2][y2];mp[x2][y2]=mp[x1][y1];mp[x1][y1]=tmp;tmp=mp[x4][y4];mp[x4][y4]=mp[x3][y3];mp[x3][y3]=mp[x2][y2];mp[x2][y2]=mp[x1][y1];mp[x1][y1]=tmp;tmp=mp[x4][y4];mp[x4][y4]=mp[x3][y3];mp[x3][y3]=mp[x2][y2];mp[x2][y2]=mp[x1][y1];mp[x1][y1]=tmp;vis[x1][y1]=true;vis[x2][y2]=true;vis[x3][y3]=true;vis[x4][y4]=true;}}}}for(int i=1;i<=n;++i){for(int j=1;j<=n;++j){cout<<mp[i][j];}cout<<'\n';}
}
D.ABA \(\text{diff }658\)
给定一个字符串 \(S\),求满足下列条件的 \((i,j,k)\) 对数
- \(i\lt j\lt k\)
- \(S_i,S_j,S_k\) 能构成回文
\(|S| \le 2\times 10^5\)
其实也就是找 \(S_i=S_k\) 的对数,然后中间填什么都行
注意到两个相等的数,如果位置分别为 \(i,j\),那么它们中间夹着的任何一个数都能造成贡献,也即贡献为 \(j-i-1\)
如果枚举到一个数 \(S_i\) 时,有 \(k\) 个数(\(S_{j1},S_{j2}\cdots\))与其相等,则对答案的贡献应该是 \((S_i-S_{j1}-1)+(S_i-S_{j2}-1)\cdots=k(S_i-1)-\sum S_{j}\)
因此维护前缀和即可
#include<bits/stdc++.h>
using namespace std;
#define int long long
int possum[27];
int poscnt[27];
string x;
int ans=0;
signed main(){cin>>x;for(int i=0;i<=(int)x.length()-1;++i){ans+=max(0ll,(i-1)*poscnt[x[i]-'A']-possum[x[i]-'A']);poscnt[x[i]-'A']++;possum[x[i]-'A']+=i;}cout<<ans;
}
E.3 Team Division \(\text{diff }1424\)
有三支队伍和 \(N\) 个人,每支队伍初始有一定人数,每个人有权值 \(B_i\),可以移动任意队伍中的人,问三支队伍的权值和相等所需的最小移动次数
\(N\le 100,\sum B_i\le 1500\)
首先注意到
- \(\sum B_i \not \equiv 0\pmod 3\) 无解
考虑换种思路做,清空所有队里的人,然后一个一个往里填,如果当前填的人本来就在这个队里,则对移动次数无贡献,否则有 \(1\) 的贡献
然后,注意到最终每个队的总和不超过 \(\frac{\sum B_i}{3}=500\),因此设计 \(f_{i,j,k}\) 表示放到第 \(i\) 个人,第一队有 \(j\) 的权值和,第二队有 \(k\) 的权值和的最少交换次数(第三队可以用 \(\sum B_i-j-k\) 算出来)
直接枚举每个人填数即可
最后从 \(f_{n,\frac{\sum}{3},\frac{\sum}{3},\frac{\sum}{3}}\) 统计答案
#include<bits/stdc++.h>
using namespace std;
#define int long long
int n;
int a[101],b[101];
int f[101][501][501];
int sum[5];
int tot;
signed main(){cin>>n;for(int i=1;i<=n;++i){cin>>a[i]>>b[i];sum[a[i]]+=b[i];tot+=b[i];}if(tot%3!=0){cout<<-1;return 0;}memset(f,0x3f,sizeof f);f[0][0][0]=0;for(int i=1;i<=n;++i){for(int j=0;j<=tot/3;++j){for(int k=0;k<=tot/3;++k){if(j>=b[i]) f[i][j][k]=min(f[i][j][k],f[i-1][j-b[i]][k]+(a[i]!=1));if(k>=b[i]) f[i][j][k]=min(f[i][j][k],f[i-1][j][k-b[i]]+(a[i]!=2));if(tot-j-k>=b[i]) f[i][j][k]=min(f[i][j][k],f[i-1][j][k]+(a[i]!=3));}}}if(f[n][tot/3][tot/3]>=0x3f3f3f3f3f3f3f3f){cout<<-1;return 0;}cout<<f[n][tot/3][tot/3]<<endl;
}
F.Road Blocked \(\text{diff }1546\)
无向图中实现如下操作
1 i
删除第 \(i\) 条边2 x y
查询 \(x,y\) 两点的最短路径长度\(N\le 300,Q\le 2\times 10^5\)
有删边,应该倒过来处理询问
然后考虑加边怎么处理变化
用 DP 的思路,加了边之后,只有经过这条边两个端点的路径才能有可能改变,因为 \(N\) 比较小,因此考虑将这条新加的边的两个端点当中转点,枚举所有点对的路径来更新
类似弗洛伊德(其实就是弗洛伊德)
#include<bits/stdc++.h>
using namespace std;
#define int long long
int f[301][301];
int n,m,q;
struct edge{int l,r,w;
}e[300*300+1];
struct opera{int op;int id;int x,y;
}op[200001];
bool vis[300*300+1];
vector<int>ans;
signed main(){ios::sync_with_stdio(false);cin>>n>>m>>q;for(int i=1;i<=m;++i){int x,y,z;cin>>x>>y>>z;e[i]={x,y,z};}for(int i=1;i<=q;++i){cin>>op[i].op;if(op[i].op==1){cin>>op[i].id;vis[op[i].id]=true;}else{cin>>op[i].x>>op[i].y;}}memset(f,0x3f,sizeof f);for(int i=1;i<=m;++i){if(!vis[i]){f[e[i].l][e[i].r]=min(f[e[i].l][e[i].r],e[i].w);f[e[i].r][e[i].l]=min(f[e[i].r][e[i].l],e[i].w);}}for(int k=1;k<=n;++k){for(int i=1;i<=n;++i){for(int j=1;j<=n;++j){f[i][j]=min(f[i][j],f[i][k]+f[k][j]);}}}for(int t=q;t>=1;--t){if(op[t].op==1){f[e[op[t].id].l][e[op[t].id].r]=min(f[e[op[t].id].l][e[op[t].id].r],e[op[t].id].w);f[e[op[t].id].r][e[op[t].id].l]=min(f[e[op[t].id].r][e[op[t].id].l],e[op[t].id].w);// cout<<"add "<<e[op[t].id].l<<" "<<e[op[t].id].r<<endl;for(int i=1;i<=n;++i){for(int j=1;j<=n;++j){f[i][j]=min(f[i][j],f[i][e[op[t].id].l]+f[e[op[t].id].l][j]);}}for(int i=1;i<=n;++i){for(int j=1;j<=n;++j){f[i][j]=min(f[i][j],f[i][e[op[t].id].r]+f[e[op[t].id].r][j]);}}}if(op[t].op==2){if(f[op[t].x][op[t].y]>=0x3f3f3f3f3f3f3f3f){ans.push_back(-1);}else ans.push_back(f[op[t].x][op[t].y]);}}std::reverse(ans.begin(),ans.end());for(int i:ans){cout<<i<<"\n";}
}
G.Road Blocked 2 \(\text{diff }2014\)
对给定无向图判断:
- 对所有 \(i\in[1,m]\),如果将第 \(i\) 条边删去,\(1\) 至 \(N\) 的最短路径是否发生变化
\(N,M\le 2\times 10^5\)
看的第二篇题解
判断这个其实就是判断一条边是不是会被所有 \(1\) 到 \(N\) 的最短路经过
我们可以根据 \(dis\) 判一下边 \((x,y)\) 是否在最短路上:
或者反着来也行,但是这样我们就能判断其是否在最短路上
然后我们需要判断一条边是否为所有最短路的必经边,判断一条边被几条最短路径经过是不容易的,但是判断 \(1\) 到一个节点有几条最短路径是容易的
因此我们可以这样判断:
- 如果 \(1\) 到 \(x\) 有 \(c_x\) 条最短路径,到 \(N\) 有 \(c_N\) 条最短路径,\(N\) 到 \(y\) 有 \(c_y\) 条最短路径,则如果 \(c_x\times c_y=c_N\),那么就说明边 \((x,y)\) 是最短路必经边
这么做有一个好处就是分别拿 \(1,N\) 跑两边最短路就能统计答案
#include <bits/stdc++.h>
using namespace std;
#define int long long
int n,m;
const int p=1e9+7;
struct __edge{int to,w;
};
vector<__edge>e[200001];
struct ___edge{int from,to,w;
};
___edge edge[200001];
int dis[2][200001];
bool vis[2][200001];
int cnt[2][200001];
struct node{int id,dis;bool operator <(const node &A)const{return dis>A.dis;}
};
priority_queue<node>q;
void dij(int s,int id){memset(dis[id],0x3f,sizeof dis[id]);memset(vis[id],0,sizeof vis[id]);memset(cnt[id],0,sizeof cnt[id]);dis[id][s]=0;cnt[id][s]=1;q.push({s,dis[id][s]});while(!q.empty()){node u=q.top();q.pop();if(dis[id][u.id]<u.dis) continue;if(vis[id][u.id]) continue;vis[id][u.id]=true;for(__edge i:e[u.id]){if(dis[id][i.to]>dis[id][u.id]+i.w){dis[id][i.to]=dis[id][u.id]+i.w;q.push({i.to,dis[id][i.to]});cnt[id][i.to]=cnt[id][u.id];}else if(dis[id][i.to]==dis[id][u.id]+i.w){cnt[id][i.to]=(cnt[id][i.to]+cnt[id][u.id])%p;}}}
}
signed main(){cin>>n>>m;for (int i=1;i<=m;i++){int u,v,w;cin>>u>>v>>w;e[u].push_back({v,w});e[v].push_back({u,w});edge[i]={u,v,w};}dij(1,0);dij(n,1);for(int i=1;i<=m;++i){int u=edge[i].from,v=edge[i].to,w=edge[i].w;if((dis[0][u]+w+dis[1][v]==dis[0][n] and cnt[0][u]*cnt[1][v]%p==cnt[0][n])or (dis[0][v]+w+dis[1][u]==dis[0][n] and cnt[0][v]*cnt[1][u]%p==cnt[0][n])){cout<<"Yes\n";}else{cout<<"No\n";}}
}