l = qgis.utils.iface.activeLayer()
p = l.dataProvider()
for f in p.getFeatures(QgsFeatureRequest()):
    g = f.geometry()
    if len(g.asPolygon()) > 0 and len(g.asPolygon()[0]) == 4:
        l.select(f.id())