分治——P3612 [USACO17JAN]Secret Cow Code S+平面最近点对

传送门：[USACO17JAN]Secret Cow Code S - 洛谷

题意:对于一个长度为n的字符串s翻倍长度时将s[n-1]放到s[n]上去，随后将s[0]~s[n-2]的字符全部放到s[n+1]到s[2*n-1]上去，可以无限翻倍，要求第m个位置上的字符是什么

思路: 将翻倍出来的新字符串先看成两个独立的字符串s和s1,比如样例中的 cow 和 wco，

可以看见在s1中s1[i-1]==s[i],当i-1==0时特判s1[0]等于s[n-1],根据这个性质

在s+s1的字符串里面要求第m个位置上的字符(假设m小于等于n*2)，需要先将m-=n，减去一个串的长度如果此时m==0则m变成n,其他的就m--，这两种情况对应上面的==0和!=0的情况

代码：

#include<bits/stdc++.h>
using namespace std;
#define int long long
typedef pair<int,int> PII;
const int N=2e5+10;
signed main(){
    string s;
    int n,m;
    cin>>s>>m;
    int len=s.length();
    for(n=len;n<=m;n<<=1);
    n>>=1;
    m-=1;
    while(n>=len){
        if(m>=n){
            m-=n;
            if(m==0) m=n-1;
            else m-=1;
        }
        n>>=1;
    }
    cout<<s[m%len];
return 0;
}

传送门1: 平面最近点对（加强版） - 洛谷(浮点数)

传送门2:平面最近点对（加强加强版） - 洛谷

思路：分治的思想,将当前一个的平面按照x坐标分成左右两个部分，求左半边的最短距离和右半边的最短距离，再求左边的点与右边的点形成的点对的最短距离(关键步骤)，不断分下去，直到一个部分只有两个点就直接返回两个点的距离，关键在求左右两边的点的最短距离上面

ps:两道题目的代码只有在最后输出有点区别，其余都一样，如果把最短距离变成全局变量可以快一点，如下图所示 

代码1:

#include<bits/stdc++.h>
using namespace std;
#define int long long
struct node{
    int x,y;
}point[400005];
int mpt[400005];
int dis(node a,node b){
return  (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
}
bool cmp1(node a,node b)
{
	if(a.x==b.x)return a.y<b.y;
	return a.x<b.x;
}
bool cmp2(int a,int b)
{
	return point[a].y<point[b].y;
}
int run(int left,int right)
{
	int d=1e18;
	if(left==right)return d;
	if(left+1==right)return dis(point[left],point[right]);
	int mid=(left+right)>>1;
	int d1=run(mid,right); //左边的最短距离
	int d2=run(left,mid-1); //右边的最短距离
	int minn=min(d1,d2);
	int k=0;
	for(int i=left;i<=right;i++)
	{
		if(pow(point[mid].x-point[i].x,2)<=minn)//求有可能凑出比当前最短距离还要短距离的点
		{
			mpt[++k]=i;
		}
	}
	sort(mpt+1,mpt+k+1,cmp2); 
	for(int i=1;i<=k;i++)
	{
		for(int j=i+1;j<=k && pow(point[mpt[j]].y-point[mpt[i]].y,2)<=minn;j++)
		{
			if(minn>=dis(point[mpt[i]],point[mpt[j]]))
				minn=dis(point[mpt[i]],point[mpt[j]]);
		}
	}
	return minn;
}
signed main()
{
	int n;
	scanf("%lld",&n);
	for(int i=1;i<=n;i++)
	{
		scanf("%lld%lld",&point[i].x,&point[i].y);
	}
	sort(point+1,point+n+1,cmp1);
	printf("%.4lf",sqrt((double)run(1,n)));
	return 0;
}

 代码2：

#include<bits/stdc++.h>
using namespace std;
#define int long long
struct node{
    int x,y;
}point[400005];
int mpt[400005];
int dis(node a,node b){
return  (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
}
bool cmp1(node a,node b)
{
	if(a.x==b.x)return a.y<b.y;
	return a.x<b.x;
}
bool cmp2(int a,int b)
{
	return point[a].y<point[b].y;
}
int run(int left,int right)
{
	int d=1e18;
	if(left==right)return d;
	if(left+1==right)return dis(point[left],point[right]);
	int mid=(left+right)>>1;
	int d1=run(mid,right);
	int d2=run(left,mid-1);
	int minn=min(d1,d2);
	int k=0;
	for(int i=left;i<=right;i++)
	{
		if(pow(point[mid].x-point[i].x,2)<=minn)
		{
			mpt[++k]=i;
		}
	}
	sort(mpt+1,mpt+k+1,cmp2);
	for(int i=1;i<=k;i++)
	{
		for(int j=i+1;j<=k && pow(point[mpt[j]].y-point[mpt[i]].y,2)<=minn;j++)
		{
			if(minn>=dis(point[mpt[i]],point[mpt[j]]))
				minn=dis(point[mpt[i]],point[mpt[j]]);
		}
	}
	return minn;
}
signed main()
{
	int n;
	scanf("%lld",&n);
	for(int i=1;i<=n;i++)
	{
		scanf("%lld%lld",&point[i].x,&point[i].y);
	}
	sort(point+1,point+n+1,cmp1);
	printf("%lld",run(1,n));
	return 0;
}