// Slower C++ solution:
#include <bits/stdc++.h>

using namespace std;

typedef long long   signed int LL;

typedef long long unsigned int LU;

#define incID(i, l, r) for(LL i = (l)    ; i <  (r); ++i)

#define incII(i, l, r) for(LL i = (l)    ; i <= (r); ++i)

#define decID(i, l, r) for(LL i = (r) - 1; i >= (l); --i)

#define decII(i, l, r) for(LL i = (r)    ; i >= (l); --i)

#define inc(i, n)  incID(i, 0, n)

#define inc1(i, n) incII(i, 1, n)

#define dec(i, n)  decID(i, 0, n)

#define dec1(i, n) decII(i, 1, n)

#define inID(v, l, r) ((l) <= (v) && (v) <  (r))

#define inII(v, l, r) ((l) <= (v) && (v) <= (r))

#define PB push_back

#define EB emplace_back

#define MP make_pair

#define FI first

#define SE second

#define  ALL(v)  v.begin(),  v.end()

#define RALL(v) v.rbegin(), v.rend()

template<typename T> bool setmin  (T & a, T b) { if(b <  a) { a = b; return true; } else { return false; } }

template<typename T> bool setmax  (T & a, T b) { if(b >  a) { a = b; return true; } else { return false; } }

template<typename T> bool setmineq(T & a, T b) { if(b <= a) { a = b; return true; } else { return false; } }

template<typename T> bool setmaxeq(T & a, T b) { if(b >= a) { a = b; return true; } else { return false; } }

LL mo(LL a, LL b) { assert(b > 0); a %= b; if(a < 0) { a += b; } return a; }

LL fl(LL a, LL b) { assert(b > 0); return (a > 0 ? a / b : (a - b + 1) / b); }

LL ce(LL a, LL b) { assert(b > 0); return (a < 0 ? a / b : (a + b - 1) / b); }

template<typename T> T gcd(T a, T b) { return (b == 0 ? a : gcd(b, a % b)); }

template<typename T> T lcm(T a, T b) { return a / gcd(a, b) * b; }

#define bit(b, i) (((b) >> (i)) & 1)

#define BC __builtin_popcountll

#define SC static_cast

#define SI(v) SC<int>(v.size())

#define SL(v) SC<LL >(v.size())

#define RF(e, v) for(auto & e: v)

#define ef else if

#define UR assert(false)



// ---- ----



template<LL M> class ModInt {

private:

	LL v = 0;

public:

	ModInt() { }

	ModInt(LL vv) { setval(vv); }

	ModInt & setval(LL vv) { v = vv % M; if(v < 0) { v += M; } return (*this); }

	LL getval() const { return v; }

	ModInt & operator+=(const ModInt & b)       { return setval(v + b.v); }

	ModInt & operator-=(const ModInt & b)       { return setval(v - b.v); }

	ModInt & operator*=(const ModInt & b)       { return setval(v * b.v); }

	ModInt & operator/=(const ModInt & b)       { return setval(v * b.inv()); }

	ModInt & operator^=(            LU b)       { return setval(ex(v, b)); }

	ModInt   operator+ (                ) const { return ModInt(+v); }

	ModInt   operator- (                ) const { return ModInt(-v); }

	ModInt   operator+ (const ModInt & b) const { return ModInt(v + b.v); }

	ModInt   operator- (const ModInt & b) const { return ModInt(v - b.v); }

	ModInt   operator* (const ModInt & b) const { return ModInt(v * b.v); }

	ModInt   operator/ (const ModInt & b) const { return ModInt(v * b.inv()); }

	ModInt   operator^ (            LU b) const { return ModInt(ex(v, b)); }

	LL inv() const {

		LL x = (ex_gcd(v, M).FI + M) % M;

		assert(v * x % M == 1);

		return x;

	}

	LL ex(LL a, LU b) const {

		LL D = 64, x[64], y = 1;

		inc(i, D) { if((b >> i) == 0) { D = i; break; } }

		inc(i, D) { x[i] = (i == 0 ? a : x[i - 1] * x[i - 1]) % M; }

		inc(i, D) { if((b >> i) & 1) { (y *= x[i]) %= M; } }

		return y;

	}

	pair<LL, LL> ex_gcd(LL a, LL b) const {

		if(b == 0) { return MP(1, 0); }

		auto p = ex_gcd(b, a % b);

		return MP(p.SE, p.FI - (a / b) * p.SE);

	}

};

template<LL M> ModInt<M> operator+(LL a, const ModInt<M> & b) { return  b + a; }

template<LL M> ModInt<M> operator-(LL a, const ModInt<M> & b) { return -b + a; }

template<LL M> ModInt<M> operator*(LL a, const ModInt<M> & b) { return  b * a; }

template<LL M> ModInt<M> operator/(LL a, const ModInt<M> & b) { return  a * b.inv(); }

template<LL M> istream & operator>>(istream & is, ModInt<M> & b) { LL v; is >> v; b.setval(v); return is; }

template<LL M> ostream & operator<<(ostream & os, const ModInt<M> & b) { return (os << b.getval()); }



// ---- ----



typedef ModInt<1'000'000'007> MI;



int n;

string s;



const int K = 2;

MI h[K][5001], b[K] = { 100000000, 12345678 }; //, 87654321, 7777777 };

MI e[K][5000];



bool eq(int i, int j, int v) {

	inc(k, K) {

		MI x = (h[k][i + v] - h[k][i]) * e[k][j - i];

		MI y = (h[k][j + v] - h[k][j]);

		if(x.getval() != y.getval()) { return false; }

	}

	return true;

}



bool f(int v) {

	incID(j, v, n) {

	incII(i, 0, j - v) {

		if(! (i + v <= j && j + v <= n)) { continue; }

		if(eq(i, j, v)) { return true; }

	}

	}

	

	return false;

}



int main() {

	cin >> n >> s;

	

	inc(k, K) {

	inc(i, n) {

		h[k][i + 1] = h[k][i] + (s[i] * (b[k] ^ i));

		e[k][i] = (b[k] ^ i);

	}

	}

	

	int L = 0, H = n / 2 + 1;

	while(H - L > 1) {

		int M = (L + H) / 2;

		(f(M) ? L : H) = M;

	}

	

	cout << L << endl;

	

	return 0;

}


// Faster C++ solution for the same problem (measured speedup: 1.1609524205507287x):
#include <bits/stdc++.h>

using namespace std;

typedef long long   signed int LL;

typedef long long unsigned int LU;

#define incID(i, l, r) for(LL i = (l)    ; i <  (r); ++i)

#define incII(i, l, r) for(LL i = (l)    ; i <= (r); ++i)

#define decID(i, l, r) for(LL i = (r) - 1; i >= (l); --i)

#define decII(i, l, r) for(LL i = (r)    ; i >= (l); --i)

#define inc(i, n)  incID(i, 0, n)

#define inc1(i, n) incII(i, 1, n)

#define dec(i, n)  decID(i, 0, n)

#define dec1(i, n) decII(i, 1, n)

#define inID(v, l, r) ((l) <= (v) && (v) <  (r))

#define inII(v, l, r) ((l) <= (v) && (v) <= (r))

#define PB push_back

#define EB emplace_back

#define MP make_pair

#define FI first

#define SE second

#define  ALL(v)  v.begin(),  v.end()

#define RALL(v) v.rbegin(), v.rend()

template<typename T> bool setmin  (T & a, T b) { if(b <  a) { a = b; return true; } else { return false; } }

template<typename T> bool setmax  (T & a, T b) { if(b >  a) { a = b; return true; } else { return false; } }

template<typename T> bool setmineq(T & a, T b) { if(b <= a) { a = b; return true; } else { return false; } }

template<typename T> bool setmaxeq(T & a, T b) { if(b >= a) { a = b; return true; } else { return false; } }

LL mo(LL a, LL b) { assert(b > 0); a %= b; if(a < 0) { a += b; } return a; }

LL fl(LL a, LL b) { assert(b > 0); return (a > 0 ? a / b : (a - b + 1) / b); }

LL ce(LL a, LL b) { assert(b > 0); return (a < 0 ? a / b : (a + b - 1) / b); }

template<typename T> T gcd(T a, T b) { return (b == 0 ? a : gcd(b, a % b)); }

template<typename T> T lcm(T a, T b) { return a / gcd(a, b) * b; }

#define bit(b, i) (((b) >> (i)) & 1)

#define BC __builtin_popcountll

#define SC static_cast

#define SI(v) SC<int>(v.size())

#define SL(v) SC<LL >(v.size())

#define RF(e, v) for(auto & e: v)

#define ef else if

#define UR assert(false)



// ---- ----



int n, dp[5001][5001];

string s;



int main() {

	cin >> n >> s;

	

	dec(i, n) {

	dec(j, n) {

		if(s[i] == s[j]) { dp[i][j] = dp[i + 1][j + 1] + 1; }

	}

	}

	

	int ans = 0;

	inc(i, n) {

	inc(j, i) {

		setmax(ans, min(dp[i][j], SC<int>(i - j)));

	}

	}

	

	cout << ans << endl;

	

	return 0;

}