#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long ll; typedef unsigned int ui; const ll mod = 1000000007; const ll INF = (ll)1000000007 * 1000000007; typedef pair P; #define stop char nyaa;cin>>nyaa; #define rep(i,n) for(int i=0;i=0;i--) #define Rep(i,sta,n) for(int i=sta;i=sta;i--) #define rep1(i,n) for(int i=1;i<=n;i++) #define per1(i,n) for(int i=n;i>=1;i--) #define Rep1(i,sta,n) for(int i=sta;i<=n;i++) typedef long double ld; typedef complex Point; const ld eps = 1e-8; const ld pi = acos(-1.0); typedef pair LP; struct uf { private: vector par,ran; public: vector edge,size; void init(int n) { par.resize(n,0); ran.resize(n,0); size.resize(n,1); edge.resize(n,0); rep(i,n) { par[i]=i; } } int find(int x) { if (par[x]==x)return x; else return par[x]=find(par[x]); } void unite(int x, int y) { x=find(x),y=find(y); if (x==y){ edge[x]++; return; } if (ran[x]> w >> h; rep(i,h){ rep(j,w){ cin >> m[i][j]; } } uf u;u.init(h*w); rep(i,h){ rep(j,w-1){ if (m[i][j]==m[i][j+1]) u.unite(i*w+j,i*w+j+1); } } rep(j,w){ rep(i,h-1){ if (m[i][j]==m[i+1][j]) u.unite(i*w+j,(i+1)*w+j); } } rep(i,h){ rep(j,w){ //cout << i << " " << j << endl; //cout << "size:" << u.size[u.find(i*w+j)] << " edge:" << u.edge[u.find(i*w+j)] << endl; if (u.size[u.find(i*w+j)]<=u.edge[u.find(i*w+j)]) { cout << "possible" << endl; return; } } } cout << "impossible" << endl; } int main(){ ios::sync_with_stdio(false); cin.tie(0); cout << fixed << setprecision(50); solve(); }