题目大意:一个字符串找出有多少 “ 一个半回文串 ”
一个半回文串的定义为,两个回文串恰好a的右边界是b的回文中心,而b的左边界恰好是a的回文中心。
题目思路:有两种做法,主席树和树状数组,主席树慢些但是题目并不想卡主席树。
分析题目可以知道,只要求出有多少对回文串相互包含对方的回文中心就可以了。
做法一:主席树
主席树的做法比较裸,先跑一边马拉车,对于每一个回文中心,查询他到他的右端点这个区间有多少个回文中心的左端点在当前回文中心的左边就可以了。
#include<bits/stdc++.h>
#define ll long long
using namespace std;
const ll MAXN = 5e5+5;
int b[MAXN],a[MAXN];
char s[MAXN],Ma[MAXN*2];
int Mp[MAXN*2];
void Manacher(char s[],int len) /// Manacher
{
int l=0;
Ma[l++]='$';
Ma[l++]='#';
for(int i=0;i<len;i++){
Ma[l++]=s[i];
Ma[l++]='#';
}
Ma[l]=0;
int mx=0,id=0;
for(int i=0;i<l;i++){
Mp[i]=mx>i?min(Mp[2*id-i],mx-i):1;
while(Ma[i+Mp[i]]==Ma[i-Mp[i]])Mp[i]++;
if(i+Mp[i]>mx){
mx=i+Mp[i];
id=i;
}
}
}
ll sum[MAXN*30];
int L[MAXN*30],R[MAXN*30],T[MAXN];
int cnt=0;
int build(int l,int r) ///初始化主席树
{
int rt=++cnt;
sum[rt]=0;
int mid=(l+r)/2;
if(l<r){
L[rt]=build(l,mid);
R[rt]=build(mid+1,r);
}
return rt;
}
int update(int pre,int l,int r,int x) ///建主席树
{
int rt=++cnt;
L[rt]=L[pre];R[rt]=R[pre];sum[rt]=sum[pre]+1;
int mid=(l+r)/2;
if(l<r){
if(x<=mid)L[rt]=update(L[pre],l,mid,x);
else R[rt]=update(R[pre],mid+1,r,x);
}
return rt;
}
ll query(int u,int v,int l,int r,int k) ///查询区间大于k的有多少个
{
ll ret=0;
int mid=(l+r)/2;
if(k<=l)return sum[v]-sum[u];
if(k<=mid)ret += query(L[u],L[v],l,mid,k);
ret+=query(R[u],R[v],mid+1,r,k);
return ret;
}
int main()
{
ll ans;
int t;
scanf("%d",&t);
while(t--)
{
cnt=0;
scanf("%s",s);
int len=strlen(s);
Manacher(s,len);
ans=0;
int n=len;
for(int i=1;i<2*len+2;i++){
if(Ma[i]!='#'){
int pos=i/2,now=(Mp[i]-1)/2;
b[pos]=min(pos+now,n);a[pos]=max(1,pos-now);
}
}
T[0]=0;
for(int i=1;i<=n;i++){
T[i]=update(T[i-1],1,n,b[i]);
}
for(int i=1;i<=n;i++){
ans+=query(T[a[i]-1],T[i-1],1,n,i);
}
printf("%lld\n",ans);
}
}
做法二:树状数组
做法还是很巧妙的,问题就出在怎么解决两个相互包含的情况,怎么解决相互这个问题,其实和排序降维的思想很相似,我们可以创建一个vector[i]存下当前点为左端点的回文中心编号,然后顺序扫一遍vector,加入树状数组,累计答案只需要查询当前回文中心到左端点这个区间的和就好了。
#include<bits/stdc++.h>
#define ll long long
using namespace std;
const ll MAXN = 5e5+5;
int b[MAXN],a[MAXN];
char s[MAXN],Ma[MAXN*2];
int Mp[MAXN*2],c[MAXN];
void add(int x,int k)
{
for(;x<=MAXN;x+=x&-x)
c[x]+=k;
}
ll ask(int x)
{
ll ans=0;
for(;x;x-=x&-x)
ans+=c[x];
return ans;
}
void Manacher(char s[],int len) /// Manacher
{
int l=0;
Ma[l++]='$';
Ma[l++]='#';
for(int i=0;i<len;i++){
Ma[l++]=s[i];
Ma[l++]='#';
}
Ma[l]=0;
int mx=0,id=0;
for(int i=0;i<l;i++){
Mp[i]=mx>i?min(Mp[2*id-i],mx-i):1;
while(Ma[i+Mp[i]]==Ma[i-Mp[i]])Mp[i]++;
if(i+Mp[i]>mx){
mx=i+Mp[i];
id=i;
}
}
}
ll sum[MAXN*30];
int L[MAXN*30],R[MAXN*30],T[MAXN];
int cnt=0;
vector<int>v[MAXN];
int main()
{
ll ans;
int t;
scanf("%d",&t);
while(t--)
{
cnt=0;
scanf("%s",s);
memset(c,0,sizeof(c));
int len=strlen(s);
Manacher(s,len);
ans=0;
int n=len;
for(int i=1;i<=n;i++){
v[i].clear();
}
for(int i=1;i<2*len+2;i++){
if(Ma[i]!='#'){
int pos=i/2,now=(Mp[i]-1)/2;
b[pos]=min(pos+now,n);
a[pos]=max(1,pos-now);
v[a[pos]].push_back(pos);
}
}
for(int i=1;i<=n;i++){
int len=v[i].size();
for(int j=0;j<len;j++){
int now=v[i][j];
add(now,1);
}
ans+=ask(b[i])-ask(i);
}
printf("%lld\n",ans);
}
}
B - K-th Number
题目大意:n个数字,所有长度大于k的区间都把第k大扔进另一个数组,求另一个数组的第m大。
题目思路:
首先想到的肯定是计算每个数作为第二个数组中的数的贡献次数,找第m就行了。但是每个作为第k大的贡献并不好算,想到一个n方的算法就是枚举每个数字,只要比他大的就可能造成贡献,所以把比它小的改成0,比它大的改成1,尺取一下可以得到有多少比它大的数字作为b数组的贡献值。结果发现是具有二分性的,如果选的这个数字越小,那么他对b数组第m大的贡献数就越多。
#include<bits/stdc++.h>
#define ll long long
using namespace std;
ll n,m,k;
const ll MAXN = 1e5+5;
ll a[MAXN],b[MAXN];
bool check(ll x)
{
ll sum = 0;
ll l = 1,r = 1,num = 0;
while(l<=r&&r<=n){
num += (a[r] >= x);
while(num>=k && l<=r){
num -= (a[l] >= x);
l++;
}
sum += l-1;
r++;
}
if(sum >= m)return 1;
else return 0;
}
int main()
{
ll T;
scanf("%lld",&T);
while(T--){
scanf("%lld%lld%lld",&n,&k,&m);
for(ll i=1;i<=n;i++){
scanf("%lld",&a[i]);
b[i] = a[i];
}
sort(b+1,b+1+n);
ll l = 1,r = n,ans;
while(l<=r){
ll mid = (l+r)/2;
if(check(b[mid])){
l = mid+1;
ans = mid;
}
else r = mid-1;
}
printf("%lld\n",b[ans]);
}
}
//
M - Geometry Problem
题目大意:n个点,求一个圆心和半径满足,至少有一半的点在圆上。
题目思路:
由于题目保证存在这样的点了,那么三个点可以确定一个圆,只要每次从所给点中取出三个点判断是否满足题意就可以了,每次取成功的概率大概是1 / 8,期望是8次可以成功,但是也有脸黑的情况,不过脸再黑,100次总要出来的。
#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int MAXN = 1e5+5;
const double eps = 1e-6;
int n;
struct Point
{
double x,y;
Point(){}
Point(double a,double b){
x = a;y = b;
}
}a[MAXN];
double dis(Point a,Point b)
{
return sqrt( (a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y) );
}
Point Getcir(Point a,Point b,Point c){//给予三个点,求圆心。
double t1=a.x*a.x+a.y*a.y;
double t2=b.x*b.x+b.y*b.y;
double t3=c.x*c.x+c.y*c.y;
double temp=a.x*b.y+b.x*c.y+c.x*a.y-a.x*c.y-b.x*a.y-c.x*b.y;
double x=(t2*c.y+t1*b.y+t3*a.y-t2*a.y-t3*b.y-t1*c.y)/temp/2;
double y=(t3*b.x+t2*a.x+t1*c.x-t1*b.x-t2*c.x-t3*a.x)/temp/2;
return Point(x,y);
}
bool check(Point x,double dist)
{
int cnt = 0;
for(int i=1;i<=n;i++){
if(fabs(dis(a[i],x) - dist) <eps)cnt++;
}
if(cnt>=ceil(n/2.0))return 1;
else return 0;
}
int main()
{
int T;
scanf("%d",&T);
while(T--)
{
scanf("%d",&n);
for(int i=1;i<=n;i++){
scanf("%lf%lf",&a[i].x,&a[i].y);
}
if(n==1){
printf("%.10lf %.10lf %.10lf\n",a[1].x+1,a[1].y,1.0);
}
else if(n<=4){
printf("%.10lf %.10lf %.10lf\n",(a[1].x+a[2].x)/2.0,(a[1].y+a[2].y)/2.0,dis(a[1],a[2])/2.0);
}
else{
Point ans ;
double r;
int _=0;
for(int i=1;i<100;i++){
srand(i);
int x = rand()%(n+1),y = rand()%(n+1),z = rand()%(n+1);
if( (a[x].x-a[y].x)*(a[y].y-a[z].y) - (a[x].y-a[y].y)*(a[y].x-a[z].x) <eps ){
ans = Getcir(a[x],a[y],a[z]);
double dist = dis(a[x],ans);
if(check(ans,dist)){
r = dist;
break;
}
}
}
printf("%.10lf %.10lf %.10lf\n",ans.x,ans.y,r);
}
}
}