#pragma GCC target("avx2") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #include using namespace std; #define rep(i,n) for(int i = 0; i < int(n); ++i) #define rrep(i,n) for(int i = int(n)-1; i >= 0; --i) #define all(x) (x).begin(), (x).end() #define rall(x) (x).rbegin(), (x).rend() template void chmax(T& a, const T& b) {a = max(a, b);} template void chmin(T& a, const T& b) {a = min(a, b);} using ll = long long; using P = pair; using VI = vector; using VVI = vector; using VL = vector; using VVL = vector; using ull = unsigned long long; ull e[2000][32], et[2000][32], now[2000][32], tmp[2000][32]; int main() { ios::sync_with_stdio(false); cin.tie(0); int v; ull d; cin >> v >> d; assert(v <= 2000); rep(i, v) rep(j, v) { char t; cin >> t; e[i][j >> 6] |= ull(t - '0') << (j & 63); } rep(i, v) { now[i][i >> 6] = 1 << (i & 63); } while(d) { rep(i, 2000) rep(j, 32) et[i][j] = 0; rep(i, 2000) rep(j, 2000) { et[j][i >> 6] |= (e[i][j >> 6] >> (j & 63) & 1) << (i & 63); } if (d & 1) { rep(i, 2000) rep(j, 32) tmp[i][j] = now[i][j], now[i][j] = 0; rep(i, 2000) rep(j, 2000) { /* rep(k, 32) { if (tmp[i][k] & et[j][k]) { now[i][j >> 6] |= 1 << (j & 63); break; } } */ ull x = 0; rep(k, 32) x |= tmp[i][k] & et[j][k]; now[i][j >> 6] |= bool(x) << (j & 63); } } rep(i, 2000) rep(j, 32) tmp[i][j] = e[i][j], e[i][j] = 0; rep(i, 2000) rep(j, 2000) { /* rep(k, 32) { if (tmp[i][k] & et[j][k]) { e[i][j >> 6] |= 1 << (j & 63); break; } } */ ull x = 0; rep(k, 32) x |= tmp[i][k] & et[j][k]; e[i][j >> 6] |= bool(x) << (j & 63); } d >>= 1; } bool ok = true; rep(i, v) rep(j, v) if (!(now[i][j >> 6] >> (j & 63) & 1)) ok = false; cout << (ok ? "Yes\n" : "No\n"); }