#include<bits/stdc++.h>
#define st first
#define nd second
#define pb push_back
#define ppb pop_back
#define ii pair<int,int>
#define ll long long
#define orta ((bas+son)>>1)
#define sz(x) ((int)x.size())
#define all(x) x.begin(),x.end()
#define inf 1000000000
#define N 2019
#define MOD 998244353
using namespace std;
#define greater bos
struct frac {
ll u,d;
};
frac mul(frac a,frac b) {
ll u=(ll)a.u*b.u;
ll d=(ll)a.d*b.d;
ll g=__gcd(abs(u),abs(d));
return {u/g,d/g};
}
frac dvd(frac a,frac b) {
swap(b.u,b.d);
return mul(a,b);
}
frac sub(frac a,frac b) {
ll u=(ll)a.u*b.d-(ll)b.u*a.d;
ll d=(ll)a.d*b.d;
ll g=__gcd(abs(u),abs(d));
return {u/g,d/g};
}
frac add(frac a,frac b) {
ll u=(ll)a.u*b.d+(ll)b.u*a.d;
ll d=(ll)a.d*b.d;
ll g=__gcd(abs(u),abs(d));
return {u/g,d/g};
}
bool eq(frac a,frac b) {
return (ll)a.u*b.d==(ll)b.u*a.d;
}
bool greateq(frac a,frac b) {
return (ll)a.u*b.d>=(ll)b.u*a.d;
}
bool greater(frac a,frac b) {
return (ll)a.u*b.d>(ll)b.u*a.d;
}
frac gmin(frac a,frac b) {
return greater(a,b)?b:a;
}
int gnext(frac a) {
return a.u/a.d+1;
}
int n,l;
int v[N][N],u[N],ptr[N];
frac req[N];
vector<frac> cp[N];
int main() {
scanf("%d %d",&n,&l);
for(int i=1;i<=n;i++) {
req[i]={0,n};
for(int j=1;j<=l;j++) {
scanf("%d",&v[i][j]);
req[i].u+=v[i][j];
}
}
for(int i=1;i<=n;i++) {
//cerr<<i<<"\n";
frac cur={0,1};
while(sz(cp[i])<n) {
//cerr<<cur.u<<" "<<cur.d<<"\n";
cp[i].pb(cur);
frac rq=req[i];
for(int j=gnext(cur);j<=l;j++) {
//cerr<<j<<" "<<cur.u<<" "<<cur.d<<"\n";
frac df=sub({j,1},cur);
frac sum=mul(df,{v[i][j],1});
if(greateq(sum,rq)) {
frac go=dvd(rq,{v[i][j],1});
cur=add(cur,go);
break ;
}
else {
rq=sub(rq,sum);
cur={j,1};
}
}
}
cp[i].pb({l,1});
}
for(int i=1;i<=n;i++) {
//cerr<<"--------------------------\n";
//for(auto x:cp[i]) printf("%d %d\n",x.u,x.d);
}
//cerr<<"--------------------------\n";
//cerr<<"OK\n";
vector<frac> res;
vector<int> ans;
for(int i=1;i<=n;i++) {
frac mn={inf,1};
int tut=-1;
for(int j=1;j<=n;j++) {
if(u[j]) continue ;
while(sz(res) && ptr[j]<n && greater(res.back(),cp[j][ptr[j]])) {
++ptr[j];
}
if(ptr[j]<n) {
mn=gmin(mn,cp[j][ptr[j]+1]);
if(eq(mn,cp[j][ptr[j]+1])) {
tut=j;
}
}
}
ans.pb(tut);
u[tut]=1;
res.pb(mn);
}
res.ppb();
for(auto x:res) printf("%lld %lld\n",x.u,x.d);
for(auto x:ans) printf("%d ",x);
}
Compilation message
naan.cpp: In function 'int main()':
naan.cpp:102:7: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
scanf("%d %d",&n,&l);
~~~~~^~~~~~~~~~~~~~~
naan.cpp:110:9: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
scanf("%d",&v[i][j]);
~~~~~^~~~~~~~~~~~~~~
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
2 ms |
376 KB |
Output is correct |
2 |
Correct |
3 ms |
376 KB |
Output is correct |
3 |
Correct |
3 ms |
376 KB |
Output is correct |
4 |
Correct |
3 ms |
376 KB |
Output is correct |
5 |
Correct |
3 ms |
376 KB |
Output is correct |
6 |
Correct |
2 ms |
248 KB |
Output is correct |
7 |
Correct |
3 ms |
376 KB |
Output is correct |
8 |
Correct |
3 ms |
376 KB |
Output is correct |
9 |
Correct |
2 ms |
376 KB |
Output is correct |
10 |
Correct |
3 ms |
376 KB |
Output is correct |
11 |
Correct |
3 ms |
376 KB |
Output is correct |
12 |
Correct |
3 ms |
376 KB |
Output is correct |
13 |
Correct |
3 ms |
376 KB |
Output is correct |
14 |
Correct |
3 ms |
376 KB |
Output is correct |
15 |
Correct |
3 ms |
376 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
3 ms |
376 KB |
Output is correct |
2 |
Correct |
3 ms |
376 KB |
Output is correct |
3 |
Correct |
3 ms |
376 KB |
Output is correct |
4 |
Correct |
4 ms |
376 KB |
Output is correct |
5 |
Correct |
4 ms |
504 KB |
Output is correct |
6 |
Correct |
4 ms |
376 KB |
Output is correct |
7 |
Correct |
3 ms |
376 KB |
Output is correct |
8 |
Correct |
3 ms |
376 KB |
Output is correct |
9 |
Correct |
4 ms |
376 KB |
Output is correct |
10 |
Correct |
4 ms |
376 KB |
Output is correct |
11 |
Correct |
4 ms |
376 KB |
Output is correct |
12 |
Correct |
2 ms |
376 KB |
Output is correct |
13 |
Correct |
3 ms |
376 KB |
Output is correct |
14 |
Correct |
4 ms |
424 KB |
Output is correct |
15 |
Correct |
4 ms |
376 KB |
Output is correct |
16 |
Correct |
5 ms |
504 KB |
Output is correct |
17 |
Correct |
6 ms |
376 KB |
Output is correct |
18 |
Correct |
5 ms |
380 KB |
Output is correct |
19 |
Correct |
5 ms |
504 KB |
Output is correct |
20 |
Correct |
4 ms |
376 KB |
Output is correct |
21 |
Correct |
5 ms |
504 KB |
Output is correct |
22 |
Correct |
5 ms |
504 KB |
Output is correct |
23 |
Correct |
2 ms |
376 KB |
Output is correct |
24 |
Correct |
4 ms |
376 KB |
Output is correct |
25 |
Correct |
3 ms |
376 KB |
Output is correct |
26 |
Correct |
3 ms |
376 KB |
Output is correct |
27 |
Correct |
4 ms |
504 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
2 ms |
376 KB |
Output is correct |
2 |
Correct |
3 ms |
376 KB |
Output is correct |
3 |
Correct |
3 ms |
376 KB |
Output is correct |
4 |
Correct |
3 ms |
376 KB |
Output is correct |
5 |
Correct |
3 ms |
376 KB |
Output is correct |
6 |
Correct |
2 ms |
248 KB |
Output is correct |
7 |
Correct |
3 ms |
376 KB |
Output is correct |
8 |
Correct |
3 ms |
376 KB |
Output is correct |
9 |
Correct |
2 ms |
376 KB |
Output is correct |
10 |
Correct |
3 ms |
376 KB |
Output is correct |
11 |
Correct |
3 ms |
376 KB |
Output is correct |
12 |
Correct |
3 ms |
376 KB |
Output is correct |
13 |
Correct |
3 ms |
376 KB |
Output is correct |
14 |
Correct |
3 ms |
376 KB |
Output is correct |
15 |
Correct |
3 ms |
376 KB |
Output is correct |
16 |
Correct |
3 ms |
376 KB |
Output is correct |
17 |
Correct |
3 ms |
376 KB |
Output is correct |
18 |
Correct |
3 ms |
376 KB |
Output is correct |
19 |
Correct |
4 ms |
376 KB |
Output is correct |
20 |
Correct |
4 ms |
504 KB |
Output is correct |
21 |
Correct |
4 ms |
376 KB |
Output is correct |
22 |
Correct |
3 ms |
376 KB |
Output is correct |
23 |
Correct |
3 ms |
376 KB |
Output is correct |
24 |
Correct |
4 ms |
376 KB |
Output is correct |
25 |
Correct |
4 ms |
376 KB |
Output is correct |
26 |
Correct |
4 ms |
376 KB |
Output is correct |
27 |
Correct |
2 ms |
376 KB |
Output is correct |
28 |
Correct |
3 ms |
376 KB |
Output is correct |
29 |
Correct |
4 ms |
424 KB |
Output is correct |
30 |
Correct |
4 ms |
376 KB |
Output is correct |
31 |
Correct |
5 ms |
504 KB |
Output is correct |
32 |
Correct |
6 ms |
376 KB |
Output is correct |
33 |
Correct |
5 ms |
380 KB |
Output is correct |
34 |
Correct |
5 ms |
504 KB |
Output is correct |
35 |
Correct |
4 ms |
376 KB |
Output is correct |
36 |
Correct |
5 ms |
504 KB |
Output is correct |
37 |
Correct |
5 ms |
504 KB |
Output is correct |
38 |
Correct |
2 ms |
376 KB |
Output is correct |
39 |
Correct |
4 ms |
376 KB |
Output is correct |
40 |
Correct |
3 ms |
376 KB |
Output is correct |
41 |
Correct |
3 ms |
376 KB |
Output is correct |
42 |
Correct |
4 ms |
504 KB |
Output is correct |
43 |
Correct |
474 ms |
11000 KB |
Output is correct |
44 |
Incorrect |
3032 ms |
46352 KB |
X_i is not increasing |
45 |
Halted |
0 ms |
0 KB |
- |