sicily 1023. Funny Contest

很久之前做了，但是弄了好久才弄懂，现在看一下后缀数组又懵了。。。我弱爆了。。

题目就是按照后缀数组的倍增算法，因为倍增算法用倍增算法排序后，刚好可以统计这次得到的名次。

不过原始的后缀数组不是一个循环的字符串，我们把它改成循环就行了

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;


const int maxn = 110000;
struct node{
    int w;
    int loca;
}data[maxn];
bool ccmp(node a,node b){
    return a.w<b.w;
}
int r[maxn];
int sa[maxn];
int c[maxn];
int k[maxn];
int tt[maxn];
int n;
int cmp(int *r,int a,int b,int l)
{return r[a]==r[b]&&r[(a+l)%n]==r[(b+l)%n];}
int wa[maxn],wb[maxn],wv[maxn],wws[maxn];
void da(int *r,int *sa,int n,int m)
{
    memset(c,0,sizeof(c));
    ///测试数据aabaaaab
    int i,j,p,*x=wa,*y=wb,*t;
    for(i=0;i<m;i++) wws[i]=0;
    for(i=0;i<n;i++) wws[x[i]=r[i]]++;
    for(i=1;i<m;i++) wws[i]+=wws[i-1];
    for(i=n-1;i>=0;i--) sa[--wws[x[i]]]=i;
    memset(tt,0,sizeof(tt));
    tt[1] = 1;
    for(p=1,k[sa[0]]=0,i=1;i<n;i++){
        k[sa[i]]=r[sa[i-1]]==r[sa[i]]?p-1:p++;
        tt[ k[ sa[i] ]+1 ]++;
    }
    //for(int i = 0;i<n;i++)cout <<k[i]<<" ";
    //cout <<endl;
//  for(i = 0;i<n;i++)cout <<k[i]<<" ";
    //cout <<endl;
//  cout <<tt[0]<<" ";
    for(int i = 2;i<n;i++)tt[i] +=tt[i-1];
    //cout <<endl;
    for(i = 0;i<n;i++){
        c[i]+=tt[k[i]];
    /// cout <<c[i]<<" ";
    }


    ///这里sa得到的是只按第一个关键字排序的结果。。
    sa[i]表示排第i的是原来的字符串"aabaaaab"的哪个位置
    如 此时 0 1 3 4 5 6 2 7
    //for(i = 0;i<n;i++)printf("%d ",sa[i]);
    //cout <<endl;
    for(j=1,p=1;j*2<=n;j*=2,m=p)
    {
        按第二个关键字进行排序
        for(p=0,i=n-j;i<n;i++) y[p++]=i;/长度不及2j的第二个，第二个关键字就是没有，那么全部应该排前
        p = 0;
         for(i=0;i<n;i++)  y[p++]=(sa[i]-j+n)%n;
         ///此时sample 结果  y ： 7 0 2 3 4 5 1 6


    //   for(int i = 0;i<n;i++)cout <<y[i]<<" ";
    //   cout <<endl;
         //system("pause");


         //x数组保存的是rank的结果
         for(i=0;i<n;i++) wv[i]=x[y[i]];
         for(i=0;i<m;i++) wws[i]=0;
         for(i=0;i<n;i++) wws[wv[i]]++;
         for(i=1;i<m;i++) wws[i]+=wws[i-1];
         for(i=n-1;i>=0;i--) sa[--wws[wv[i]]]=y[i];
         //for(int i = 0;i<n;i++)cout <<x[i]<<" ";
        // cout <<endl;
         memset(tt,0,sizeof(tt));
         tt[1] = 1;
         for(t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++){
             x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
             tt[x[sa[i]]+1]++;
        //  cout <<p<<endl;;
         }
         for(int i = 2;i<n;i++)tt[i] +=tt[i-1];
         for(i = 0;i<n;i++)c[i]+=tt[x[i]];
        /*for(int i = 0;i<n;i++)cout <<sa[i]<<" ";
        cout <<endl;*/
    //  for(int i = 0;i<n;i++)cout <<x[i]<<" ";
    //  cout <<endl;
        //system("pause");
    }   
        printf("%d",c[0]);
        for( i = 1;i<n;i++)printf(" %d",c[i]);
        printf("\n");
        return;
}


int main(){
    
    
    
    while(scanf("%d",&n)!=EOF){
        
        //int n;
        //scanf("%d",&n);
        for(int i = 0;i<n;i++){
            scanf("%d",&data[i].w);
            data[i].loca = i;
        }
        sort(data,data+n,ccmp);
        int cou =1;
        r[data[0].loca] = 1;
        for(int i = 1;i<n;i++)if(data[i].w==data[i-1].w)r[data[i].loca] = cou;
        else {
            r[data[i].loca] = ++cou;
        }
        da(r,sa,n,n+1);
        ///int l = n;
        //for(int i = 0;i<l;i++)printf("%d ",sa[i]);
    }
return 0;
}                                 

1023. Funny Contest