1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
#include <iostream>
#include <vector>
#include <algorithm>

using namespace std;

struct pkt
{
    int x, y;
};

int dx[] = {-1, 0, 1, 0};
int dy[] = {0, 1, 0, -1};
long long odp = 0, n, k;
const long long MOD = 1000000007;

bool Obwodowy(vector<string>& p, pkt a)
{
    for (int i = 0; i < 4; ++i)
    {
        pkt nowy = {a.x + dx[i], a.y + dy[i]};
        if (nowy.x >= 0 && nowy.x < n && nowy.y >= 0 && nowy.y < n)
            if (p[nowy.x][nowy.y] == '#')
                return true;
    }
    return false;
}

bool Ok(vector<int> a, vector<pkt>& puste, vector<string>& p)
{
    sort(a.begin(), a.end());
    do
    {
        bool ok = false;
        for (int i = 0; i < k; ++i)
        {
            bool przylega = false;
            for (int j = 0; j < 4; ++j)
            {
                pkt nowy = {puste[a[i]].x + dx[j], puste[a[i]].y + dy[j]};
                if (nowy.x >= 0 && nowy.x < n && nowy.y >= 0 && nowy.y < n)
                    if (p[nowy.x][nowy.y] == '#')
                        przylega = true;
            }
            if (przylega)
                p[puste[a[i]].x][puste[a[i]].y] = '#';
            else
                break;
            if (i == k - 1)
                ok = true;
        }
        for (int i = 0; i < k; ++i)
            p[puste[a[i]].x][puste[a[i]].y] = '.';
        if (ok)
            return true;
    } while (next_permutation(a.begin(), a.end()));
    return false;
}

void rec(vector<int> a, int nr, int x, vector<pkt>& puste, vector<string>& p)
{
    a[nr] = x;
    if (nr == k - 1)
    {
        if (Ok(a, puste, p))
            ++odp;
        return;
    }
    for (int i = x + 1; i < puste.size(); ++i)
        rec(a, nr + 1, i, puste, p);
}

int main()
{
    ios_base::sync_with_stdio(0);
    cin >> n >> k;
    vector<string> p(n);
    vector<pkt> puste;
    for (int i = 0; i < n; ++i)
    {
        cin >> p[i];
        for (int j = 0; j < n; ++j)
            if (p[i][j] == '.')
                puste.push_back({i, j});
    }
    int obw = 0;
    for (int i = 0; i < n; ++i)
        for (int j = 0; j < n; ++j)
            if (p[i][j] == '.' && Obwodowy(p, {i, j}))
                ++obw;
    if (k == 1)
    {
        cout << obw;
        return 0;
    }
    if (k == 2)
    {
        long long x = 0, y = 0;
        for (int i = 0; i < n; ++i)
        {
            for (int j = 0; j < n; ++j)
            {
                if (p[i][j] == '.' && Obwodowy(p, {i, j}))
                {
                    long long nowe = 0;
                    for (int k = 0; k < 4; ++k)
                    {
                        pkt nowy = {i + dx[k], j + dy[k]};
                        if (nowy.x >= 0 && nowy.x < n && nowy.y >= 0 && nowy.y < n)
                            if (p[nowy.x][nowy.y] == '.' && !Obwodowy(p, nowy))
                                ++nowe;
                    }
                    x += obw - 1;
                    y += nowe;
                }
            }
        }
        cout << (x / 2 + y) % MOD;
        return 0;
    }
    vector<int> a(k);
    for (int i = 0; i < puste.size(); ++i)
        rec(a, 0, i, puste, p);
    cout << odp;
    return 0;
}