#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<map>
#include<vector>
#include<queue>
#include<functional>
#include<cstring>
#include<cstdlib>
#include<complex>
using namespace std;
typedef long long ll;
typedef complex < double > Point;
/*
    Pointの使い方
    Point AにおけるX座標...real(A) or A.real()
                   Y座標...imag(A) or A.imag() 状況によってちょいちょい使い分けとか

    Point A,BにおいてA+B,A-B,(double)k * Aは定義済み
    原点からの距離はsqrt(norm(A))
    二点間の距離はsqrt(norm(A-B))
    norm(A)はX^2+Y^2を返す
*/
typedef vector < Point > Polygon;
namespace std {
    bool operator < (Point a,Point b) {
        //ここで座標の比較における"大きい"を定義
        return real(a) != real(b) ? real(a) < real(b) : imag(a) < imag(b);
    }
}

const int INF = 10000000;



struct Line : Polygon{
    Line(Point p,Point q){ push_back(p);push_back(q); }
};
struct Circle{
    Point p;
    double r;
    Circle(Point a, double s){ p = a;r = s; }
};

double cross(Point p,Point q){
    return imag(conj(p)*q);
}

double dot(Point p,Point q){
    return real(conj(p)*q);
}

int ccw(Point a,Point b,Point c){
    b-=a,c-=a;
    if(cross(b,c)>0)return 1; //反時計回り
    if(cross(b,c)<0)return -1; //時計回り
    if(dot(b,c)<0)return 2; //直線上に c - a - b
    if(norm(b)<norm(c))return -2; // 直線上に a - b - c
    return 0; // 直線上に a - c - b
}

Point normalize(Point P, double l){
    double a = real(P), b = imag(P);
    double y = sqrt(b*b/(a*a+b*b)*l*l);
    double x = (a/b)*y;

    return Point(x,y);
}

double heron(Point A, Point B, Point C){
    double a = sqrt(norm(A-B));
    double b = sqrt(norm(B-C));
    double c = sqrt(norm(C-A));
    double s = (a+b+c)/2;

    return sqrt(s*(s-a)*(s-b)*(s-c));
}

//↓2円の交点二つを返す（ポインタ）
void intersectCC(Circle A, Circle B, Point *p, Point *q){

    Point P,Q;
    if(B.p<A.p)swap(A.p,B.p);
    Point dv = B.p - A.p;//ベクトル


    double L = sqrt(norm(dv));
    //距離がLの場合
    if( L==(A.r + B.r) ){
        //交点1つ
        //虚数軸に沿って並んでいるとバグる？
        dv = normalize(dv,A.r);
        P = Q = Point(real(dv)+real(A.p),imag(dv)+imag(A.p));
    }
    else if( L<(A.r + B.r) && L > fabs(A.r - B.r) ){
        //交点2つ
        /*
        交点の一方Qからベクトルdvに垂線を下ろしたところをptとする
        ふたつの直角三角形p0-pt-Q、p1-pt-Qができる
        p0-pt間の長さをtとすると、これらの三角形のQ-ptの長さは等しいため、

        r0*r0 - t*t = r1*r1 - (L-t)*(L-t)
        r0*r0 - t*t = r1*r1 - L*L - t*t + 2*L*t
        t = (r0*r0 - r1*r1 + L*L)/(2*L)

        tから、p0-Qベクトルとdvベクトルとの角度aを求め、
        dvベクトルから、+a,-aの座標を求めれば交点となる
        */

        //tを求める
        double t = 0.5*( A.r*A.r - B.r*B.r + L*L)/L;

        //p0-Qベクトルとdvベクトルの角度
        double a = acos( t/A.r );

        //dvベクトルの角度
        double a0 = atan2( imag(dv), real(dv) );

        //交点
        P = Point( real(A.p) + A.r*cos(a0+a), imag(A.p) + A.r*sin(a0+a) );
        Q = Point( real(A.p) + A.r*cos(a0-a), imag(A.p) + A.r*sin(a0-a) );

    }
    else{
        //交点なし
    }

    *p = P;
    *q = Q;

    return;
}

int main(){
    int n;
    while(cin>>n,n){
        vector < Point > crosspoint;
        Point p[100];
        for(int i = 0;i < n;i++){
            scanf("%lf,%lf",&real(p[i]),&imag(p[i]));
        }
        for(int i = 0;i < n;i++){
            for(int j = i + 1;j < n;j++){
                Point P,Q;
                intersectCC(Circle(p[i],1),Circle(p[j],1),&P,&Q);
                crosspoint.push_back(P);
                if(P!=Q)crosspoint.push_back(Q);
            }
        }
        int ans = 1;
        for(int i = 0;i < crosspoint.size();i++){
            int cnt = 0;
            for(int j = 0;j < n;j++){
                if(sqrt(norm(crosspoint[i] - p[j]))<=1.0)cnt++;
            }
            if(cnt > ans)ans = cnt;
        }
        cout<<ans<<endl;
    }
    return 0;
}