#include <cstdio>
#include<math.h>
#include<queue>
using namespace std;
//这道题我的思路就是,怎么体现广度搜索的话,算是把第一种情况能分成的情况全都罗列出来
//这算是第一层,然后把第一层全都放到队列中,并且他们的步数是step,然后再开始判断从第一层出来的第二层,
//操作最难的还是在于怎么进行分层和记录步数,应该是需要一个结构体数组,他里面包含了一个数组,就是每个容量还有就是步骤
int a,b,c,k,mmin;//算是每个的容量吧
struct M{
int x,y,z;
int step;
}Node,Temp,mark;//这个结构体是用来进行分配的
//最大的难点还是在于怎么分配上面吧
bool f[100][100][100]={false};//初始化是没有访问过的
//再开始写BFS算法
int BFS(M x){
//创建一个队列,然后把这个点入栈
queue<M> q;
q.push(x);
while(!q.empty()){
//取出队头元素,再删除
Temp=q.front();
q.pop();
//然后开始判断这个位置是不是有结果
if(Temp.x==k||Temp.y==k||Temp.z==k)
return Temp.step;//直接返回结果
//如果不是的话,就需要开始浩大的分类 V1到v2 v1到v3 等等
//首先是v1到v2 也就是:现在v1的值大于0(也就是能拿出去) v2内的东西小于V2的容量
if(Temp.x>0&&b>Temp.y){
//那就是选择最小能移动的值,从v1的当前值和v2容量减去v2当前值之间进行比较
mmin=min(Temp.x,b-Temp.y);
//然后再用一个新结构体对他进行操作
mark.x=Temp.x-mmin;
mark.y=Temp.y+mmin;
mark.z=Temp.z;
mark.step=Temp.step+1;
//并且要对这个位置进行访问判定
if(!f[mark.x][mark.y][mark.z]){
//没有被访问过,那就设置为访问过了,并且进栈
f[mark.x][mark.y][mark.z]=true;
q.push(mark);
}
}
//然后是v1到v3
if(Temp.x>0&&c>Temp.z){
//那就是选择最小能移动的值,从v1的当前值和v2容量减去v2当前值之间进行比较
mmin=min(Temp.x,c-Temp.z);
//然后再用一个新结构体对他进行操作
mark.x=Temp.x-mmin;
mark.y=Temp.y;
mark.z=Temp.z+mmin;
mark.step=Temp.step+1;
//并且要对这个位置进行访问判定
if(!f[mark.x][mark.y][mark.z]){
//没有被访问过,那就设置为访问过了,并且进栈
f[mark.x][mark.y][mark.z]=true;
q.push(mark);
}
}
//然后是v2到v3
if(Temp.y>0&&c>Temp.z){
//那就是选择最小能移动的值,从v1的当前值和v2容量减去v2当前值之间进行比较
mmin=min(Temp.y,c-Temp.z);
//然后再用一个新结构体对他进行操作
mark.x=Temp.x;
mark.y=Temp.y-mmin;
mark.z=Temp.z+mmin;
mark.step=Temp.step+1;
//并且要对这个位置进行访问判定
if(!f[mark.x][mark.y][mark.z]){
//没有被访问过,那就设置为访问过了,并且进栈
f[mark.x][mark.y][mark.z]=true;
q.push(mark);
}
}
//然后是v2到v1
if(Temp.y>0&&a>Temp.x){
//那就是选择最小能移动的值,从v1的当前值和v2容量减去v2当前值之间进行比较
mmin=min(Temp.y,a-Temp.x);
//然后再用一个新结构体对他进行操作
mark.x=Temp.x+mmin;
mark.y=Temp.y-mmin;
mark.z=Temp.z;
mark.step=Temp.step+1;
//并且要对这个位置进行访问判定
if(!f[mark.x][mark.y][mark.z]){
//没有被访问过,那就设置为访问过了,并且进栈
f[mark.x][mark.y][mark.z]=true;
q.push(mark);
}
}
//v3到v1
if(Temp.z>0&&a>Temp.x){
//那就是选择最小能移动的值,从v1的当前值和v2容量减去v2当前值之间进行比较
mmin=min(Temp.z,a-Temp.x);
//然后再用一个新结构体对他进行操作
mark.x=Temp.x+mmin;
mark.y=Temp.y;
mark.z=Temp.z-mmin;
mark.step=Temp.step+1;
//并且要对这个位置进行访问判定
if(!f[mark.x][mark.y][mark.z]){
//没有被访问过,那就设置为访问过了,并且进栈
f[mark.x][mark.y][mark.z]=true;
q.push(mark);
}
}
//最后是v3到v2
if(Temp.z>0&&b>Temp.y){
//那就是选择最小能移动的值,从v1的当前值和v2容量减去v2当前值之间进行比较
mmin=min(Temp.z,b-Temp.y);
//然后再用一个新结构体对他进行操作
mark.x=Temp.x;
mark.y=Temp.y+mmin;
mark.z=Temp.z-mmin;
mark.step=Temp.step+1;
//并且要对这个位置进行访问判定
if(!f[mark.x][mark.y][mark.z]){
//没有被访问过,那就设置为访问过了,并且进栈
f[mark.x][mark.y][mark.z]=true;
q.push(mark);
}
}
}
return -1;
}
int main(){
//因为是多次输入的
while(scanf("%d%d%d%d",&a,&b,&c,&k)!=EOF){
//对第一个位置进行初始化
Node.x=a;
Node.y=0;
Node.z=0;
Node.step=0;
f[a][b][c]=true;//先标记为访问过了
int step=BFS(Node);
if(step>0){
printf("YES\n");
printf("%d\n",step);
}
else
printf("No\n");
}
}