记录编号 161053 评测结果 AAAAAAAAAAAAAAAAAAAA
题目名称 [ZJOI 2015]诸神眷顾的幻想乡 最终得分 100
用户昵称 GravatarAsm.Def 是否通过 通过
代码语言 C++ 运行时间 1.937 s
提交时间 2015-04-30 21:47:41 内存使用 113.01 MiB
显示代码纯文本
  1. /***********************************************************************/
  2. /**********************By Asm.Def-Wu Jiaxin*****************************/
  3. /***********************************************************************/
  4. #include <cstdio>
  5. #include <cstring>
  6. #include <cstdlib>
  7. #include <cctype>
  8. #include <algorithm>
  9. using namespace std;
  10. #define SetFile(x) ( freopen(#x".in", "r", stdin), freopen(#x".out", "w", stdout) );
  11. #define UseFREAD
  12. #ifdef UseFREAD
  13. #define getc() *(file_ptr++)
  14. #define FreadLenth 5000000
  15. char CHARPOOL[FreadLenth], *file_ptr = CHARPOOL;
  16. #else
  17. #define getc() getchar()
  18. #endif
  19. #ifdef DEBUG
  20. #include <ctime>
  21. #endif
  22. template<class T>inline void getd(T &x){
  23. char ch = getc();bool neg = false;
  24. while(!isdigit(ch) && ch != '-')ch = getc();
  25. if(ch == '-')ch = getc(), neg = true;
  26. x = ch - '0';
  27. while(isdigit(ch = getc()))x = x * 10 - '0' + ch;
  28. if(neg)x = -x;
  29. }
  30. /***********************************************************************/
  31. const int maxn = 100003, maxs = 1000000;
  32. typedef long long LL;
  33. int N, c, *adj[maxn], d[maxn], col[maxn], leaf[22], lcnt;
  34.  
  35. struct SAM{
  36. SAM *son[10], *link;
  37. int len;
  38. LL cnt;
  39. }Node[maxs<<1], *iter = Node + 1;
  40.  
  41. inline SAM * insert(SAM *last, int ch){
  42. SAM *cur = iter++, *it, *s, *clone;
  43. cur->len = last->len + 1;
  44. it = last;
  45. while(it && !it->son[ch])it->son[ch] = cur, it = it->link;
  46. if(!it)cur->link = Node;
  47. else{
  48. if(it->len + 1 == (s = it->son[ch])->len)cur->link = s;
  49. else{
  50. clone = iter++;*clone = *s;clone->len = it->len + 1;
  51. cur->link = s->link = clone;
  52. do{
  53. it->son[ch] = clone;
  54. it = it->link;
  55. }while(it && it->son[ch] == s);
  56. }
  57. }
  58. return cur;
  59. }
  60.  
  61. void Build(SAM *last, int cur, int p){
  62. SAM *tmp = insert(last, col[cur]);
  63. int son;
  64. for(int i = 0;i < d[cur];++i)if((son = adj[cur][i]) != p)
  65. Build(tmp, adj[cur][i], cur);
  66. }
  67.  
  68. inline void init(){
  69. getd(N), getd(c);
  70. size_t intsize = sizeof(int);
  71. int i, a, b, u[maxn], v[maxn], it[maxn] = {0};
  72. for(i = 1;i <= N;++i)getd(col[i]);
  73. for(i = 1;i < N;++i){
  74. getd(u[i]), getd(v[i]);
  75. ++d[u[i]], ++d[v[i]];
  76. }
  77. for(i = 1;i <= N;++i){
  78. adj[i] = (int*)malloc(intsize * (d[i] + 1));
  79. if(d[i] == 1)leaf[lcnt++] = i;
  80. }
  81. for(i = 1;i < N;++i){
  82. a = u[i], b = v[i];
  83. adj[a][it[a]++] = b;
  84. adj[b][it[b]++] = a;
  85. }
  86. for(i = 0;i < lcnt;++i)Build(Node, leaf[i], 0);
  87. }
  88.  
  89. int main(){
  90.  
  91. #ifdef DEBUG
  92. freopen("test.txt", "r", stdin);
  93. #elif !defined ONLINE_JUDGE
  94. SetFile(zjoi15_substring);
  95. #endif
  96.  
  97. #ifdef UseFREAD
  98. fread(file_ptr, 1, FreadLenth, stdin);
  99. #endif
  100.  
  101. init();
  102. LL Ans = 0;
  103. for(SAM *it = Node + 1;it != iter;++it)Ans += it->len - it->link->len;
  104. printf("%lld\n", Ans);
  105.  
  106. #ifdef DEBUG
  107. printf("\n%.3lf sec \n", (double)clock() / CLOCKS_PER_SEC);
  108. #endif
  109. return 0;
  110. }
  111.