#include "bits/stdc++.h"
using namespace std;
using ll = long long;
const ll INF = 0x3f3f3f3f3f3f3f3f;
const int N = 2e5+10;

struct MF {
    struct edge {
        int v, nxt, cap, flow;
    } e[N];

    int fir[N], cnt = 0;

    int n, S, T;
    ll maxflow = 0;
    int dep[N], cur[N];

    void init(int x) {
        memset(fir, -1, sizeof fir);
        maxflow = 0;
        cnt = 0;
    }

    void addedge(int u, int v, int w) {
        e[cnt] = {v, fir[u], w, 0};
        fir[u] = cnt++;
        e[cnt] = {u, fir[v], 0, 0};
        fir[v] = cnt++;
    }

    bool bfs() {
        queue<int> q;
        memset(dep, 0, sizeof(int) * (n + 1));

        dep[S] = 1;
        q.push(S);
        while (q.size()) {
            int u = q.front();
            q.pop();
            for (int i = fir[u]; ~i; i = e[i].nxt) {
                int v = e[i].v;
                if ((!dep[v]) && (e[i].cap > e[i].flow)) {
                    dep[v] = dep[u] + 1;
                    q.push(v);
                }
            }
        }
        return dep[T];
    }

    int dfs(int u, int flow) {
        if ((u == T) || (!flow)) return flow;

        int ret = 0;
        for (int& i = cur[u]; ~i; i = e[i].nxt) {
            int v = e[i].v, d;
            if ((dep[v] == dep[u] + 1) &&
                (d = dfs(v, min(flow - ret, e[i].cap - e[i].flow)))) {
                ret += d;
                e[i].flow += d;
                e[i ^ 1].flow -= d;
                if (ret == flow) return ret;
            }
        }
        return ret;
    }

    void dinic() {
        while (bfs()) {
            memcpy(cur, fir, sizeof(int) * (n + 1));
            maxflow += dfs(S, INF);
        }
    }
} mf;

int mp[300][300];

int n, m;
ll c,d;
void dis(){
    for(int i = 1; i <=n ;++i){
        for(int j = 1; j<=m;++j){
            cout<<mp[i][j]<<' ';
        }cout<<'\n';
    }cout<<'\n';
}

int ctx[300],cty[300];

int main() {

    cin >> n >> m>>c>>d;
    ll cst = d*n*m;
    for(int i=1;i<=n;++i){
        for(int j = 1; j<=m;++j){
            char c;cin>>c;
            if(c=='*'){
                mp[i][j] = 0;
                cst -=d;
            }else{
                mp[i][j] = 1;
            }
        }
    }

    for(int i = 1; i <= n ;++i){
        for(int j = 1; j <=m;++j){
            ctx[i]+=mp[i][j];
        }
    }

    for(int j = 1 ; j<=m;++j){
        for(int i = 1; i <= n ;++i){
            cty[j] += mp[i][j];
        }
    }

    ll ans =INF;
    ans = min(ans,cst);

    for(int k =  1;  k <= max(n,m);++k){
        cst+=c;

        mf.init(1);
        mf.n = n+m+2;
        mf.S = n+m+1;
        mf.T = n+m+2;


        for(int i = 1; i<= n ;++i)
            for(int j = 1; j<=m;++j)
                if(mp[i][j])
                    mf.addedge(i,n+j,1);

        for(int i = 1; i<=n;++i)mf.addedge(mf.S,i,min(k,ctx[i]));
        for(int j = 1; j<=m;++j)mf.addedge(n+j,mf.T,min(k,cty[j]));

        mf.dinic();
        ans = min(ans,cst - mf.maxflow*d);
    }

    cout<<ans;

    return 0;
}