Measurement of the Y(1S) pair production cross section and search for resonances decaying to Y(1S)μ+μ− in proton-proton collisions at √s=13TeV