QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#105216 | #6395. Equation Discovering | Sommohito# | WA | 2ms | 3728kb | C++20 | 4.2kb | 2023-05-13 17:07:36 | 2023-05-13 17:07:40 |
Judging History
answer
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef long double ld;
const int N = 2e5 +5;
typedef pair<int,int> pii;
//int cnt(int comp)
//{
// ll ans=0;
//
// if(comp-2>=0)
// {
// int baki=comp-2;
// for(int l=0;l<=baki;l++)
// {
// int L=cnt(l);
// int R=cnt(baki-l);
// ans+=L*R*4;
// }
// }
// if(comp-1>=0)
// {
// ans+=cnt(comp-1)*2;
// }
// ans++;
//
// return ans;
//}
class S;
class S
{
public:
string op;
S * lft;
S * right;
ld eval(ld x)
{
if(op=="+") return lft->eval(x)+right->eval(x);
if(op=="-") return lft->eval(x)-right->eval(x);
if(op=="*") return lft->eval(x)*right->eval(x);
if(op=="/")
{
if(right->eval(x)>0.01) return lft->eval(x)/right->eval(x);
else return -1e500;
}
if(op=="sin") return sin(lft->eval(x));
if(op=="cos") return cos(lft->eval(x));
if(op=="x") return x;
}
string make()
{
if(op=="+") return "("+lft->make()+"+"+right->make()+")";
if(op=="-") return "("+lft->make()+"-"+right->make()+")";
if(op=="*") return "("+lft->make()+"*"+right->make()+")";
if(op=="/") return "("+lft->make()+"/"+right->make()+")";
if(op=="sin") return "sin("+lft->make()+")";
if(op=="cos") return "cos("+lft->make()+")";
if(op=="x") return "x";
}
};
vector<ld> x,y;
vector<S*> save[10];
int dp[10];
int n;
vector<S*> cnt(int comp)
{
if(dp[comp])
return save[comp];
dp[comp]=true;
vector<S*> ans;
if(comp-2>=0)
{
int baki=comp-2;
for(int l=0; l<=baki; l++)
{
auto L=cnt(l);
auto R=cnt(baki-l);
for(auto &l:L)
{
for(auto &r:R)
{
{
S * notun = new S();
notun->op="+";
notun->lft=l;
notun->right=r;
ans.push_back(notun);
}
{
S * notun = new S();
notun->op="-";
notun->lft=l;
notun->right=r;
ans.push_back(notun);
}
{
S * notun = new S();
notun->op="*";
notun->lft=l;
notun->right=r;
ans.push_back(notun);
}
{
S * notun = new S();
notun->op="/";
notun->lft=l;
notun->right=r;
ans.push_back(notun);
}
}
}
}
}
if(comp-1>=0)
{
auto L=cnt(comp-1);
for(auto &l:L)
{
{
S * notun = new S();
notun->op="sin";
notun->lft=l;
ans.push_back(notun);
}
{
S * notun = new S();
notun->op="cos";
notun->lft=l;
ans.push_back(notun);
}
}
}
S * notun = new S();
notun->op="x";
ans.push_back(notun);
return save[comp]=ans;
}
void TEST_CASES()
{
cin>>n;
// n=2;
x.resize(n);
y.resize(n);
for(int i=0; i<n; i++)
{
cin>>x[i]>>y[i];
// x[i]=1e3-(i+1);
// y[i]=sin(x[i])*cos(x[i]);
}
for(int comp=0;comp<=9;comp++)
{
auto got=cnt(comp);
for(auto r:got)
{
bool f=1;
for(int i=0;i<n;i++)
{
if(abs(r->eval(x[i])-y[i])>1e-3) f=0;
}
if(f)
{
cout<<r->make();
exit(0);
}
}
}
}
/*
*/
int32_t main()
{
#ifndef APURBA
ios_base::sync_with_stdio(false);
cin.tie(nullptr);
#endif
//freopen("input->txt","r",stdin);
//freopen("out1->txt","w",stdout);
int t=1;
//cin>>t;
while(t--)
{
TEST_CASES();
}
return 0;
}
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 2ms
memory: 3516kb
input:
3 1.000000 1.000000 2.000000 4.000000 3.000000 9.000000
output:
(x*x)
result:
ok great!!
Test #2:
score: 0
Accepted
time: 2ms
memory: 3728kb
input:
3 0.618000 1.517072 0.314000 3.132637 1.414000 0.494016
output:
(sin(x)/(x*x))
result:
ok great!!
Test #3:
score: 0
Accepted
time: 2ms
memory: 3612kb
input:
5 77.685777 233.057331 -66.445083 -199.335249 79.966717 239.900151 84.982130 254.946390 -31.528900 -94.586700
output:
(x+(x+x))
result:
ok great!!
Test #4:
score: -100
Wrong Answer
time: 2ms
memory: 3648kb
input:
5 25.032427 -0.100652 38.727324 1.658518 27.684334 -0.669555 64.282391 8.275303 52.640700 -0.962660
output:
sin((x/(x-x)))
result:
wrong answer div a small number