結果

問題 No.2626 Similar But Different Name
ユーザー TairitsuMeowTairitsuMeow
提出日時 2024-02-09 22:51:48
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,130 ms / 3,000 ms
コード長 2,785 bytes
コンパイル時間 1,488 ms
コンパイル使用メモリ 174,852 KB
実行使用メモリ 358,440 KB
最終ジャッジ日時 2024-09-28 16:03:36
合計ジャッジ時間 25,615 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 115 ms
326,284 KB
testcase_01 AC 116 ms
324,284 KB
testcase_02 AC 116 ms
324,320 KB
testcase_03 AC 116 ms
324,272 KB
testcase_04 AC 116 ms
324,400 KB
testcase_05 AC 116 ms
324,568 KB
testcase_06 AC 116 ms
324,324 KB
testcase_07 AC 116 ms
326,388 KB
testcase_08 AC 118 ms
326,428 KB
testcase_09 AC 118 ms
324,624 KB
testcase_10 AC 119 ms
330,508 KB
testcase_11 AC 117 ms
330,636 KB
testcase_12 AC 115 ms
330,480 KB
testcase_13 AC 117 ms
330,636 KB
testcase_14 AC 114 ms
330,632 KB
testcase_15 AC 118 ms
330,504 KB
testcase_16 AC 119 ms
330,632 KB
testcase_17 AC 118 ms
330,636 KB
testcase_18 AC 1,112 ms
358,440 KB
testcase_19 AC 558 ms
352,824 KB
testcase_20 AC 562 ms
349,532 KB
testcase_21 AC 573 ms
349,400 KB
testcase_22 AC 1,113 ms
355,548 KB
testcase_23 AC 1,130 ms
355,552 KB
testcase_24 AC 1,111 ms
355,296 KB
testcase_25 AC 1,112 ms
353,888 KB
testcase_26 AC 1,118 ms
355,540 KB
testcase_27 AC 1,120 ms
358,272 KB
testcase_28 AC 1,108 ms
355,660 KB
testcase_29 AC 1,098 ms
355,744 KB
testcase_30 AC 1,102 ms
354,876 KB
testcase_31 AC 1,115 ms
356,976 KB
testcase_32 AC 1,123 ms
353,392 KB
testcase_33 AC 1,111 ms
353,888 KB
testcase_34 AC 1,106 ms
355,560 KB
testcase_35 AC 1,123 ms
354,868 KB
testcase_36 AC 1,109 ms
352,652 KB
testcase_37 AC 1,121 ms
353,400 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

// Problem: No.2626 Similar But Different Name
// Contest: yukicoder
// URL: https://yukicoder.me/problems/no/2626
// Memory Limit: 512 MB
// Time Limit: 3000 ms

#include<bits/stdc++.h>
#define debug(x) cerr<<(#x)<<" "<<(x)<<endl
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
#define pii pair<ll,ll>
#define rep(i,a,b) for(ll i=(a);i<=(b);++i)
#define per(i,a,b) for(ll i=(a);i>=(b);--i)
using namespace std;
bool Mbe; 
ll read(){
	ll x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
void write(ll x){
	if(x<0)putchar('-'),x=-x;
	if(x>9)write(x/10);
	putchar(x%10+'0');
}
const ll N=4e6+9;
const ull B=131;
ll n,m,k,ans[N],r[N],f[N],g[N];
char s[N],t[N];
ull hs[N],ht[N],pw[N],hs2[N];
struct Comp{
	ld x,y;
	Comp(ld a=0,ld b=0){x=a,y=b;}
};
inline Comp operator+(const Comp&a,const Comp&b){return Comp(a.x+b.x,a.y+b.y);}
inline Comp operator-(const Comp&a,const Comp&b){return Comp(a.x-b.x,a.y-b.y);}
inline Comp operator*(const Comp&a,const Comp&b){return Comp(a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x);}
Comp a[N],b[N]; 
const ld PI=acos(-1.0);
void FFT(Comp*a,ll p,ll sgn){
	rep(i,0,p-1){
		if(i<r[i])swap(a[i],a[r[i]]);
	}
	for(ll mid=2;mid<=p;mid<<=1){
		Comp wn=Comp(cos(2*PI/mid),sin(2*PI/mid)*sgn);
		ll R=mid>>1;
		for(ll j=0;j<p;j+=mid){
			Comp w=Comp(1.0,0.0);
			rep(k,j,j+R-1){
				Comp x=w*a[R+k];
				a[R+k]=a[k]-x,a[k]=a[k]+x;
				w=w*wn;
			}
		}
	}
}
void MUL(ll*f,ll*g,ll n,ll m){
	ll c=0,p=1;
	while(p<=n+m)p*=2,c++;
	rep(i,0,p-1)r[i]=(r[i>>1]>>1)|((i&1)?(p>>1):0);
	rep(i,0,n-1)a[i].x=f[i],a[i].y=0;
	rep(i,n,p-1)a[i].x=0,a[i].y=0;
	rep(i,0,m-1)b[i].x=g[i],b[i].y=0;
	rep(i,m,p-1)b[i].x=0,b[i].y=0;
	FFT(a,p,1);
	FFT(b,p,1);
	rep(i,0,p-1)a[i]=a[i]*b[i];
	FFT(a,p,-1);
	rep(i,0,n+m-1)f[i]=(ll)(a[i].x/p+0.5);
}
ull Hs(ll l,ll r){
	return hs[r]-hs[l-1]*pw[r-l+1];
}
ull Hs2(ll l,ll r){
	return hs2[r]-hs2[l-1]*pw[r-l+1];
}
map<ull,ll>mp;
bool Med;
int main(){
	cerr<<fabs(&Med-&Mbe)/1048576.0<<"MB\n";
	n=read(),m=read(),k=read();
	scanf("%s%s",s,t);
	rep(i,0,n-1)f[i]=islower(s[i])?1:0;
	rep(i,0,m-1)g[i]=islower(t[i])?1:0;
	reverse(g,g+m);
	MUL(f,g,n,m);
	rep(i,m-1,n-1)ans[i]+=f[i];
	memset(f,0,sizeof(f));
	memset(g,0,sizeof(g));
	rep(i,0,n-1)f[i]=isupper(s[i])?1:0;
	rep(i,0,m-1)g[i]=isupper(t[i])?1:0;
	reverse(g,g+m);
	MUL(f,g,n,m);
	pw[0]=1;
	rep(i,1,n)pw[i]=pw[i-1]*B;
	rep(i,m-1,n-1)ans[i]+=f[i];
	hs[0]=tolower(s[0]);
	rep(i,1,n-1)hs[i]=hs[i-1]*B+tolower(s[i]);
	hs2[0]=(islower(s[0]));
	rep(i,1,n-1)hs2[i]=hs2[i-1]*B+(islower(s[i]));
	ht[0]=tolower(t[0]);
	rep(i,1,m-1)ht[i]=ht[i-1]*B+tolower(t[i]);
	ll res=0;
	rep(i,m-1,n-1){
		if(Hs(i-m+1,i)!=ht[m-1])continue;
		if(ans[i]!=m&&ans[i]+k>=m)res++;
	}
	write(res);
	return 0;
}
0