#include <bits/stdc++.h>
using namespace std;

#define rep(i,n) for(int i=1;i<=n;++i)
#define all(x) begin(x), end(x)
#define sz(x) (int)(x).size()
#define pb push_back
typedef pair<int, int> pii;
typedef vector<int> vi;
typedef long double ld;
typedef long long ll;
const ld eps=1e-12;
struct P{
    ld x,y;
    P operator+(P p){return (P){x+p.x,y+p.y};}
    P operator*(ld c){return (P){x*c,y*c};}
    P operator-(P p){return (P){x-p.x,y-p.y};}
    ld dis(){return hypot(x,y);}
};
int n,m;
ld d;
vector<P> ans;
void sol(){
    P p=(P){0,0},T=(P){1.L*n,1.L*m};
    for(int i=0;i<=n;++i){
        int stp;
        if(i&1)stp=-1;else stp=1;
        for(int j=0;j<=1;++j){
            P nxt=p;
            if(j==0)nxt.y+=stp*m;
            else nxt.x+=1;
            if((p-nxt).dis()+(nxt-T).dis()>d+eps){
                ans.pb(p);
                ld l=0,r=1;
                rep(tim,256){
                    ld mid=(l+r)/2;
                    P pos=p*(1-mid)+nxt*mid;
                    if((p-pos).dis()+(pos-T).dis()>d)r=mid;else l=mid;
                }
                P pos=p*(1-l)+nxt*l;
                if((p-pos).dis()>eps)ans.pb(pos);
                ans.pb(T);
                return;
            }else{
                ans.pb(p);
                d-=(p-nxt).dis();
                p=nxt;

            }
        }
    }
}
int main() {
    scanf("%d%d%Lf",&n,&m,&d);
    sol();
    printf("%d\n",ans.size());
    for(P p:ans)printf("%.15Lf %.15Lf\n",p.x,p.y);
    return 0;
}