QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#184963#5114. Cells Coloringzwp1234TL 1ms3540kbC++174.1kb2023-09-21 15:04:162023-09-21 15:04:16

Judging History

你现在查看的是最新测评结果

  • [2023-09-21 15:04:16]
  • 评测
  • 测评结果:TL
  • 用时:1ms
  • 内存:3540kb
  • [2023-09-21 15:04:16]
  • 提交

answer

//# pragma GCC optimize(2)
#include <bits/stdc++.h>
//#define int long long
#define SZ(x) (int)(x).size()
#define fs first
#define sc second
#define PII pair<int,int>
#define ls(u) (u)<<1
#define rs(u) (u)<<1|1
#define pb push_back
#define eb emplace_back
#define vi vector<int>
#define vvi vector<vi >
#define vI vector<PII>
#define db double
#define all(a) (a).begin(),(a).end()
#define die cout << "???" << endl
//#define endl '\n'
using namespace std;
const int mod = 998244353;
//const int mod = 1169996969;
//const int mod = 1e9+7;
//const int inf = LLONG_MAX;
typedef long long ll;
const int inf = 0x3f3f3f3f;
const ll INF = 0x3f3f3f3f3f3f3f3f;
//const int inf = 1e9;
//const db PI = acos(-1.0);
const db eps = 1e-5;
typedef unsigned long long  ull;
const int N = 3e2+5,M = 2e5+5,K = 605;

inline int rd() {
    int f = 0; int x = 0; char ch = getchar();
    for (; !isdigit(ch); ch = getchar()) f |= (ch == '-');
    for (; isdigit(ch); ch = getchar()) x = (x << 1) + (x << 3) + ch - '0';
    if (f) x = -x;
    return x;
}

void write(int a) {if(a>=10)write(a/10);putchar(a%10+'0');}
void wt(int a) {if(a < 0) {putchar('-');a = -a;}write(a);}
void wwt(int a){wt(a);putchar('\n');}

// void write(ll a) {if(a>=10)write(a/10);putchar(a%10+'0');}
// void wt(ll a) {if(a < 0) {putchar('-');a = -a;}write(a);}
// void wwt(ll a){wt(a);putchar('\n');}

int lowbit(int x) {return (x&(-x));}
int n,m,k;
int S,T;

char s[N][N];
int col[N],row[N];

int h[K],e[M],ne[M],f[M],tot;

void add(int u,int v, int x)
{
    e[tot] = v,ne[tot] = h[u],f[tot] = x,h[u] = tot++;
    e[tot] = u,ne[tot] = h[v],f[tot] = 0,h[v] = tot++;
}

int dep[K],cur[K],gap[K],q[K];
int nod;

void bfs()
{
    int hd = 0, tl = 0;
    q[++tl] = T, ++gap[dep[T] = 1];

    while(hd < tl)
    {
        int u = q[++hd];

        for(int i = h[u];~i;i = ne[i])
        {
            int v = e[i];
            if(!dep[v] && f[i^1])
            {
                dep[v] = dep[u] + 1;
                gap[dep[v]]++;
                q[++tl] = v;
            }
        }
    }
}

int dfs(int u,int a)
{
    if(u == T || a == 0) return a;

    int flow = 0,tmp;

    for(int& i = cur[u];~i;i = ne[i])
    {
        int v = e[i];
        if(dep[v] == dep[u] - 1 && (tmp = dfs(v,min(a,f[i]))))
        {
            flow += tmp;
            f[i] -= tmp;
            f[i^1] += tmp;
            a -= tmp;

            if(a == 0)
            {
                return flow;
            }
        }
    }

    --gap[dep[u]];
    if(!gap[dep[u]]) gap[S] = nod + 1;
    ++gap[++dep[u]];

    return flow;
}

int ISAP()
{
    static int res = 0;
    bfs();
    while(dep[S] <= nod)
    {
        memcpy(cur,h,sizeof(*cur)*nod);
        res += dfs(S,inf);
    }

    return res;
}

signed main()
{
    ios::sync_with_stdio(0);
    cin.tie(0),cout.tie(0);
    #ifndef ONLINE_JUDGE
    freopen("in.txt", "r+", stdin);
    freopen("out.txt", "w+", stdout);
    #endif

    int c,d;
    
    cin >> n >> m >> c >> d;

    int mxk = 0;
    int sum = 0;

    S = 0;
    T = m + n + 1;
    nod = T + 1;

    for(int i = S;i<=T;i++)
    {
        h[i] = -1;
    }

    for(int i = 1;i<=n;i++)
    {
        cin >> (s[i] + 1);
        for(int j = 1;j<=m;j++) s[i][j] = s[i][j] == '.',row[i] += s[i][j],col[j] += s[i][j],mxk = max(mxk,col[j]),sum += s[i][j];
        mxk = max(mxk,row[i]);
    }

    for(int i = 1;i<=n;i++) add(S,i,0);
    for(int i = 1;i<=m;i++) add(i + n ,T,0);

    for(int i = 1;i<=n;i++)
    {
        for(int j = 1;j<=m;j++)
        {
            if(s[i][j])
            {
                add(i,j+n,1);
            }
        }
    }


    ll ans = INF;  

    for(int k = 0;k<=mxk;k++)
    {
        if(k)
        {
            memset(gap,0,sizeof (*gap)*nod);
            memset(dep,0,sizeof (*dep)*nod);
            for(int i = 0;i<n;i++) ++f[i*2];
            for(int i = 0;i<m;i++) ++ f[(i+n)*2];
        }

        ans = min(ans,1ll*c*k + 1ll*d*(sum - ISAP()));
    }

    cout << ans;

    return 0;

}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 1ms
memory: 3484kb

input:

3 4 2 1
.***
*..*
**..

output:

4

result:

ok 1 number(s): "4"

Test #2:

score: 0
Accepted
time: 0ms
memory: 3540kb

input:

3 4 1 2
.***
*..*
**..

output:

2

result:

ok 1 number(s): "2"

Test #3:

score: -100
Time Limit Exceeded

input:

250 250 965680874 9042302
..**.*****..**..**.*****..***..***.**......*.***.*...***.*....*.**.*.**.*.*.****...*.******.***.************....**.*..*..***.*******.*.***.*..**..****.**.*.*..***.****..**.....***.....*.****...*...*.***..****..**.*.*******..*.*******.*.*.*.****.*.***
....**.*******.*.******...

output:


result: