ちょっと考えたが、先人達のコードも有り、
案外　あっさりと解けた。

ここから、縮めてくる連中は出てくるだろうか。。。


【初稿】: 91 バイト

x#y|x^2>=y='*'|1>0=' '
main=do n<-readLn;mapM putStrLn[[n#(x^2+y^2)|x<-[-n..n]]|y<-[-n..n]]

【修正後1】: 88 バイト

main=do n<-readLn
        mapM putStrLn[["* "!!(0^(max(1+n^2-x^2-y^2)0))|x<-[-n..n]]|y<-[-n..n]]

(※)実際は、readLn と mapM は ; でつないでます。

【その他】: 90 バイト

main=interact$\\
(\n->unlines$[["* "!!(0^(max(1+n^2-x^2-y^2)0))|x<-[-n..n]]|y<-[-n..n]]).read

(※)\\ は、行をつないでいる。

【修正後2】: 84 バイト

main=do n<-readLn
        mapM putStrLn[[" *"!!(0^max(x^2+y^2-n^2)0)|x<-[-n..n]]|y<-[-n..n]]

(※)実際は、readLn と mapM は ; でつないでます。

【修正後3】: 82 バイト

main=do n<-readLn
        mapM putStrLn[["* "!!div(x^2+y^2)(n^2+1)|x<-[-n..n]]|y<-[-n..n]]

(※)実際は、readLn と mapM は ; でつないでます。