Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #include <fstream>
- #include <iostream>
- #define LL long long
- using namespace std;
- /// ecuatia: a1*p4+a2*p3+a3*p2+a4*p1=k*x => a1*p4+a2*p3=k*x-(a3*p2+a4*p1)
- /// https://www.youtube.com/watch?v=ws5dDS1ByNk
- LL div(LL divizor,LL n)
- {
- if(n<0)
- n=-n;
- LL d;
- if(n==1)
- return n;
- else
- for(d=divizor;d<=n;++d)
- if(n%d==0)
- break;
- return d;
- }
- void term_a4(LL d,LL &a1,LL &a2,LL &a3,LL &a4)
- {
- LL A1=a1,A2=a2,A3=a3,A4=a4;
- A1=d+A1;A2=A1*d+A2;A3=A2*d+A3;A4=A3*d+A4;
- cout<<A1<<' '<<A2<<' '<<A3<<' '<<A4<<'\n';
- if(A4==0)
- a1=A1,a2=A2,a3=A3,a4=A4;
- }
- void term_a3(LL d,LL &a1,LL &a2,LL &a3)
- {
- LL A1=a1,A2=a2,A3=a3;
- A1=d+A1;A2=A1*d+A2;A3=A2*d+A3;
- cout<<A1<<' '<<A2<<' '<<A3<<'\n';
- if(A3==0)
- a1=A1,a2=A2,a3=A3;
- }
- int main()
- {
- int n,b,x,k;
- LL a1,a2,a3,a4;
- ifstream f("recc.in");
- f>>n>>b>>x>>k>>a1>>a2>>a3>>a4;
- f.close();
- LL t,divizor=2; /// t^4-a1*t^3-a2*t^2-a3*t-a4=0
- a1=-a1;a2=-a2;a3=-a3;a4=-a4;
- if(a4==1)
- divizor=1;
- else
- divizor=2;
- while(a4)
- {
- while(a4%divizor)
- ++divizor;
- cout<<"DIVIZ= "<<divizor<<endl,term_a4(divizor,a1,a2,a3,a4),++divizor;
- }
- if(a3==1)
- divizor=1;
- else
- divizor=2;
- while(a3)
- {
- while(a3%divizor)
- ++divizor;
- cout<<"DIVIZ= "<<divizor<<endl,term_a3(divizor,a1,a2,a3),++divizor;
- }
- int rez=0;
- ofstream g("recc.out");
- g<<rez<<'\n';
- g.close();
- return 0;
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement