#include<bits/stdc++.h>
using namespace std;
using ll=long long;
const int N=2e5+10;
int n;
char a[2][N];
struct SA{
    static const int N=::N*2,K=__lg(N)+1;
    char a[N];
    int sa[N],rk[N],h[K][N];
    void getsa(int n,int m=127){
        static int old[N*2],id[N],p[N],cnt[N];
        fill(old,old+1+n+n,0);
        for(int i=1;i<=n;i++)cnt[rk[i]=a[i]]++;
        for(int i=1;i<=m;i++)cnt[i]+=cnt[i-1];
        for(int i=n;i>=1;i--)sa[cnt[rk[i]]--]=i;
        fill(cnt,cnt+1+m,0);
        for(int len=1,k;len==1||m^n;len<<=1){
            k=0;
            for(int i=n-len+1;i<=n;i++)p[++k]=i;
            for(int i=1;i<=n;i++)if(sa[i]>len)p[++k]=sa[i]-len;
            for(int i=1;i<=n;i++)cnt[id[i]=rk[p[i]]]++;
            for(int i=1;i<=m;i++)cnt[i]+=cnt[i-1];
            for(int i=n;i>=1;i--)sa[cnt[id[i]]--]=p[i];
            copy(rk,rk+1+n,old),fill(cnt,cnt+1+m,0),m=0;
            for(int i=1;i<=n;i++){
                rk[sa[i]]=old[sa[i]]==old[sa[i-1]]&&old[sa[i]+len]==old[sa[i-1]+len]?m:++m;
            }
        }
        for(int i=1,x=0;i<=n;i++){
            for(x-=!!x;max(i,sa[rk[i]-1])+x<=n&&a[i+x]==a[sa[rk[i]-1]+x];x++);
            h[0][rk[i]]=x;
        }
        for(int i=1;(1<<i)<=n;i++){
            for(int j=1;j+(1<<i)-1<=n;j++){
                h[i][j]=min(h[i-1][j],h[i-1][j+(1<<i-1)]);
            }
        }
    }
    int LCP(int l,int r){
        if(l=rk[l],r=rk[r],l>r)swap(l,r);
        int k=__lg(r-l++);
        return min(h[k][l],h[k][r-(1<<k)+1]);
    }
}A,B;
int LCP(int a,int b,int x,int y){
    if(min(b,y)<1||max(b,y)>n)return 0;
    return min({A.LCP(a*n+b,x*n+y),n-b+1,n-y+1});
}
int LCS(int a,int b,int x,int y){
    if(min(b,y)<1||max(b,y)>n)return 0;
    return min({B.LCP(n+n+1-(a*n+b),n+n+1-(x*n+y)),b,y});
}
int solve(int op){
    for(int i=1;i<=n;i++){
        A.a[i]=a[0][i],A.a[i+n]=a[1][i];
    }
    for(int i=1;i<=n+n;i++)B.a[i]=A.a[n+n-i+1];
    A.getsa(n+n),B.getsa(n+n);
    for(int k=(n+1)/2;k>=1;k--){
        for(int l=0,r=k;r<=n;l+=k,r+=k){
            int x=LCS(0,l,0,r);
            int y=LCP(0,l+1,0,r+1);
            int z=LCP(0,l+1+y,1,r+y);
            if(x+y+z>=k)return k;
        }
        for(int l=0,r=k;r<=n;l+=k,r+=k){
            int x=LCP(0,l+1,1,r);
            int y=LCS(0,l,1,r-1);
            int z=LCS(0,l-y,0,r-y);
            // if(op&&k==3&&l==3){
            //     fprintf(stderr,"%d %d %d %d %d %s %s %d %d\n",l,r,x,y,z,a[0]+1,a[1]+1,l-y,r-y);
            //     fprintf(stderr,"%s\n",B.a+1);
            // }
            if(x+y+z>=k)return k;
        }
    }
    return 0;
}
int main(){
    scanf("%d%s%s",&n,a[0]+1,a[1]+1);
    int ans=solve(0);
    swap(a[0],a[1]);
    reverse(a[0]+1,a[0]+1+n);
    reverse(a[1]+1,a[1]+1+n);
    ans=max(ans,solve(1));
    cout<<ans*2<<endl;
    return 0;
}