哦 今天一时兴起想刷提高课 不知道是提高课的题太水了 还是讲的很细
大概两个多小时就能刷十多道吧
P1967 [NOIP2013 提高组] 货车运输
这道题思路挺简单的吧 不过我好像是dfs 那里写错了一点调了蛮久的
这个大概可以一眼看出最大生成树
可是我最开始想的就是dfs来找最大载重 时间复杂度大概是O（q*n）1e8感觉可以通过的啊
可是题解里面都是lca 那没办法 那我们就用倍增lca吧 时间复杂度是O（qlogn+nlog）
这里fa就是lca里的fa
W和fa的定义差不多 该点距2^j的点的最大限重
注意我们更新时就因该把k放在外层 因为是倍增嘛
dfs的时候记得把根节点初始成inf 其实不初始也能ac
最后还要注意的就是dfs时不要和我一样
fa[j][0]的爹是x;
depth[j]=depth[x]+1;
#include <bits/stdc++.h>
using namespace std;
const int N = 5e4+10;
const int M = 5e4+10;
const int mod = 998244353;
#define int long long
#define endl '\n'
#define Endl '\n'
#define _ 0
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}
int n,m,p[N],h[N],e[M],ne[M],w[M],idx,W[N][21],fa[N][21],depth[N];
bool vis[N];
struct Edge{
int a,b,w;
bool operator < (const Edge&W)const{
return w>W.w;
}
}edge[M];
void add(int a,int b,int c){
e[idx]=b;
w[idx]=c;
ne[idx]=h[a];
h[a]=idx++;
}
int find(int x){
if(p[x]!=x)p[x]=find(p[x]);
return p[x];
}
void kruskal(){
sort(edge,edge+m);
for(int i=1;i<=n;i++)p[i]=i;
for(int i=0;i<m;i++){
int a=edge[i].a,b=edge[i].b,w=edge[i].w;
if(find(a)!=find(b)) {
p[find(a)]=find(b);
add(a,b,w),add(b,a,w);
}
}
}
void dfs(int x){
vis[x]=1;
for(int i=h[x];~i;i=ne[i]){
int j=e[i];
if(!vis[j]){
fa[j][0]=x;
depth[j]=depth[x]+1;
W[j][0]=w[i];
for(int k=1;k<=20;k++)fa[j][k]=fa[fa[j][k-1]][k-1];
dfs(j);
}
}
}
int lca(int a,int b){
if(depth[a]<depth[b])swap(a,b);
int ans=inf;
for(int k=20;k>=0;k--){
if(depth[fa[a][k]]>=depth[b]){
ans=min(W[a][k],ans);//
a=fa[a][k];
}
}
if(a==b)return ans;
for(int k=20;k>=0;k--){
if(fa[a][k]!=fa[b][k]){
ans=min(ans,min(W[a][k],W[b][k]));
a=fa[a][k];
b=fa[b][k];
}
}
ans=min(ans,min(W[a][0],W[b][0]));//
return ans;
}
signed main(){
cin>>n>>m;
memset(h,-1,sizeof h);
for(int i=0;i<m;i++){
int a,b,c;cin>>a>>b>>c;
edge[i]={a,b,c};
}
kruskal();
for(int i=1;i<=n;i++){
if(!vis[i]){
depth[i]=1;
dfs(i);
W[i][0]=inf;
}
}
for(int k=1;k<=20;k++){
for(int i=n;i>=1;i--){
W[i][k]=min(W[i][k-1],W[fa[i][k-1]][k-1]);
}
}
int q;cin>>q;
while(q--){
int a,b;cin>>a>>b;
if(find(a)!=find(b))cout<<-1<<endl;
else cout<<lca(a,b)<<endl;
}
return ~~(0^_^0);
}下面基本都很简单就随便贴一贴把
AcWing 1015. 摘花生
#include <bits/stdc++.h>
using namespace std;
const int N = 110;
const int M = 5e4+10;
const int mod = 998244353;
#define int long long
#define endl '\n'
#define Endl '\n'
#define _ 0
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}
int f[N][N],a[N][N];
signed main(){
int t;cin>>t;
while(t--){
int n,m;cin>>n>>m;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
cin>>a[i][j];
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
f[i][j]=max(f[i-1][j],f[i][j-1])+a[i][j];
cout<<f[n][m]<<endl;
}
return ~~(0^_^0);
}AcWing 1018. 最低通行费
#include <bits/stdc++.h>
using namespace std;
const int N = 110;
const int M = 5e4+10;
const int mod = 998244353;
#define int long long
#define endl '\n'
#define Endl '\n'
#define _ 0
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}
int f[N][N],a[N][N];
signed main(){
int n;cin>>n;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
cin>>a[i][j];
memset(f,0x3f3f,sizeof f);
f[1][0]=f[0][1]=0;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
f[i][j]=min(f[i-1][j],f[i][j-1])+a[i][j];
cout<<f[n][n]<<endl;
return ~~(0^_^0);
}AcWing 1027. 方格取数
这道题还有点难度 但也不是特别难
主要是以前没见过 就是一个四维的吧 然后他简化成了三维
#include <bits/stdc++.h>
using namespace std;
const int N = 15;
const int M = 5e4+10;
const int mod = 998244353;
#define int long long
#define endl '\n'
#define Endl '\n'
#define _ 0
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}
int f[N<<1][N][N],a[N][N];
signed main(){
fast
int n;cin>>n;
int x,y,w;
while(cin>>x>>y>>w,x||y||w)a[x][y]=w;
for(int k=2;k<=n*2;k++){
for(int i1=1;i1<=n;i1++){
for(int i2=1;i2<=n;i2++){
int j1=k-i1,j2=k-i2;
if(j1<1||j1>n||j2<1||j2>n)continue;
if(i1==i2){
f[k][i1][i2]=max(max(max(f[k-1][i1][i2],f[k-1][i1-1][i2]),f[k-1][i1][i2-1]),f[k-1][i1-1][i2-1])+a[i1][k-i1];
}else{
f[k][i1][i2]=max(max(max(f[k-1][i1][i2],f[k-1][i1-1][i2]),f[k-1][i1][i2-1]),f[k-1][i1-1][i2-1])+a[i1][k-i1]+a[i2][k-i2];
}
}
}
}
cout<<f[2*n][n][n]<<endl;
return ~~(0^_^0);
}AcWing 275. 传纸条
和上面一题一模一样虽然这个多了不能相交的条件
这个随便想一下就知道吧
相交肯定不是最优解
#include <bits/stdc++.h>
using namespace std;
const int N = 55;
const int M = 5e4+10;
const int mod = 998244353;
#define int long long
#define endl '\n'
#define Endl '\n'
#define _ 0
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}
int f[N<<1][N][N],a[N][N];
signed main(){
fast
int n,m;cin>>n>>m;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
cin>>a[i][j];
for(int k=2;k<=n+m;k++){
for(int i1=1;i1<=n;i1++){
for(int i2=1;i2<=n;i2++){
int j1=k-i1,j2=k-i2;
if(j1<1||j1>m||j2<1||j2>m)continue;
if(i1==i2)f[k][i1][i2]=max(max(max(f[k-1][i1][i2],f[k-1][i1-1][i2]),f[k-1][i1][i2-1]),f[k-1][i1-1][i2-1])+a[i1][k-i1];
else f[k][i1][i2]=max(max(max(f[k-1][i1][i2],f[k-1][i1-1][i2]),f[k-1][i1][i2-1]),f[k-1][i1-1][i2-1])+a[i1][k-i1]+a[i2][k-i2];
}
}
}
cout<<f[n+m][n][n]<<endl;
return ~~(0^_^0);
}AcWing 1017. 怪盗基德的滑翔翼
#include <bits/stdc++.h>
using namespace std;
const int N = 110;
const int M = 5e4+10;
const int mod = 998244353;
#define int long long
#define endl '\n'
#define Endl '\n'
#define _ 0
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}
int f[N];
signed main(){
fast
int t;cin>>t;
while(t--){
int n;cin>>n;
int a[N];
for(int i=1;i<=n;i++)cin>>a[i];
for(int i=1;i<=n;i++)f[i]=1;
int ans=-2e9;
for(int i=1;i<=n;i++){
for(int j=1;j<i;j++){
if(a[i]>a[j])f[i]=max(f[i],f[j]+1);
ans=max(ans,f[i]);
}
}
for(int i=1;i<=n;i++)f[i]=1;
for(int i=n;i>=1;i--){
for(int j=i+1;j<=n;j++){
if(a[i]>a[j])f[i]=max(f[i],f[j]+1);
ans=max(ans,f[i]);
}
}
cout<<ans<<Endl;
}
return ~~(0^_^0);
}AcWing 1014. 登山
#include <bits/stdc++.h>
using namespace std;
const int N = 1010;
const int M = 5e4+10;
const int mod = 998244353;
#define int long long
#define endl '\n'
#define Endl '\n'
#define _ 0
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}
int n,f1[N],f[N];
signed main(){
fast
int n;cin>>n;
int a[N];
for(int i=1;i<=n;i++)cin>>a[i];
for(int i=1;i<=n;i++)f[i]=f1[i]=1;
for(int i=1;i<=n;i++)
for(int j=1;j<i;j++)
if(a[i]>a[j])f[i]=max(f[i],f[j]+1);
for(int i=n;i>=1;i--)
for(int j=i+1;j<=n;j++)
if(a[i]>a[j])f1[i]=max(f1[i],f1[j]+1);
int ans=-2e9;
for(int i=1;i<=n;i++)ans=max(ans,f[i]+f1[i]);
cout<<ans-1<<endl;
return ~~(0^_^0);
}AcWing 482. 合唱队形
#include <bits/stdc++.h>
using namespace std;
const int N = 1010;
const int M = 5e4+10;
const int mod = 998244353;
#define int long long
#define endl '\n'
#define Endl '\n'
#define _ 0
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}
int n,f1[N],f[N];
signed main(){
fast
int n;cin>>n;
int a[N];
for(int i=1;i<=n;i++)cin>>a[i];
for(int i=1;i<=n;i++)f[i]=f1[i]=1;
for(int i=1;i<=n;i++)
for(int j=1;j<i;j++)
if(a[i]>a[j])f[i]=max(f[i],f[j]+1);
for(int i=n;i>=1;i--)
for(int j=i+1;j<=n;j++)
if(a[i]>a[j])f1[i]=max(f1[i],f1[j]+1);
int ans=-2e9;
for(int i=1;i<=n;i++)ans=max(ans,f[i]+f1[i]);
cout<<n-ans+1<<endl;
return ~~(0^_^0);
}AcWing 1012. 友好城市
#include <bits/stdc++.h>
using namespace std;
const int N = 5010;
const int M = 5e4+10;
const int mod = 998244353;
#define int long long
#define endl '\n'
#define Endl '\n'
#define _ 0
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}
pair<int,int>p[N];
int f[N];
signed main(){
fast
int n;cin>>n;
for(int i=1;i<=n;i++)cin>>p[i].first>>p[i].second;
sort(p+1,p+n+1);
int a[N];
for(int i=1;i<=n;i++)a[i]=p[i].second,f[i]=1;
int ans=0;
for(int i=1;i<=n;i++){
for(int j=1;j<i;j++){
if(a[i]>a[j]){
f[i]=max(f[i],f[j]+1);
ans=max(ans,f[i]);
}
}
}
cout<<ans<<Endl;
return ~~(0^_^0);
}AcWing 1016. 最大上升子序列和
#include <bits/stdc++.h>
using namespace std;
const int N = 1010;
const int M = 5e4+10;
const int mod = 998244353;
#define int long long
#define endl '\n'
#define Endl '\n'
#define _ 0
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}
int f[N];
signed main(){
fast
int n;cin>>n;
int a[N];
int ans=-2e9;
for(int i=1;i<=n;i++)cin>>a[i],f[i]=a[i],ans=max(ans,a[i]);
for(int i=1;i<=n;i++){
for(int j=1;j<i;j++){
if(a[i]>a[j]){
f[i]=max(f[i],f[j]+a[i]);
ans=max(ans,f[i]);
}
}
}
cout<<ans<<endl;
return ~~(0^_^0);
}AcWing 1010. 拦截导弹
这个以前我做的时候还是绿题来着 刚刚去看了变黄了
确实该变变了 不过这种时候确实发现自己成长挺多的
#include <bits/stdc++.h>
using namespace std;
const int N = 1010;
const int M = 5e4+10;
const int mod = 998244353;
#define int long long
#define endl '\n'
#define Endl '\n'
#define _ 0
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}
int a[N],f[N],f1[N];
signed main(){
fast
int n=1;
while(cin>>a[n])n++;
n--;
int ans=-2e9;
for(int i=1;i<=n;i++)f1[i]=f[i]=1;
for(int i=1;i<=n;i++){
for(int j=1;j<i;j++){
if(a[j]>=a[i]){
f[i]=max(f[i],f[j]+1);
ans=max(ans,f[i]);
}
}
}
if(ans!=-2e9)cout<<ans<<endl;
else cout<<1<<endl;
ans=-2e9;
for(int i=1;i<=n;i++){
for(int j=1;j<i;j++){
if(a[i]>a[j]){
f1[i]=max(f1[i],f1[j]+1);
ans=max(ans,f1[i]);
}
}
}
if(ans!=-2e9)cout<<ans<<endl;
else cout<<1<<endl;
return ~~(0^_^0);
}AcWing 187. 导弹防御系统
这道题可以爆搜 可以迭代加深 因为答案很浅
然后我们用了个 类似以前的log级别的LIS模型
f[i]表示第i个数的最小值？
差不多吧50范围加个基础剪枝就过了
#include <bits/stdc++.h>
using namespace std;
const int N = 55;
const int M = 5e4+10;
const int mod = 998244353;
#define int long long
#define endl '\n'
#define Endl '\n'
#define _ 0
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}
int ans,down[N],up[N],n,a[N];
void dfs(int u,int s,int d){
if(s+d>=ans)return;
if(u==n+1){
ans=s+d;
return;
}
int k=0;
while(k<s&&up[k]>=a[u])k++;
int t=up[k];
up[k]=a[u];
if(k==s)dfs(u+1,s+1,d);
else dfs(u+1,s,d);
up[k]=t;
k=0;
while(k<d&&down[k]<=a[u])k++;
t=down[k];
down[k]=a[u];
if(k==d)dfs(u+1,s,d+1);
else dfs(u+1,s,d);
down[k]=t;
}
signed main(){
fast
while(cin>>n,n){
for(int i=1;i<=n;i++)cin>>a[i];
ans=2e9;
dfs(1,0,0);
cout<<ans<<endl;
}
return ~~(0^_^0);
}
扫一扫
238
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
