AtCoder Beginner Contest 254

关注

原创最新推荐文章于 2026-09-14 13:36:30 发布 403 阅读 AI 写同款· GEO 优化 ›

ITK-snap进行分割区域标注、三维显示及改变显示颜色在进行医学图像器官分割时，无论我们是用传统方法还是深度学习，最后肯定要用到手动分割的区域作为金标准来验证我们方法的精度。但是手动标注器官的轮廓是十分耗时的，这里我介绍一个比较简单的标注方法，以肺分割作为例子。我们使用ITK-snap作为分割标注工具。当切片间距离比较小时，比如1mm，一个患者的完整HRCT序列通常在200张以上，如果让医生一一进行分割区域标注是非常耗时乏味的。因此，我们可以使用最简单的阈值法比如OSTU，加上基本的形态学操作对肺进行粗分割得到粗略的分割掩码，在此基础上再手动的对分割结果进阅读详情

D - Together Square

题意：
计算有对少对小于等于 $\leq n \leq 2×10^5)$ 的 $i, j$ ，满足 $i \times j$ 是平方数
思路：
如果 $i$ 是平方数，那么 $j$ 也是平方数才行。
否则， $j$ 应该是平方数乘以 $i$ 中出现次数为奇数的质因子

#include<bits/stdc++.h>
#define pb push_back
typedef long long ll;
using namespace std;
typedef pair<int,int> PII;
const int maxn=2e5+7;
int n,vis[maxn];
int main()
{
	scanf("%d",&n);
	vector<int>v;
	for(int i=1;;i++)
	{
		if(i*i>n) break;
		vis[i*i]=1;v.pb(i*i);
	}
	ll ans=0;
	for(int i=1;i<=n;i++)
	{
		if(vis[i]) ans+=v.size();
		else
		{
			int nn=i;
			for(int j=1;j<v.size();j++)
			{
				while(nn%v[j]==0)
				{
					nn/=v[j];
				}
			}
			for(int j=0;j<v.size();j++)
			{
				if(v[j]*nn>n) break;
				ans++;
			}
		}
	}
	printf("%lld\n",ans);
}

F - Rectangle GCD

题意：
有长度为 $n(1≤n≤2×10^5)$ 的数组 $A,B(1≤A _i,B_i≤10^9)$ ，以及一个 $n \times n$ 的正方形，第 $i$ 行第 $j$ 列的数为 $A_i+B_j$
有 $q(1≤q≤2×10^5)$ 次询问，每次询问一个矩形区域内所有数字的 $\gcd$
思路：
矩形内数字为
$a_i+b_j \quad a_{i}+b_{j+1} \quad a_{i}+b_{j+2} \quad a_{i}+b_{j+3} \quad ...\\ a_{i+1}+b_{j} \quad a_{i+1}+b_{j+1} \quad a_{i+1}+b_{j+3} \quad a_{i+1}+b_{j+2} \quad ...\\ ......$
已知性质
$\gcd(a,b,c,d)=\gcd(a,b-a,c-b,d-c)$
所以问题转化为 $\gcd$ 数组 $B, A$ 的差分数组
用线段树维护即可

#include<bits/stdc++.h>
#define pb push_back
typedef long long ll;
using namespace std;
typedef pair<int,int> PII;
const int maxn=2e5+7;
struct gcd_tree
{
	int nums[maxn];
	struct tree
	{
		int l,r,sum;
		int mid()
		{
		return (l+r)/2;
		}
	}tre[maxn<<2];
	void pushup(int num)
	{
		tre[num].sum=__gcd(tre[2*num].sum,tre[2*num+1].sum);
	}
	void build(int num,int l,int r)
	{
		tre[num].l=l,tre[num].r=r;
		if(l==r)
		{
			tre[num].sum=nums[l];return;
		}
		int mid=tre[num].mid();
		build(2*num,l,mid);
		build(2*num+1,mid+1,r);
		pushup(num);
	}
	int query(int num,int l,int r)
	{
		if(tre[num].l==l && tre[num].r==r)
		return tre[num].sum;
		int mid=tre[num].mid();
		if(l>mid)
		return query(2*num+1,l,r);
		else if(r<=mid)
		return query(2*num,l,r);
		else
		return __gcd(query(2*num,l,mid),query(2*num+1,mid+1,r));
	}
}da,db;
int n,m;
int a[maxn],b[maxn];
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	for(int i=1;i<=n;i++) scanf("%d",&b[i]);
	for(int i=1;i<=n;i++) da.nums[i]=abs(a[i]-a[i-1]);
	da.build(1,1,n);
	for(int i=1;i<=n;i++) db.nums[i]=abs(b[i]-b[i-1]);
	db.build(1,1,n);
	while(m--)
	{
		int h1,h2,w1,w2;scanf("%d%d%d%d",&h1,&h2,&w1,&w2);
		int ans=a[h1]+b[w1];
		if(h2>h1) ans=__gcd(ans,da.query(1,h1+1,h2));
		if(w2>w1) ans=__gcd(ans,db.query(1,w1+1,w2));
		printf("%d\n",ans);
	}
}