请用PASCAL编程:如何判定一个三位数为水仙花数。
program lianxi06(input,output);    Var  n,a,b,c:integer;;    Begin     Write('n=');readln(n);        a:=n div 100;        b:=(n div 10) -a*10 mod 10;        c:=n -a*100-b*10; if n=c*c*c+b*b*b+a*a*a then wirteln(‘yes’) else writeln(‘no’); End.