| ||||||||||
| Online Judge | Problem Set | Authors | Online Contests | User | ||||||
|---|---|---|---|---|---|---|---|---|---|---|
| Web Board Home Page F.A.Qs Statistical Charts | Current Contest Past Contests Scheduled Contests Award Contest | |||||||||
为什么会一直超时?????????#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <cstdio>
#include <cstring>
using namespace std;
const int INF = 1 << 15;
const int MAXN = 200 + 25;
struct Edge{
int from,to,cap,flow,cost;
Edge(){};
Edge(int _from,int _to,int _cap,int _flow)
:from(_from),to(_to),cap(_cap),flow(_flow){};
};
vector<Edge> edges;
vector<int> G[MAXN];
int cur[MAXN],d[MAXN];
bool vst[MAXN];
int N,M,src,sink;
int tt,ss;
int in[MAXN];
int B[MAXN][MAXN],C[MAXN][MAXN];
int ans[MAXN][MAXN];
bool flag;
void init(){
ss = N + M + 1; tt = ss + 1;
src = tt + 1; sink = src + 1;
for(int i = 0;i <= sink;++i)
G[i].clear();
edges.clear();
}
inline void addEdge(int from,int to,int cap){
edges.push_back(Edge(from,to,cap,0));
edges.push_back(Edge(to,from,0,0));
int sz = edges.size();
G[from].push_back(sz - 2);
G[to].push_back(sz - 1);
}
inline void input(int x,int y,char *ope,int c){
int x1,x2,y1,y2;
x1 = x2 = x;
y1 = y2 = y;
if(!x) x1 = 1,x2 = N;
if(!y) y1 = 1,y2 = M;
for(int i = x1;i <= x2;++i){
for(int j = y1;j <= y2;++j){
if(ope[0] == '=')
B[i][j] = C[i][j] = c;
else if(ope[0] == '>')
B[i][j] = max(B[i][j],c + 1);
else
C[i][j] = min(C[i][j],c - 1);
if(C[i][j] < B[i][j])
flag = false;
}
}
}
bool BFS(){
memset(vst,0,sizeof(vst));
queue<int> Q;
Q.push(src);
d[src] = 0;
vst[src] = 1;
while(!Q.empty()){
int x = Q.front(); Q.pop();
for(int i = 0;i < (int)G[x].size();++i){
Edge& e = edges[G[x][i]];
if(!vst[e.to] && e.cap > e.flow){
d[e.to] = d[x] + 1;
vst[e.to] = 1;
Q.push(e.to);
}
}
}
return vst[sink];
}
int DFS(int x,int a){
if(x == sink||a == 0)
return a;
int flow = 0,f;
for(int& i = cur[x];i < (int)G[x].size();++i){
Edge& e = edges[G[x][i]];
if(d[e.to] == d[x] + 1&&(f = DFS(e.to,min(a,e.cap - e.flow))) > 0){
e.flow += f;
edges[G[x][i]^1].flow -= f;
flow += f;
a -= f;
if(a == 0) break;
}
}
return flow;
}
int maxFlow(){
int flow = 0;
while(BFS()){
memset(cur,0,sizeof(cur));
flow += DFS(src,INF);
}
return flow;
}
int main()
{
// freopen("Input.txt","r",stdin);
int T;
scanf("%d",&T);
while(T--){
char ope[3];
int x,y,c;
scanf("%d%d",&N,&M);
init();
memset(in,0,sizeof(in));
for(int i = 1;i <= N;++i){
for(int j = 1;j <= M;++j){
B[i][j] = 0;
C[i][j] = INF;
}
}
for(int i = 1;i <= N;++i){ //行
scanf("%d",&c);
addEdge(ss,i,0);
in[ss] -= c;
in[i] += c;
}
for(int j = 1;j <= M;++j){ //列
scanf("%d",&c);
addEdge(j+N,tt,0);
in[j + N] -= c;
in[tt] += c;
}
int cas;
scanf("%d",&cas);
while(cas--){
scanf("%d%d%s%d",&x,&y,ope,&c);
input(x,y,ope,c);
}
if(flag){
puts("IMPOSSIBLE\n");
continue;
}
//建图
for(int i = 1;i <= N;++i){
for(int j = 1;j <= M;++j){
addEdge(i,j + N,C[i][j] - B[i][j]);
in[i] -= B[i][j];
in[j+N] += B[i][j];
}
}
int sum = 0;
for(int i = 1;i <= tt;++i){
if(in[i] > 0){
sum += in[i];
addEdge(src,i,in[i]);
}
if(in[i] < 0){
addEdge(i,sink,-in[i]);
}
}
addEdge(tt,ss,INF); int flow = maxFlow();
if(flow != sum){
puts("IMPOSSIBLE");
} else {
maxFlow();
for(int i = 1;i <= N;++i){
for(int j = 0;j < (int)G[i].size();++j){
Edge& e = edges[G[i][j]];
if(N < e.to && e.to < ss)
ans[e.from][e.to - N] = e.flow + B[e.from][e.to - N];
}
}
for(int i = 1;i <= N;++i){
for(int j = 1;j <= M;++j){
printf("%d ",ans[i][j]);
}
puts("");
}
}
puts("");
}
return 0;
}
Followed by: Post your reply here: |
All Rights Reserved 2003-2013 Ying Fuchen,Xu Pengcheng,Xie Di
Any problem, Please Contact Administrator