Entries

スポンサーサイト (この記事を編集する[管理者用])

上記の広告は1ヶ月以上更新のないブログに表示されています。
新しい記事を書く事で広告が消せます。

コメント

コメントの投稿

コメントの投稿
管理者にだけ表示を許可する

トラックバック

トラックバック URL
http://rsujskf.blog32.fc2.com/tb.php/1111-f372946f
この記事にトラックバックする(FC2ブログユーザー)

Beta Round #7 D問題 - Palindrome Degree (この記事を編集する[管理者用])

Source

Codeforces Beta Round #7 D問題
Problem description
Beta Round #7の自分の参加記録

問題概要

回文Aの次数を,Aの前半と後半が次数k-1の回文だったら次数kとする.前半と後半が回文じゃなければ次数1.
与えられた文字列の全てのprefixの次数の和を求める問題.
文字列の長さは5000000以下.

解法

全てのプレフィックスが回文かどうかの判定ができれば後はDPするだけ.
判定はrolling hashを使えばできる.
Manacherとかでもできるらしいけど,自分はよく知らない.

C言語のスパゲッティなコード

このコードは実際に本番でサブミットしたコードです.見にくいと思われます.

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define REP(i,a,b) for(i=a;i<b;i++)
#define rep(i,n) REP(i,0,n)

#define ull unsigned long long
#define D 12009

ull  pw[6666666], hs1[6666666], hs2[6666666];
char in[6666666];
int  dp[6666666];

int main(){
  int i,j,k,l,m,n;
  int a,b,len;
  long long res=0;

  pw[0]=1; REP(i,1,6666666) pw[i]=pw[i-1]*D;
  scanf("%s",in);
  for(i=0;;i++) if(in[i]<' ') break; n=i;

  hs1[0]=0; rep(i,n) hs1[i+1]=hs1[i]*D+in[i];
  hs2[0]=0; rep(i,n) hs2[i+1]=hs2[i]*D+in[n-1-i];

  dp[0]=0;
  REP(i,1,n+1){
    len = i/2;
    a=0; b=i-len;
    if( hs1[a+len]-hs1[a]*pw[len] != hs2[n-b]-hs2[n-b-len]*pw[len] ){dp[i]=0; continue;}
    dp[i] = dp[len]+1;
    res += dp[i];
  }

  printf("%I64lld\n",res);
  return 0;
}

コメント

コメントの投稿

コメントの投稿
管理者にだけ表示を許可する

トラックバック

トラックバック URL
http://rsujskf.blog32.fc2.com/tb.php/1111-f372946f
この記事にトラックバックする(FC2ブログユーザー)

Appendix

Recent Articles

ブログ内検索

Ads


(プライバシーポリシー)
上記広告は1ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。