2022.9.27日调试笔记

前言

前几天调了几个很离谱的bug,在此记录一下坑点。

1:P1016 [NOIP1999 提高组] 旅行家的预算

本来以为只是一道很简单的贪心,但是发现代码确实很不好写。

下面给出我的代码:

#include<iostream>
#include<cstdio>
using namespace std;
double a[10],b[10],c,k,oil,ans;
int n;
int main() {
	cin>>a[0]>>c>>k>>b[0]>>n;
	swap(a[0],a[n+1]);
	for(int i=1;i<=n;i++)
		cin>>a[i]>>b[i];
	for(int i=1;i<=n+1;i++)
		if(a[i]-a[i-1]>c*k) {
			cout<<"No Solution"<<endl;
			return 0;
		}
	for(int i=0,t;i<=n;i=t) {
		for(int j=i+1;a[j]-a[i]<=c*k&&j<=n+1;j++) 
			if(b[j]<b[i]) {
				ans+=b[i]*((a[j]-a[i])/k-min(oil,(a[j]-a[i])/k));
				oil-=min(oil,(a[j]-a[i]/k));
				t=j;
				goto flag;
			}
		{
		int mint=i+1;
		for(int j=i+2;a[j]-a[i]<=c*k&&j<=n+1;j++) 
			if(b[j]<=b[mint]) 
				mint=j;
		ans+=(c-oil)*b[i];
		oil=c-(a[mint]-a[i])/k;
		t=mint;
		}
		flag:;
	}
	printf("%.2lf",ans);
	return 0;
}

主要思路就是模拟+贪心,贪心策略如下:
在加满油可以抵达的范围内,如果存在油价比所处加油站更低的加油站,那么行驶到此加油站,所需加油的数量只需要刚好能够抵达此加油站即可。
如果不存在这种加油站,那么就加满油,行驶到后面能抵达范围内油价相对最低的加油站。加满油是为了到了后面那个加油站可以少加一点,省钱。

题解区也有很多不同的做法,但这些做法正确的原因都是遵循了以下原则:
无论什么时候行驶,都尽量使用最便宜的油。

代码中有几句不容易理解:

对于第一种情况:

ans+=b[i]*((a[j]-a[i])/k-min(oil,(a[j]-a[i])/k));
oil-=min(oil,(a[j]-a[i]/k));

如果存在油价比所处加油站更低的加油站,那么应该尽量少加油,只需要能够行驶到此加油站即可。
这两句话的语义是:
第一句话,需要花费的钱数=原来已经花的钱数+加油站的油价*行驶所需新购买的油量。
其中,新购买的油量=行驶所需油量-行驶所需以前购买的油的数量。显然应该先烧之前购买的油,这样新购买的油量会少一些,比较省。
其中,行驶所需的以前购买的油量=min(以前购买的油量,行驶所需油量)。

第二句话,抵达下一个加油站之后的油量=现有油量-行驶所需油量+新购买的油量=现有油量-(行驶所需的以前购买的油量+新购买的油量)+新购买的油量=现有油量-行驶所需的以前购买的油量。

对于第二种情况:

ans+=(c-oil)*b[i];
oil=c-(a[mint]-a[i])/k;

这两句话的语义是:
第一句话,需要花的钱数=原来已经花的钱数+加油站的油价*加满油箱所需要的购买的油量。
其中,加满油箱所需要购买的油量=油箱容量-现有油量。

第二句话,抵达指定加油站之后的油量=油箱容量(因为油箱要先加满)-行驶到此加油站所消耗的油量。

其他类型代码调试:

机房同学给我一份75分代码:

#include<bits/stdc++.h>
using namespace std;
double d[10],jg[10],cha[10],as,ans;
int f; 
double maxx(double a,double b){
	return a>b?a:b;
}
int main(){
	double sy=0,m,v,w,p;
	int n;
	ios::sync_with_stdio(0);
	cin.tie(0);
	cout.tie(0);
	cin>>m>>v>>w>>p>>n;
	as=v*w;
	for(int i=1;i<=n;i++){
		cin>>d[i]>>jg[i];
		cha[i]=d[i]-d[i-1];
		if(cha[i]>as){
			cout<<"No Solution";
			return 0;
		}
	}
	jg[0]=p;
	d[n+1]=m;
	jg[n+1]=0;
	d[0]=0;
	if(m-d[n]>as){
    	cout<<"No Solution";
		return 0;
	}
	for(int i=0;i<=n;i++){
		if(i>=1){
			sy-=(d[i]-d[i-1])/w;
		}
		f=0;
		for(int j=i+1;j<=n+1;j++){
			if(d[j]-d[i]>as){
				break;
			}
			if(jg[j]<jg[i]){
				f=j;
				break;
			}
		}
		if(f==0){
			ans+=(v-sy)*jg[i];
			sy=v;
		}else{
			if((d[f]-d[i])/w-sy<0){
				continue;
			}
			ans+=jg[i]*((d[f]-d[i])/w-sy);
			sy+=(d[f]-d[i])/w;
		}
	}
	printf("%.2lf",ans);
	return 0;
}

接下来对这个代码进行一些分析。
首先发现源文件内声明了一个从来没有用过的maxx函数。
其次不建议使用同时使用和时关闭流同步。(事实上由于输出操作只有一次,所以可以通过本题。)

同学提供了如下hack数据:
输入:

87.75 13.03 5.75 7.29 3 
22.10 7.38 
24.21 6.81 
82.08 6.96 

标准输出:

105.95

实际输出103.45。经过调试发现在i=1时进入了语句:

if(i>=1){
	sy-=(d[i]-d[i-1])/w;
	}

但是sy并没有如预期中减小。经检查变量的值均无误。
推测是精度缺失,将式中变量声明为long double,结果不变。
仍然推测是精度缺失,将上文中的赋值语句修改如下:

sy-=(d[i]-d[i-1])*100.0/w/100.0

成功通过该题。
将式中变量的类型改回double,无法通过该题。

小细节

1.对于long double类型,注意printf要对应修改为%llf。(听说有的地方不支持%llf,也可以在最后一步将ans强制类型转换为double)
2.题目中要求四舍五入,恰好printf就是四舍五入,因此不需要额外设置。

修改之后的代码

#include<bits/stdc++.h>
using namespace std;
long double d[10],jg[10],cha[10],as,ans;
int f; 
int main(){
	long double sy=0,m,v,w,p;
	int n;
	cin>>m>>v>>w>>p>>n;
	as=v*w;
	for(int i=1;i<=n;i++){
		cin>>d[i]>>jg[i];
		cha[i]=d[i]-d[i-1];
		if(cha[i]>as){
			cout<<"No Solution";
			return 0;
		}
	}
	jg[0]=p;
	d[n+1]=m;
	jg[n+1]=0;
	d[0]=0;
	if(m-d[n]>as){
    	cout<<"No Solution";
		return 0;
	}
	for(int i=0;i<=n;i++){
		if(i>=1){
			sy-=(d[i]-d[i-1])*100.0/w/100.0;
//			cout<<"***i="<<i<<":"<<(d[i]-d[i-1])/w<<"***";
		}
//		cout<<"i="<<i<<" now oil="<<sy<<" ";
		f=0;
		for(int j=i+1;j<=n+1;j++){
			if(d[j]-d[i]>as){
				break;
			}
			if(jg[j]<jg[i]){
				f=j;
				break;
			}
		}
		if(f==0){
			ans+=(v-sy)*jg[i];
			sy=v;
		}else{
			if((d[f]-d[i])/w-sy<=0){
				continue;
			}
//			cout<<"*f="<<f<<"*";
			ans+=jg[i]*((d[f]-d[i])/w-sy);
//			cout<<"sy:"<<sy<<"to";
			sy+=(d[f]-d[i])/w;
//			cout<<sy<<" ";
		}
//		cout<<"ans="<<ans<<endl;
	}
//	for(int i=0;i<=n+1;i++) {
//		cout<<i<<":"<<d[i]<<endl;
//	}
	printf("%.2llf",ans);
	return 0;
}

2.AtCoder Beginner Contest 270 ,F - Transportation

主要思路:
将飞机场、港口和桥(我将road称为桥)进行处理,飞机场和港口分别建立超级源点,桥直接连边,跑最小生成树即为答案。
但是这道题会有特殊情况,比如其实不需要建飞机场,那么就不应该向代表飞机场的超级源点连边,但是因为最小生成树算法的关系,还是不得不向其连一条边,因此需要进行特判。

实际的解决方案:
特判太麻烦了,考虑到时间限制为4s,且情况只有以下七种,考虑手动枚举这七种情况,判断图连通性后答案取最小值:

  1. 飞机
  2. 港口
  3. 桥、飞机
  4. 桥、港口
  5. 飞机、港口
  6. 桥、飞机、港口

如下代码可以通过90%数据(通过36组):

#include<iostream>
#include<queue>
#include<cstring>
#include<vector>
using namespace std;
typedef pair<long long,long long> dir;
typedef pair<long long,pair<long long,long long>> edge;
vector<vector<dir>> a;
bool vis[200005];
long long n,m;
long long Prim() {
	memset(vis,0,sizeof vis);
	priority_queue<edge,vector<edge>,greater<edge>> q;
	long long sum=0;
	q.push(edge{0,pair<long long,long long>{0,1}});
	while(!q.empty()) {
		long long fa=q.top().second.first;
		long long u=q.top().second.second;
		long long w=q.top().first;
		q.pop();
		if(vis[u]) continue;
		sum+=w;
		vis[u]=true;
		for(auto&i:a[u])
			if(!vis[i.second]) 
				q.push(edge{i.first,pair<long long,long long>{u,i.second}});
	}
	bool flag=true;
	for(long long i=1;i<=n;i++)
		flag=flag&&vis[i];
	if(!flag) return 1e18;
	return sum;
}

long long solve(vector<vector<edge>> b) {
	a.clear();
	for(auto&j:b)
		for(auto&i:j)
			a[i.second.first].push_back(dir{i.first,i.second.second});
	return Prim();
}
int main() {
	cin>>n>>m;
	for(long long i=0;i<=n+2;i++)
		a.push_back(vector<dir>{});
	vector<edge> sky;
	for(long long i=1;i<=n;i++) {
		long long w;
		cin>>w;
		sky.push_back(edge{w,pair<long long,long long>{i,n+2}});
		sky.push_back(edge{w,pair<long long,long long>{n+2,i}});
	}
	vector<edge> harbor;
	for(long long i=1;i<=n;i++) {
		long long w;
		cin>>w;
		harbor.push_back(edge{w,pair<long long,long long>{i,n+1}});
		harbor.push_back(edge{w,pair<long long,long long>{n+1,i}});
	}
	vector<edge> road;
	for(long long i=1;i<=m;i++) {
		long long u,v,w;
		cin>>u>>v>>w;
		road.push_back(edge{w,pair<long long,long long>{u,v}});
		road.push_back(edge{w,pair<long long,long long>{v,u}});
	}
	long long ans=solve(vector<vector<edge>>{road});
	ans=min(ans,solve(vector<vector<edge>>{sky}));
	ans=min(ans,solve(vector<vector<edge>>{harbor}));
	ans=min(ans,solve(vector<vector<edge>>{road,sky}));
	ans=min(ans,solve(vector<vector<edge>>{road,harbor}));
	ans=min(ans,solve(vector<vector<edge>>{sky,harbor}));
	ans=min(ans,solve(vector<vector<edge>>{sky,harbor,road}));
	cout<<ans<<endl;
	return 0;
}

未通过全部数据的原因是,solve函数中vector容器a在清空之后未重新指定大小。
经过修改,如下代码可以通过全部数据:

#include<iostream>
#include<queue>
#include<cstring>
#include<vector>
using namespace std;
typedef pair<long long,long long> dir;
typedef pair<long long,pair<long long,long long>> edge;
vector<vector<dir>> a;
bool vis[200005];
long long n,m;
long long Prim() {
	memset(vis,0,sizeof vis);
	priority_queue<edge,vector<edge>,greater<edge>> q;
	long long sum=0;
	q.push(edge{0,pair<long long,long long>{0,1}});
	while(!q.empty()) {
		long long fa=q.top().second.first;
		long long u=q.top().second.second;
		long long w=q.top().first;
		q.pop();
		if(vis[u]) continue;
		sum+=w;
		vis[u]=true;
		for(auto&i:a[u])
			if(!vis[i.second]) 
				q.push(edge{i.first,pair<long long,long long>{u,i.second}});
	}
	bool flag=true;
	for(long long i=1;i<=n;i++)
		flag=flag&&vis[i];
	if(!flag) return 1e18;
	return sum;
}

long long solve(vector<vector<edge>> b) {
	a.clear();
	for(long long i=0;i<=n+2;i++)
		a.push_back(vector<dir>{});
	for(auto&j:b)
		for(auto&i:j)
			a[i.second.first].push_back(dir{i.first,i.second.second});
	return Prim();
}
int main() {
	cin>>n>>m;
	vector<edge> sky;
	for(long long i=1;i<=n;i++) {
		long long w;
		cin>>w;
		sky.push_back(edge{w,pair<long long,long long>{i,n+2}});
		sky.push_back(edge{w,pair<long long,long long>{n+2,i}});
	}
	vector<edge> harbor;
	for(long long i=1;i<=n;i++) {
		long long w;
		cin>>w;
		harbor.push_back(edge{w,pair<long long,long long>{i,n+1}});
		harbor.push_back(edge{w,pair<long long,long long>{n+1,i}});
	}
	vector<edge> road;
	for(long long i=1;i<=m;i++) {
		long long u,v,w;
		cin>>u>>v>>w;
		road.push_back(edge{w,pair<long long,long long>{u,v}});
		road.push_back(edge{w,pair<long long,long long>{v,u}});
	}
	long long ans=solve(vector<vector<edge>>{road});
	ans=min(ans,solve(vector<vector<edge>>{sky}));
	ans=min(ans,solve(vector<vector<edge>>{harbor}));
	ans=min(ans,solve(vector<vector<edge>>{road,sky}));
	ans=min(ans,solve(vector<vector<edge>>{road,harbor}));
	ans=min(ans,solve(vector<vector<edge>>{sky,harbor}));
	ans=min(ans,solve(vector<vector<edge>>{sky,harbor,road}));
	cout<<ans<<endl;
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值