wa从一点到晚上11点没停过，也不知道为什么错，第二天换了个思路做，终于过了。这题还是有点问题的，数据有点水，我看到有人贴的代码baabbaab这组数据是4 0，明显错的，但是却可以过。

下面的是我第一次写的代码，也不知道到底哪里错了，能想到的样例都过了。

正的直接后缀数组求出即可。反着的我还用height注意细节。

希望大牛们指出错误，我真的找了几个小时都找不出。

  1. #include <stdio.h>
  2. #include <string.h>
  3. #define maxn 51000
  4. int wa[maxn],wb[maxn],wv[maxn],ws[maxn];
  5. int cmp(int *r,int a,int b,int l)
  6. {return r[a]==r[b]&&r[a+l]==r[b+l];}
  7. int min(int x,int y)
  8. {return x<y?x:y;}
  9. int max(int x,int y)
  10. {return x>y?x:y;}
  11. void da(int *r,int *sa,int n,int m)
  12. {
  13.     int i,j,p,*x=wa,*y=wb,*t;
  14.     for(i=;i<m;i++) ws[i]=;
  15.     for(i=;i<n;i++) ws[x[i]=r[i]]++;
  16.     for(i=;i<m;i++) ws[i]+=ws[i-];
  17.     for(i=n-;i>=;i--) sa[--ws[x[i]]]=i;
  18.     for(j=,p=;p<n;j*=,m=p)
  19.     {
  20.         for(p=,i=n-j;i<n;i++) y[p++]=i;
  21.         for(i=;i<n;i++) if(sa[i]>=j) y[p++]=sa[i]-j;
  22.         for(i=;i<n;i++) wv[i]=x[y[i]];
  23.         for(i=;i<m;i++) ws[i]=;
  24.         for(i=;i<n;i++) ws[wv[i]]++;
  25.         for(i=;i<m;i++) ws[i]+=ws[i-];
  26.         for(i=n-;i>=;i--) sa[--ws[wv[i]]]=y[i];
  27.         for(t=x,x=y,y=t,p=,x[sa[]]=,i=;i<n;i++)
  28.             x[sa[i]]=cmp(y,sa[i-],sa[i],j)?p-:p++;
  29.     }
  30.     return;
  31. }
  32. int Rank[maxn],height[maxn];
  33. void calheight(int *r,int *sa,int n)
  34. {
  35.     int i,j,k=;
  36.     for(i=;i<=n;i++) Rank[sa[i]]=i;
  37.     for(i=;i<n;height[Rank[i++]]=k)
  38.         for(k?k--:,j=sa[Rank[i]-];r[i+k]==r[j+k];k++);
  39.     return;
  40. }
  41. int r[maxn],sa[maxn];
  42. char s[maxn],st[maxn];
  43. int main()
  44. {
  45.     int i,j,t,n;
  46.     scanf("%d",&t);
  47.     while(t--)
  48.     {
  49.         scanf("%d",&n);
  50.         if(n==)
  51.         {
  52.             printf("0 0\n");
  53.             continue;
  54.         }
  55.         scanf("%s",s);
  56.         int len=n;
  57.         for(i=;i<n;i++)
  58.             s[len++]=s[i];
  59.         for(i=;i<len;i++)
  60.             r[i]=s[i];
  61.         r[len]=;
  62.         da(r,sa,len+,);
  63.         calheight(r,sa,len);
  64.         int setone=;
  65.         for(i=;i<=len;i++)
  66.         {
  67.             if(len-sa[i]-<n)
  68.                 continue;
  69.             else if(height[i]>=n)
  70.                 setone=min(sa[i-],sa[i]);
  71.             else
  72.             {
  73.                 setone=sa[i];
  74.             }
  75.         }
  76.         setone%=n;
  77.         //printf("%d\n",setone);
  78.         int fl=;
  79.         for(i=len-;i>=;i--)
  80.             st[fl++]=s[i];
  81.         for(i=;i<len;i++)
  82.             r[i]=st[i];
  83.         r[len]=;
  84.         da(r,sa,len+,);
  85.         calheight(r,sa,len);
  86.         int settwo=;
  87.         for(i=;i<=len;i++)
  88.         {
  89.             if(len-sa[i]-<n)
  90.                 continue;
  91.             else if(height[i]>=n)
  92.                 settwo=max(sa[i-],sa[i]);
  93.             else
  94.             {
  95.                 settwo=sa[i];
  96.             }
  97.         }
  98.         //printf("%d\n",settwo);
  99.         int flag=;
  100.         for(i=;i<n;i++)
  101.         {
  102.             if(s[setone+i]==st[settwo+i])
  103.                 continue;
  104.             else if(s[setone+i]>st[settwo+i])
  105.             {
  106.                 flag=;
  107.                 break;
  108.             }
  109.             else
  110.             {
  111.                 flag=;
  112.                 break;
  113.             }
  114.         }
  115.         settwo=n-settwo%n-;
  116.         if(flag==)
  117.         {
  118.             printf("%d 0\n",setone+);
  119.         }
  120.         else if(flag==)
  121.         {
  122.             printf("%d 1\n",settwo+);
  123.         }
  124.         else
  125.         {
  126.             if(settwo>=setone)
  127.                 printf("%d 0\n",setone+);
  128.             else printf("%d 1\n",settwo+);
  129.         }
  130.     }
  131. }

下面是ac代码，反向通过先在末尾加未出现过的最大字母，再多添加0（为了使0位附0），那么就能避免上面代码出现的情况，即满足条件的在sa[]里面的中间部分，但是不能确定后面是否还有，添加字母后，就能保证能够满足，这样只要在[0，原字符串长度）里找到最大的rank的位置就好。

  1. #include <stdio.h>
  2. #include <string.h>
  3. #define maxn 101000
  4. int wa[maxn],wb[maxn],wv[maxn],ws[maxn];
  5. int cmp(int *r,int a,int b,int l)
  6. {return r[a]==r[b]&&r[a+l]==r[b+l];}
  7. int min(int x,int y)
  8. {return x<y?x:y;}
  9. int max(int x,int y)
  10. {return x>y?x:y;}
  11. void da(int *r,int *sa,int n,int m)
  12. {
  13.     int i,j,p,*x=wa,*y=wb,*t;
  14.     for(i=;i<m;i++) ws[i]=;
  15.     for(i=;i<n;i++) ws[x[i]=r[i]]++;
  16.     for(i=;i<m;i++) ws[i]+=ws[i-];
  17.     for(i=n-;i>=;i--) sa[--ws[x[i]]]=i;
  18.     for(j=,p=;p<n;j*=,m=p)
  19.     {
  20.         for(p=,i=n-j;i<n;i++) y[p++]=i;
  21.         for(i=;i<n;i++) if(sa[i]>=j) y[p++]=sa[i]-j;
  22.         for(i=;i<n;i++) wv[i]=x[y[i]];
  23.         for(i=;i<m;i++) ws[i]=;
  24.         for(i=;i<n;i++) ws[wv[i]]++;
  25.         for(i=;i<m;i++) ws[i]+=ws[i-];
  26.         for(i=n-;i>=;i--) sa[--ws[wv[i]]]=y[i];
  27.         for(t=x,x=y,y=t,p=,x[sa[]]=,i=;i<n;i++)
  28.             x[sa[i]]=cmp(y,sa[i-],sa[i],j)?p-:p++;
  29.     }
  30.     return;
  31. }
  32. int Rank[maxn],height[maxn];
  33. void calheight(int *r,int *sa,int n)
  34. {
  35.     int i,j,k=;
  36.     for(i=;i<=n;i++) Rank[sa[i]]=i;
  37.     for(i=;i<n;height[Rank[i++]]=k)
  38.         for(k?k--:,j=sa[Rank[i]-];r[i+k]==r[j+k];k++);
  39.     return;
  40. }
  41. int r[maxn],sa[maxn];
  42. char s[maxn],st[maxn];
  43. int main()
  44. {
  45.     int i,j,t,n;
  46.     scanf("%d",&t);
  47.     while(t--)
  48.     {
  49.         scanf("%d",&n);
  50.         scanf("%s",s);
  51.         memset(r,,sizeof(r));
  52.         int len=n;
  53.         for(i=;i<n;i++)
  54.             s[len++]=s[i];
  55.         for(i=;i<len;i++)
  56.             r[i]=s[i];
  57.         r[len]=;
  58.         da(r,sa,len+,);
  59.         calheight(r,sa,len);
  60.         int setone=;
  61.         int maxx=-;
  62.         for(i=;i<n;i++)
  63.         {
  64.             if(maxx<Rank[i])
  65.             {
  66.                 maxx=Rank[i];
  67.                 setone=i;
  68.             }
  69.         }
  70.         //printf("%d\n",setone);
  71.         memset(r,,sizeof(r));
  72.         int fl=;
  73.         for(i=len-;i>=;i--)
  74.             st[fl++]=s[i];
  75.         for(i=;i<len;i++)
  76.             r[i]=st[i];
  77.         r[len++]=;//添加一个最大的
  78.         r[len]=;//0位无字符
  79.         da(r,sa,len+,);
  80.         calheight(r,sa,len);
  81.         int settwo=;
  82.         maxx=-;
  83.         for(i=;i<n;i++)
  84.         {
  85.             if(maxx<Rank[i])
  86.             {
  87.                 maxx=Rank[i];
  88.                 settwo=i;
  89.             }
  90.         }
  91.         //printf("%d\n",settwo);
  92.         int flag=;
  93.         for(i=;i<n;i++)
  94.         {
  95.             if(s[setone+i]==st[settwo+i])
  96.                 continue;
  97.             else if(s[setone+i]>st[settwo+i])
  98.             {
  99.                 flag=;
  100.                 break;
  101.             }
  102.             else
  103.             {
  104.                 flag=;
  105.                 break;
  106.             }
  107.         }
  108.         settwo=n-settwo%n-;
  109.         if(flag==)
  110.         {
  111.             printf("%d 0\n",setone+);
  112.         }
  113.         else if(flag==)
  114.         {
  115.             printf("%d 1\n",settwo+);
  116.         }
  117.         else
  118.         {
  119.             if(settwo>=setone)
  120.                 printf("%d 0\n",setone+);
  121.             else printf("%d 1\n",settwo+);
  122.         }
  123.     }
  124.     return ;
  125. }

hdu5442 Favorite Donut 后缀数组 长春网赛的更多相关文章

  1. hdu 5442 Favorite Donut 后缀数组

    Favorite Donut Time Limit: 1 Sec Memory Limit: 256 MB 题目连接 http://acm.hdu.edu.cn/showproblem.php?pid ...

  2. 2013长春网赛1005 hdu 4763 Theme Section（kmp应用）

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4763 题意:给出一个字符串,问能不能在该串的前中后部找到相同的子串,输出最长的字串的长度. 分析:km ...

  3. hihocoder1236(2015长春网赛J题) Scores(bitset && 分块)

    题意:给你50000个五维点(a1,a2,a3,a4,a5),50000个询问(q1,q2,q3,q4,q5),问已知点里有多少个点(x1,x2,x3,x4,x5)满足(xi<=qi,i=1,2 ...

  4. ACM学习历程—HDU 5446 Unknown Treasure(数论)（2015长春网赛1010题）

    Problem Description On the way to the next secret treasure hiding place, the mathematician discovere ...

  5. hdu 4764 && 2013长春网赛题解

    一个组合游戏题. 解答: 从后面往前面推,首先n-1是必胜位,然后前面的k位是必败位,如此循环下去.所以题目就容易了! 代码: #include<cstdio> using namespa ...

  6. hdu 5441 Travel (2015长春网赛)

    http://acm.hdu.edu.cn/showproblem.php?pid=5441 题目大意是给一个n个城市(点)m条路线(边)的双向的路线图,每条路线有时间值(带权图),然后q个询问,每个 ...

  7. 2013长春网赛1009 hdu 4767 Bell（矩阵快速幂+中国剩余定理）

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4767 题意:求集合{1, 2, 3, ..., n}有多少种划分情况bell[n],最后结果bell[ ...

  8. 2013长春网赛1001 hdu 4759 Poker Shuffle

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4759 题意:有一堆2^n的牌,牌原先按(1,2,....k)排序,每一次洗牌都将牌分成两种情况:(1, ...

  9. 2013长春网赛1004 hdu 4762 Cut the Cake

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4762 题意:有个蛋糕,切成m块,将n个草莓放在上面,问所有的草莓放在同一块蛋糕上面的概率是多少.2 & ...

随机推荐

  1. OCP/OCA Oracle 学习001

    select * from TEST3 t select object_type, count(object_type) from user_objects group by object_type ...

  2. 把char[]数组里面的内容用MessageBox显示出来

    const char *q; q = mysql_get_client_info();  //获得一字符串指针 CString p(q); AfxMessageBox(p);

  3. linux后台执行程序相关命令

    linux下我们如果想一个任务或者程序还后台执行可以使用&,实际上linux还提供了其他任务调度的命令. bg将一个在后台暂停的命令,变成继续执行 fg将后台中的命令调至前台继续运行 jobs ...

  4. Liferay 7：Liferay DXP解决方案

    分享是美德,欢迎探讨技术 这个作者很厉害呀,写的博客都是解决很刁钻问题的.强烈推荐 http://www.liferaysolution.com/2017/06/captcha-recaptcha-w ...

  5. spring源码学习之bean的加载（二）

    这是接着上篇继续写bean的加载过程,好像是有点太多了,因为bean的加载过程是很复杂的,要处理的情况有很多,继续... 7.创建bean 常规的bean的创建时通过doCreateBean方法来实现 ...

  6. Django项目：CRM(客户关系管理系统)--41--33PerfectCRM实现King_admin编辑整张表限制

    readonly_table=False#默认表单不锁定 readonly_table=True#默认表单不锁定 # forms.py # ————————19PerfectCRM实现King_adm ...

  7. php5 常量

    <?php define("GREETING", "Welcome to w3cschool.cn!", true); echo greeting; ?& ...

  8. 正确而又严谨得ajax原生创建方式

    自己去封装一个xhr对象是一件比较麻烦的事情.其实也不麻烦,注意逻辑和一个ie6兼容方案(可无),和一个304 其他2开头的status都可以就好了 <!DOCTYPE html> < ...

  9. 严格模式下顶层箭头函数this指向的是全局对象

    我们知道普通函数调用,this在非严格模式下指向全局对象,在严格模式下是undefined.那箭头函数呢?我们知道,箭头函数没有自己的this,它的this是最近外层非箭头函数的this,那直接在顶层 ...

  10. 【洛谷P1207】双重回文数 【USACO1.2】

    P1207 [USACO1.2]双重回文数 Dual Palindromes 题目描述 如果一个数从左往右读和从右往左读都是一样,那么这个数就叫做"回文数".例如,12321就是一 ...