Comet OJ Contest 4

博客包含几类问题分析。签到题未详述;函数方面,k为奇数时函数值为1,偶数时每k + 1个出现一个0;矩阵问题通过暴力枚举横竖切刀数,用二维前缀和验证;打表题利用10k~10k + 9范围规律计算零散部分。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

  A:签到。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 1000000010
char getc(){char c=getchar();while ((c<'A'||c>'Z')&&(c<'a'||c>'z')&&(c<'0'||c>'9')) c=getchar();return c;}
int gcd(int n,int m){return m==0?n:gcd(m,n%m);}
int read()
{
	int x=0,f=1;char c=getchar();
	while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
	while (c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
	return x*f;
}
int T,a[6],cnt[6];
signed main()
{
	T=read();
	while (T--)
	{
		for (int i=1;i<=5;i++) a[i]=read();
		memset(cnt,0,sizeof(cnt));
		for (int i=1;i<=5;i++) cnt[a[i]]++;
		int mx=0;
		for (int i=1;i<=5;i++) if (cnt[i]>cnt[mx]) mx=i;
		cout<<mx<<endl;
	}
	return 0;
	//NOTICE LONG LONG!!!!!
}

  B:k是奇数时函数值均为1,k是偶数时每k+1个出现一个0。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 1000000010
char getc(){char c=getchar();while ((c<'A'||c>'Z')&&(c<'a'||c>'z')&&(c<'0'||c>'9')) c=getchar();return c;}
int gcd(int n,int m){return m==0?n:gcd(m,n%m);}
ll read()
{
	ll x=0,f=1;char c=getchar();
	while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
	while (c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
	return x*f;
}
int T;
ll l,r,k;
ll f(ll x)
{
	return x-(x+1)/(k+1);
}
signed main()
{
	T=read();
	while (T--)
	{
		l=read(),r=read(),k=read();
		if (k&1) cout<<r-l+1<<endl;
		else cout<<f(r)-f(l-1)<<endl;
	}
	return 0;
	//NOTICE LONG LONG!!!!!
}

  C:暴力枚举横竖各切多少刀,将矩阵压成一行可以得到该情况下列的划分位置,压成一列可以得到该情况下行的划分位置,然后二维前缀和暴力验证即可。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 1000000010
#define N 1010
char getc(){char c=getchar();while ((c<'A'||c>'Z')&&(c<'a'||c>'z')&&(c<'0'||c>'9')) c=getchar();return c;}
int gcd(int n,int m){return m==0?n:gcd(m,n%m);}
ll read()
{
	ll x=0,f=1;char c=getchar();
	while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
	while (c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
	return x*f;
}
int T,n,m,k,tot,a[N][N],s[N],s2[N],s3[N][N],ans[N<<1],way[N<<1];
int calc(int x,int y,int l,int r)
{
	return s3[y][r]-s3[x-1][r]-s3[y][l-1]+s3[x-1][l-1];
}
signed main()
{
	T=read();
	while (T--)
	{
		n=read(),m=read(),k=read();tot=0;
		for (int i=1;i<=n;i++)
			for (int j=1;j<=m;j++)
			{
				char c=getc();
				if (c=='0') a[i][j]=0;else a[i][j]=1,tot++;
			}
		for (int i=1;i<=n;i++)
		{
			s[i]=s[i-1];
			for (int j=1;j<=m;j++)
			if (a[i][j]) s[i]++;
		}
		for (int i=1;i<=m;i++)
		{
			s2[i]=s2[i-1];
			for (int j=1;j<=n;j++)
			if (a[j][i]) s2[i]++;
		}
		for (int i=1;i<=n;i++) 
			for (int j=1;j<=m;j++)
			s3[i][j]=s3[i-1][j]+s3[i][j-1]-s3[i-1][j-1]+1-a[i][j];
		tot=n*m-tot;
		bool isac=0;
		for (int i=1;i<=k;i++) ans[i]=2010;
		for (int i=0;i<=k;i++)
		if (tot%((i+1)*(k-i+1))==0)
		{
			int sum=tot/((i+1)*(k-i+1));
			int last=0;
			int cnt=0;
			bool flag=1;
			for (int j=1;j<n;j++)
			if (cnt<i)
			{
				if ((j-last)*m-(s[j]-s[last])==sum*(k-i+1))
				{
					way[++cnt]=j;
					last=j;
				}
				else if ((j-last)*m-(s[j]-s[last])>sum*(k-i+1)) {flag=0;break;}
			}
			if (cnt<i) flag=0;
			if (flag)
			{
				last=0;
				for (int j=1;j<m;j++)
				if (cnt<k)
				{
					if ((j-last)*n-(s2[j]-s2[last])==sum*(i+1))
					{
						way[++cnt]=n+j-1;
						last=j;
					}
					else if ((j-last)*n-(s2[j]-s2[last])>sum*(i+1)) {flag=0;break;}
				}
			}
			if (cnt<k) flag=0;
			if (flag)
			{
				bool u=1;
				for (int x=1;x<=i;x++)
				{
					for (int y=i+1;y<=k;y++)
					{
						if (calc(way[x-1]+1,way[x],(y==i+1?1:way[y-1]-n+2),way[y]-n+1)!=sum) {u=0;break;}
					}
					if (!u) break;
				}
				if (u)
				{	
					isac=1;
					for (int j=1;j<=k;j++)
					{
						if (way[j]<ans[j])
						{
							for (int x=1;x<=k;x++) ans[x]=way[x];
							break;
						}
						else if (way[j]>ans[j]) break;
					}
				}
			}
		}
		if (!isac) printf("Impossible\n");
		else
		{
			for (int i=1;i<k;i++) printf("%d ",ans[i]);
			printf("%d\n",ans[k]);
		}
	}
	return 0;
	//NOTICE LONG LONG!!!!!
}

  D:打表可知10k~10k+9范围内0~9各出现一次,于是只需要计算零散部分。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 1000000010
char getc(){char c=getchar();while ((c<'A'||c>'Z')&&(c<'a'||c>'z')&&(c<'0'||c>'9')) c=getchar();return c;}
int gcd(int n,int m){return m==0?n:gcd(m,n%m);}
ll read()
{
	ll x=0,f=1;char c=getchar();
	while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
	while (c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
	return x*f;
}
ll l,r,a[30],f[1000000];
int ff(ll n)
{
	if (n<1000000) return f[n];
	int cnt=0;
	while (n) a[++cnt]=n%10,n/=10;
	for (int j=1;j<cnt;j++) a[j]+=a[j+1];
	a[cnt]=0;
	ll t=0;
	for (int j=cnt;j>=1;j--) t=t*10+a[j]%10;
	return ff(t);
}
ll calc(ll n)
{
	int s=0;int w=ff(n/10*10);
	for (ll i=n/10*10;i<=n;i++) s+=(w+(i-n/10*10))%10;
	return n/10*45+s;
}
signed main()
{
	for (int i=0;i<=9;i++) f[i]=i;
	for (int i=10;i<=1000000;i++)
	{
		ll cnt=0,x=i;
		while (x) a[++cnt]=x%10,x/=10;
		for (int j=1;j<cnt;j++) a[j]+=a[j+1];
		a[cnt]=0;
		ll t=0;
		for (int j=cnt;j>=1;j--) t=t*10+a[j]%10;
		f[i]=f[t];
	}
	int T=read();
	while (T--)
	{
		l=read(),r=read();
		cout<<calc(r)-calc(l-1)<<endl;
	}
	return 0;
	//NOTICE LONG LONG!!!!!
}

  先咕着。

 

转载于:https://www.cnblogs.com/Gloid/p/10936526.html

内容概要:本文档详细介绍了基于MATLAB实现多目标差分进化(MODE)算法进行无人机三维路径规划的项目实例。项目旨在提升无人机在复杂三维环境中路径规划的精度、实时性、多目标协调处理能力、障碍物避让能力和路径平滑性。通过引入多目标差分进化算法,项目解决了传统路径规划算法在动态环境和多目标优化中的不足,实现了路径长度、飞行安全距离、能耗等多个目标的协调优化。文档涵盖了环境建模、路径编码、多目标优化策略、障碍物检测与避让、路径平滑处理等关键技术模块,并提供了部分MATLAB代码示例。 适合人群:具备一定编程基础,对无人机路径规划和多目标优化算法感兴趣的科研人员、工程师和研究生。 使用场景及目标:①适用于无人机在军事侦察、环境监测、灾害救援、物流运输、城市管理等领域的三维路径规划;②通过多目标差分进化算法,优化路径长度、飞行安全距离、能耗等多目标,提升无人机任务执行效率和安全性;③解决动态环境变化、实时路径调整和复杂障碍物避让等问题。 其他说明:项目采用模块化设计,便于集成不同的优化目标和动态环境因素,支持后续算法升级与功能扩展。通过系统实现和仿真实验验证,项目不仅提升了理论研究的实用价值,还为无人机智能自主飞行提供了技术基础。文档提供了详细的代码示例,有助于读者深入理解和实践该项目。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值