QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#184963 | #5114. Cells Coloring | zwp1234 | TL | 1ms | 3540kb | C++17 | 4.1kb | 2023-09-21 15:04:16 | 2023-09-21 15:04:16 |
Judging History
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;
}
詳細信息
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 ..**.*****..**..**.*****..***..***.**......*.***.*...***.*....*.**.*.**.*.*.****...*.******.***.************....**.*..*..***.*******.*.***.*..**..****.**.*.*..***.****..**.....***.....*.****...*...*.***..****..**.*.*******..*.*******.*.*.*.****.*.*** ....**.*******.*.******...