結果

問題 No.1916 Making Palindrome on Gird
コンテスト
ユーザー C
提出日時 2026-11-05 23:48:51
言語 C++23(gcc16)
(gcc 16.2.0 + boost 1.92.0 + ACL)
コンパイル:
g++-16 -O2 -fplugin=/usr/lib/yukicoder/noavx512-gcc16.so -lm -std=c++23 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 67 ms / 3,000 ms
+ 371µs
コード長 1,659 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 7,715 ms
コンパイル使用メモリ 396,784 KB
実行使用メモリ 9,924 KB
最終ジャッジ日時 2026-11-05 23:49:03
合計ジャッジ時間 10,921 ms
ジャッジサーバーID
（参考情報） 		judge2_0 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
#include <atcoder/all>
#include <boost/polygon/point_data.hpp>
#include <boost/polygon/voronoi.hpp>
#define rep(i,a,b) for(int i=(a);i<(b);i++)
#define rrep(i,a,b) for(int i=(b)-1;i>=(a);i--)
#define all(x) begin(x),end(x)
#define sz(x) (int)(x).size()
using namespace std;
using namespace atcoder;
using ll=long long;
using ull=unsigned long long;
using ld=long double;
using vll=vector<ll>;
using vvll=vector<vll>;
using vi=vector<int>;
using vvi=vector<vi>;
using pll=pair<ll,ll>;
using pi=pair<int,int>;

// using mint=modint998244353;
using mint=modint1000000007;
// using mint2=static_modint<2>;

// using BoostPoint=boost::polygon::point_data<int>;
// using Voronoi=boost::polygon::voronoi_diagram<double>;

// using mint=LazyMontgomeryModInt<998244353>;
using vm=vector<mint>;
using vvm=vector<vm>;

int main()
{
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
	//斜めi+j=tを見て素直にdpすれば良い
	int H,W;
	cin>>H>>W;
	vector<string>S(H);
	for(auto&row:S)cin>>row;
	if(H>W){
		vector<string>nS(W,string(H,' '));
		rep(i,0,H)rep(j,0,W)nS[j][i]=S[i][j];
		swap(H,W);
		swap(nS,S);
	}
	int D=H+W-2,m=D/2;
	vvm dp(H,vm(H));
	if(S[0][0]==S[H-1][W-1])dp[0][H-1]=1;
	rep(t,0,m){
		vvm ndp(H,vm(H));
		rep(a,0,H){
			rep(b,0,H){
				if(dp[a][b]==0)continue;
				rep(u,0,2){
					rep(v,0,2){
						int na=a+u,nb=b-v;
						int x=t+1-na,y=D-t-1-nb;
						if(na>=H||nb<0||x<0||x>=W||y<0||y>=W)continue;
						if(S[na][x]!=S[nb][y])continue;
						ndp[na][nb]+=dp[a][b];
					}
				}
			}
		}
		swap(dp,ndp);
	}
	mint ans=0;
	rep(a,0,H){
		ans+=dp[a][a];
		if(D%2&&a+1<H)ans+=dp[a][a+1];
	}
	cout<<ans.val()<<'\n';
}
0