牛客練習賽11

B 假的字符串
題意：給定n個字符串，互不相等，你可以任意指定字符之間的大小關系（即重定義字典序），求有多少個串可能成為字典序最小的串，並輸出它們。

tags：好題

對於一個字符串，

1】如有其它字符串是它的前綴，那肯定不可能。這個直接用字典樹處理就可以。

2】但如果以這個字符串為最小，怎麽判定其它字符串不會矛盾呢？

其實矛盾的情況詳細一點說是： 比如要以 abcd 為最小， 但又有另一個字符串 aba ，這就矛盾了。

對這種情況，在跑字典樹的時候，我們對有相同父親結點的多個兒子結點相互連邊，然後每次拓撲排序判一下環即可。

#include<bits/stdc++.h>
using
 
 namespace std;
#pragma comment(linker, "/STACK:102400000,102400000")
#define rep(i,a,b) for (int i=a; i<=b; ++i)
#define per(i,b,a) for (int i=b; i>=a; --i)
#define mes(a,b)  memset(a,b,sizeof(a))
#define INF 0x3f3f3f3f
#define MP make_pair
#define PB push_back
#define fi  first
#define se  second
typedef 
 
long long ll;
const int N = 30005, M = 300005;

int n, cnt, tr[M][26], val[M*26];
string s[N];
void Insert(int ci)
{
    int len=s[ci].size(), now=0;
    rep(j,0,len-1)
    {
        int tmp=s[ci][j]-‘a‘;
        if(tr[now][tmp]==0) tr[now][tmp]=++cnt;
        now = tr[now][tmp];
    }
    ++val[now];
}
vector
 
< int > G[26];
int in[26];
queue< int > q;
bool Toposort()
{
    rep(i,0,25) if(in[i]==0) q.push(i);
    while(!q.empty())
    {
        int u=q.front(); q.pop();
        for(int i=0; i<G[u].size(); ++i)
        {
            --in[G[u][i]];
            if(in[G[u][i]]==0) q.push(G[u][i]);
        }
    }
    rep(i,0,25) if(in[i]) return false;
    return true;
}
bool check(int ci)
{
    rep(i,0,25) G[i].clear(), in[i]=0;
    int len=s[ci].size(), now=0;
    rep(j,0,len-1)
    {
        int tmp=s[ci][j]-‘a‘;
        if(j<len-1 && val[tr[now][tmp]]) return false;
        rep(i,0,25) {
            if(tr[now][i] && i!=tmp)
                G[tmp].PB(i), ++in[i];
        }
        now = tr[now][tmp];
    }
    if(Toposort()) return true;
    return false;
}
int ans[N], ans1;
int main()
{
    scanf("%d", &n);
    rep(i,1,n)
    {
        cin>> s[i];
        Insert(i);
    }
    rep(i,1,n)
    {
        if(check(i)) ans[++ans1]=i;
    }
    printf("%d\n", ans1);
    rep(i,1,ans1)
        cout<<s[ans[i]]<<endl;

    return 0;
}

E 求最值

題意：

給你一個長為n的序列a

定義f(i,j)=(i-j)²+g(i,j)²

tags：翻譯一下，就是 (i-j)^2+(sum[i]-sum[j])^2 最小， 也就是最近點對。 分治 O(nlogn)

模板：

// 分治求平面最近點對
struct Point { double x, y; };
struct Point point[N], *px[N], *py[N];
double get_dis(Point *p1,Point *p2) {
    return sqrt((p1->x-p2->x)*(p1->x-p2->x)+(p1->y-p2->y)*(p1->y-p2->y));
}
bool cmpx(Point *p1,Point *p2) { return p1->x<p2->x; }
bool cmpy(Point *p1,Point *p2) { return p1->y<p2->y; }
double min(double a,double b) { return a<b?a:b; }
double closest(int s,int e)
{
    if(s+1==e)
        return get_dis(px[s],px[e]);
    if(s+2==e)
        return min(get_dis(px[s],px[s+1]),min(get_dis(px[s+1],px[e]),get_dis(px[s],px[e])));
    int mid=(s+e)>>1;
    double ans=min(closest(s,mid),closest(mid+1,e));//遞歸求解
    int i, j, cnt=0;
    for(i=s; i<=e; i++)//把x坐標在px[mid].x-ans~px[mid].x+ans範圍內的點取出來
    {
        if(px[i]->x>=px[mid]->x-ans && px[i]->x<=px[mid]->x+ans)
            py[cnt++]=px[i];
    }
    sort(py, py+cnt, cmpy);//按y坐標排序
    for(i=0; i<cnt; i++)
    {
        for(j=i+1;j<cnt;j++)//py數組中的點是按照y坐標升序的
        {
            if(py[j]->y-py[i]->y>=ans)  break;
            ans=min(ans,get_dis(py[i],py[j]));
        }
    }
    return ans;
}
/*void Init()
{
    for(int i=1; i<=n; ++i) {
        scanf("%lf%lf", &point[i].x, &point[i].y);
        px[i] = &point[i];
    }
    sort(px+1, px+1+n, cmpx);
}*/

// E
#include<bits/stdc++.h>
using namespace std;
#pragma comment(linker, "/STACK:102400000,102400000")
#define rep(i,a,b) for (int i=a; i<=b; ++i)
#define per(i,b,a) for (int i=b; i>=a; --i)
#define mes(a,b)  memset(a,b,sizeof(a))
#define INF 0x3f3f3f3f
#define MP make_pair
#define PB push_back
#define fi  first
#define se  second
typedef long long ll;
const int N = 100005;

struct Point { double x, y; };
struct Point point[N], *px[N], *py[N];
double get_dis(Point *p1,Point *p2) {
    return sqrt((p1->x-p2->x)*(p1->x-p2->x)+(p1->y-p2->y)*(p1->y-p2->y));
}
bool cmpx(Point *p1,Point *p2) { return p1->x<p2->x; }
bool cmpy(Point *p1,Point *p2) { return p1->y<p2->y; }
double min(double a,double b) { return a<b?a:b; }
double closest(int s,int e)
{
    if(s+1==e)
        return get_dis(px[s],px[e]);
    if(s+2==e)
        return min(get_dis(px[s],px[s+1]),min(get_dis(px[s+1],px[e]),get_dis(px[s],px[e])));
    int mid=(s+e)>>1;
    double ans=min(closest(s,mid),closest(mid+1,e));//遞歸求解
    int i, j, cnt=0;
    for(i=s; i<=e; i++)//把x坐標在px[mid].x-ans~px[mid].x+ans範圍內的點取出來
    {
        if(px[i]->x>=px[mid]->x-ans && px[i]->x<=px[mid]->x+ans)
            py[cnt++]=px[i];
    }
    sort(py, py+cnt, cmpy);//按y坐標排序
    for(i=0; i<cnt; i++)
    {
        for(j=i+1;j<cnt;j++)//py數組中的點是按照y坐標升序的
        {
            if(py[j]->y-py[i]->y>=ans)  break;
            ans=min(ans,get_dis(py[i],py[j]));
        }
    }
    return ans;
}
/*void Init()
{
    for(int i=1; i<=n; ++i) {
        scanf("%lf%lf", &point[i].x, &point[i].y);
        px[i] = &point[i];
    }
    sort(px+1, px+1+n);
}*/

int n;
double ai, sum[N];
int main()
{
    scanf("%d", &n);
    rep(i,1,n)
    {
        scanf("%lf", &ai);
        sum[i]=sum[i-1]+ai;
        point[i].x=1.0*i, point[i].y=sum[i];
        px[i] = &point[i];
    }
    sort(px+1, px+1+n, cmpx);
    double ans = closest(1, n);
    printf("%.0f\n", ans*ans);

    return 0;
}

牛客練習賽11 B trie樹+拓撲判環 E 分治求平面最近點對